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FOREWORD 


The  work  described  in  this  Annual  Technical  Report  was  performed  at  the 
United  Technologies  Research  Center  (UTRC)  under  Contract  N00014-80-C-0476 , 
Modification  P0Q002,  entitled  "Study  of  Flow  Distribution  Control  Characteristics 
in  Marine  Gas  Turbine  Waste-Heat  Recovery  Systems",  for  the  Office  of  Naval 
Research  (ONR),  This  report  summarizes  results  obtained  for  the  Phase  II  (second 
year)  study  on  flow  distribution  control  characteristics  in  waste-heat  steam 
generators  which  was  preceded  by  the  first-year  study  on  diffusers.  Dr.  Simion  C. 
Kuo  is  the  Principal  Investigator  for  this  contract  program,  and  Dr.  Ho-Tien  Shu 
is  the  major  contributor  to  this  phase  of  the  study.  The  computer  program  used  in 
analyzing  the  steam-generator  was  derived  from  an  existing  Fuel  Vaporization  Model 
originally  developed  by  Messrs.  Chiappetta  and  Szetela,  both  of  UTRC. 

The  research  contract  was  signed  by  ONR  on  July  23,  1980,  and  the  Scientific 
Officer  is  Mr.  M.  Keith  Ell ingsworth.  Mechanics  Division,  ONR,  Arlington,  Virginia. 
Valuable  guidance  and  comments  received  from  Mr.  Ellingsworth  are  gratefully 
apprec iated. 
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SUMMARY 


The  objective  of  this  study  was  to  investigate  the  effect  of  flow 
distribution  control  on  the  design  and  performance  of  marine  gas  turbine 
waste-heat  steam  generators.  The  applicable  steam  generator  design  concepts 
and  general  design  consideration  were  reviewed  and  critical  problems  associated 
with  the  design  of  marine  waste-heat  steam  generators  were  identified.  A 
once-t hrough  counter  crossflow  heat  exchanger  was  selected  as  the  candidate 
waste-heat  steam  generator  for  recovering  the  waste  heat  from  the  exhaust  of  a 
marine  gas  turbine.  A  two-dimensional  heat  exchanger  model  suitable  for  the 
study  objective  was  formulated  and  computerized.  Parametric  performance 
analyses  were  made  of  the  waste-heat  steam  generators  for  four  different 
tube  arrangements  from  which  the  most  desirable  design  was  selected  (as  base¬ 
line  waste-heat  steam  generator)  for  further  investigation.  The  effect  of  flow 
distribution  control  on  the  baseline  waste-heat  boiler  performance,  under  both 
design  and  off-design  gas  turbine  operating  conditions  were  analyzed.  It  was 
estimated  that,  at  design  condition  without  flow  distribution  control,  the 
overall  heat  transfer  rate  would  be  approximately  16  percent  less  than  that 
obtainable  based  on  uniform  flow  distribution.  With  appropriate  flow  distribu¬ 
tion  control  (using  one  flow  guide  vane  and  one  flow  injection  for  boundary 
layer  separation  control),  the  boiler  efficiency  can  be  expected  to  improve 
by  approximately  20  percent  as  compared  with  that  of  the  uncontrolled  case. 
Based  on  che  results  of  this  analytical  study,  a  suggested  experiment  program 
was  formulated  for  ONR  consideration. 

This  study  program  was  conducted  by  the  Thermal  Engineering  Group  at  UTRC 
under  Contract  N00014-80-C-0478,  Modification  P00002,  from  the  Office  of  Naval 
Research,  Mechanics  Division,  Arlington,  Virginia. 
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RESULTS  AND  CONCLUSIONS 


The  design  of  a  gas  turbine  waste-heat  boiler  or  hot-water  heater  depends  on 
the  gas  turbine  mode!  to  which  it  would  be  mated,  its  end-use,  and  space  and 
economic  criteria.  Units  designed  for  industrial  applications  have  been 
custom-built  to  fit  different  configurations  using  mostly  finned  carbon  steel 
tubes ,  ' 


For  naval  propulsion  applications,  a  once-through  forced-circulation  steam 
generator  design  should  be  selected  because  of  stability,  reliability,  compact¬ 
ness  and  lightweight  considerations.  In  order  to  achieve  maximum  performance, 
the  gas-side  pressure  loss  for  the  steam  generator  should  be  limited  to  200  mm 
water-gage,  and  the  pinch-point  temperature  should  not  be  less  than  50°F, 

,  The  analytical  model  developed  to  predict  the  waste-heat  boiler  performance 

is  based  on  the  use  of  compact  heat  exchanger  design  criteria  and  the  relaxation- 
approach  method.  The  model  is  capable  of  estimating  the  waste-heat  boiler 
performance  at  any  inlet  gas  flow  distribution. 

.  Results  of  an  extensive  parametric  performance  analysis  indicate  that  among 
the  four  candidate  tube  size  anti  arrangements  combinations,  a  circular  finned 
tube  with  the  following  dimensions  is  the  most  effective  for  the  baseline 
waste-heat  boiler  design:  tube  length  =200  ft;  outside  diameter  =0.774  inch; 
fin  diameter  =1.403  inches;  fins  per  inch  =9;  fin  arrangement:  staggered  with 
longitudinal  pitch  =1.75  inches  and  transverse  pitch  =1.557  inches. 

.  At  its  design  condition  (corresponding  to  a  50-percent  power  output  of  the  gas 
turbine),  the  baseline  waste-heat  steam  generator  with  a  uniform  gas  flow 
distribution  is  estimated  to  be  able  to  generate  approximately  28000  pound  per 
hour  of  superheated  steam  at  700°F  and  300  psia.  At  this  condition,  the 
calculated  overall  heat  transfer  rate  would  be  approximately  10000  Btu/sec;  the 
gas-side  pressure  loss  would  be  0.55  psia;  and  the  pinch-point  temperature 
would  be  approximately  75*F. 

.  When  the  water  flow  rate  of  the  baseline  waste-heat  steam  generator  is  maintained 
at  its  design  value  of  7.9  Lb/sec,  the  steam  temperatures  with  and  without  flow 
distribution  controls  are  estimated  to  be  725*F  and  450*F,  respectively.  When 
the  steam  temperature  is  maintained  at  its  design  value  of  70G*F,  the  water 
flow  rates  with  and  without  flow  distribution  control  would  be  8.1  Ib/sec  and 
6.6  Ib/sec,  respectively. 
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To  provide  adequate  technical  information  so  a  comparison  with  the  analytical 
results  can  be  made  and  to  produce  operational  experience  with  a  gas  turbine 
waste-heat  steam  generator  in  naval  propulsion  applications,  an  experimental 
program  should  be  undertaken.  A  suggested  program  consisting  of  nine  major 
tasks  would  require  approximately  sixteen  months  to  complete  and  a  level  of 
effort  of  approximately  3000  man-hours. 
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INTRODUCTION 


As  a  result  of  a  system  feasibility  study  conducted  by  NAVSEA  in  1977 
(Ref,  1)  the  Rankine  Cycle  Energy  Recovery  (RACER)  system  was  selected  as  a 
candidate  system  for  future  advanced  naval  propulsion  system.  The  RACER  system 
uses  waste  heat  recovery  from  the  exhaust  of  the  marine  gas  turbines  to 
provide  additional  propulsion  power  for  the  U.S.  Naval  combatants.  Since  then, 
the  development  of  a  reliable,  efficient  and  compact  waste  heat  steam  generator 
has  become  one  of  the  most  important  engineering  disciplines.  Critical  technology 
areas  were  defined  and  appropriate  programs  were  initiated  to  address  these  and 
thereby  to  reduce  the  risk  of  system  development  (Refs.  2  and  3).  Results  of 
these  critical-technology  programs  indicate  that  a  self-cleaning  boiler  is 
feasible,  a  low-leakage  system  can  be  demonstrated,  and  IN625  or  IN825  would  be 
the  candidate  material  for  construction  of  the  waste-heat  boiler.  Based  on 
the  general  design  objectives  outlined  by  the  Navy  and  on  the  results  of  these 
critical  technology  programs,  contracts  for  the  preliminary  design  of  the  RACER 
system  were  awarded  in  1979  and  those  for  its  development,  testing  and  evaluation 
were  awarded  in  1981  (Ref.  4). 

Although  the  results  of  system  studies  and  critical  technology  programs 
continue  to  support  the  use  of  RACER  system  for  Naval  propulsion  applications, 
some  problems  related  to  the  general  design  practices  remained  Co  be  solved  by 
the  design  engineers.  Because  each  component  of  the  RASER  system  must  be 
designed  to  satisfy  specific  system  performance  requirements,  and  particularly 
those  related  to  the  waste  heat  steam  generators,  uncertainties  related  to 
heat-transfer  and  pressure- loss  coefficients,  as  well  as  to  nonuniform  flow 
distributions  must  be  eliminated.  Although  the  degree  of  nonuni formi ty 
and  its  effect  on  the  waste-heat  boiler  performance  are  not  completely  clear, 
what  is  apparent  is  that  the  waste  heat  boiler  must  be  designed  with  care 
because  of  such  factors  as  cost,  the  space  and  weight  limitations,  and  per¬ 
formance  and  reliability  requirements.  Therefore,  an  experimental  progt am 
becomes  a  necessity.  However,  for  large  units  like  the  RACER  system,  it  would 
be  unpractical,  if  not  impossible,  to  build  a  full-scale  test  apparatus  to 
conduct  a  comprehensive  test  program.  Accordingly,  the  present  analytical 
study  was  conducted  first  to  provide  some  basic  understanding  of  the  flow 
distribution  characteristics  and  the  effect  this  flow  has  on  marine  gas  turbine 
waste-heat  boiler  performance.  Based  on  the  analytical  results,  a  desirable 
and  constructive  experiment  program  can  be  formulated  for  ONR  consideration. 

It  is  well  understood  that  any  nonuniform  flow  distribution  will  reduce 
the  heat  transfer  performance  and  at  the  same  time,  increase  the  pressure  loss 
in  a  heat  transfer  device  to  various  degrees,  depending  on  specific  design  and 
actual  operating  condition.  Several  studies  have  been  made  in  the  past  to 
investigate  the  effect  of  flow  distribution  nonuniformity  on  the  heat  exchanger 
performance  (Refs.  5  to  8).  Because  the  actual  flow  distribution  would  be 
different  from  one  design  to  another,  these  studies  were  made  based  on 
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arbitrarily  assumed  nouniform  flow  profiles  for  the  working  fluids.  Results  of 
these  studies  indicated  that  as  much  as  30  percent  reduction  in  overall  heat 
transfer  unit  (NTU)  could  be  ascribed  to  the  poor  flow  distribution  in  the  heat 
exchanger  core.  Because  the  flow  distribution  profiles  assumed  in  these 
studies  are  quite  different  from  that  of  the  marine  gas  turbine  exhaust,  and 
furthermore  the  heat  exchanger  core  considered  are  often  unsuitable  for  naval 
piopulsion  system  applications,  these  study  results  can  be  used  for  reference 
purposes  only,  but  not  suitable  for  direct  applications.  Therefore,  an  analytical 
study  of  waste-heat  boiler  performance  based  on  actual  flow  distributions 
measured  in  a  typical  marine  gas  turbine  exhaust  was  performed  and  the  results 
obtained  are  presented  in  this  report. 

The  overall  analytical  program  has  been  structured  into  two  phases.  Phase 
I  (Ref.  9)  emphasizes  the  understanding  of  the  basic  flow-distribution  phenomena 
and  its  impact  on  two-dimensional  diffuser  design  and  performance.  Results  of 
the  Phase-I  study  indicate  that  flow  distribution  in  marine  gas  turbine  exhaust 
was  highly  irregular  and  nonuniform,  and  that  this  flow  will  remain  nonuniform 
through  a  two-dimensional  diffuser  unless  proper  flow  distribution  control 
meant,  are  useu.  This  nonuniform  flow  distribution  can  be  made  more  uniform  by 
using  a  specially  designed  diffuser  which  incorporates  appropriate  guide  vanes 
and,  if  necessary,  flow  injection  at  critical  locations.  The  results  of 
Phase-I  study  were  then  used  in  this  Phase-II  study  which  emphasizes  the  effect 
of  nonuniform  flow  distribution  on  the  waste-heat  boiler  performance. 

This  report  presents  the  technical  approach  and  the  results  of  an  analytical 
study  of  flow  distribution  control  in  marine  gas  turbine  waste-heat  steam 
generators.  The  report  consists  of  three  sections  and  one  appendix.  In 
Section  I,  the  applicable  steam  generator  design  concepts  and  general  design 
considerations  are  reviewed;  the  design  data  used  by  many  manufacturers  of 
waste-heat  boiler  are  evaluated;  the  critical-problem  areas  associated  with 
design  of  marine  waste  heat  steam  generators  are  discussed;  and  a  candidate 
waste  heat  steam  generator  configuration  was  selected.  Section  II  discusses 
the  analytical  model  formulation  and  presents  the  results  of  parametric 
performance  analysis,  including  those  for  both  design  and  off-design  operations 
of  candidate  waste  heat  steam  generators.  Based  on  the  results  of  this  analytical 
study,  a  proposed  experiment  program  plan  and  schedule  was  prepared;  this  is 
presented  in  Section  III.  The  detailed  descriptions  of  the  computer  program 
developed  for  the  analytical  model  are  presented  in  Appendix  A. 
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SECTION  1 

SELECTION  OF  CANDIDATE  STEAM  GENERATOR  CONFIGURATIONS 


Gas  turbine  waste  heat  recovery  systems  have  been  designed  and  used 
with  success  for  generating  either  hot  water  or  steam  or  both  for  various 
applications  (Refs.  1.1  to  1.4).  In  addition,  results  of  technical  and  economic 
feasibility  studies  have  shown  that  the  combined  gas  and  steam  turbine  power 
system  is  attractive  for  use  in  marine  propulsion  applications  (Refs.  1.5  and 
1.6).  The  use  of  small-scale  heat  exchanger  units  for  recovering  waste  heat 
from  service  gas  turbine  generators  of  U.S.  DD963  ships  has  also  been  reported 
in  Refs.  1.7  and  1.8.  Furthermore,  the  research  and  development  efforts  leading  to 
a.i  efficient,  lightweight,  and  reliable  waste  heat  recovery  system  for  U.S. 

Navy  surface  combatant  propulsion  application  is  underway  (Refs.  1.9  and 
1.0).  Accordingly,  the  objective  of  this  task  was  to  select  a  candidate  waste  heat 
steam  generator  configuration  which  could  be  integrated  with  the  candidate  gas 
turbine  and  diffusers  investigated  in  the  Phase-1  study.  In  order  to  achieve 
this  objective,  the  applicable  steam  generator  design  concepts  and  general 
design  considerations  were  reviewed;  the  design  data,  which  include  the  system 
operating  conditions  (flow  rate,  temperature,  and  pressure),  the  boiler,  its 
efficiency,  and  the  tube  macetial,  used  by  manufacturer  of  waste-heat  boiler  or 
hot-water  heacers  were  evaluated;  and  the  critical  problem  areas  associated 
with  design  of  marine  waste  heat  boilers  were  investigated.  Based  on  this 
information,  a  candidate  waste  heat  steam  generator  configuration  was  selected. 


1.1  Design  Concepts  and  Considerations  of  Marine 
Gas  Turbine  Waste-Heat  Boilers 

In  the  process  of  specifying  a  marine  gas  turbine  waste  heat  boiler,  a 
designer  must  determine:  (1)  the  total  amount  cf  heat  that  can  be  recovered 
economically;  and  (2)  the  type  of  equipment  that  is  best  suited  to  the  available 
space  and  the  quality  of  the  steam.  Based  on  the  results  obtained,  the  designer 
proceeds  to  investigate  other  basic  design  considerations  the  most  important  of 
which  are  summarized  in  Table  1.1.  These  considerations  are  based  on  general 
design  practices  of  industrial  gas  turbine  waste  heat  boiler  designs,  or  on  the 
general  constraints  and  requirements  of  naval  ship  operations. 

1.1,1  Design-Point  Performance  Considerations 

Capacity  sizing:  in  practice,  the  design  of  marine  gas  turbine  waste  heat 
boiler  is  conducted  on  one  of  two  approaches.  The  first  is  to  design  the 
system  for  a  ship  which  would  operate  at  full  load  for  long  periods  of  time 
(such  as  commercial  marine  and  naval  auxiliary  ships);  the  other  is  to  design 
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the  system  for  efficient  operation  at  cruise,  but  still  taking  into  consideration 
the  need  to  operate  for  intermittent  periods  at  full-power  (such  as  the  naval 
combat  ships).  The  method  of  integrating  the  waste  heat  steam  generator  with 
the  marine  gas  turbine  will  depend  on  both  the  ship  type  (maximum  installed 
power  and  duty  cycle)  and  the  gas  turbine  engine  selected.  If  two  gas  turbines 
were  needed  to  power  one  propeller,  it  would  be  desirable  to  have  the  exhaust 
systems  of  these  two  turbines  directed  through  a  single  waste-heat  boiler  for 
the  steam  turbine,  thereby  reducing  weight  of  the  waste-heat  recovery  system. 

In  tnis  case,  the  maximum  heat  which  could  be  recovered  will  depend  on  the 
performance  characteristics  of  the  gas  turbines  and  their  operating  profiles. 

The  off-design  performance  of  the  steam  cycle  will  depend  on  whether  it  is 
designed  for  cruise-  or  full-power  operation. 

Flow  Parameters:  Because  the  ratio  of  gas  to  liquid  (water)  flow  rates  in 
the  waste-heat  boiler  or  hot  water  heater  are  inherently  high,  externally  extended 
(finned)  tubes  are  more  desirable  than  bare  tubes.  Many  studies  (Ref.  1.2,  1.7, 
1.11  and  1.12)  indicate  that  when  the  gas  flows  across  the  finned  side  of  the 
tubes,  the  heat  transfer  will  be  maximized,  and  therefore,  the  designers  of 
most  gas  turbine  recovery  systems  have  adopted  this  cross  flow  pattern.  To  handle 
the  relatively  large  amount  of  gas  flow  at  low  pressure  and  high  temperature  and 
to  satisfy  the  low-pressure-drop  requirement,  the  flow  area  on  the  gas  side  must 
be  adequate.  When  external  finned  tubes  are  considered,  gas-side  pressure  drop 
through  the  heat  recovery  system  may  impose  significant  penalties  on  the  operation 
of  the  gas  turbine.  In  industrial  waste  heat  boilers,  the  pressure  drop  is 
normally  limited  to  approximately  15  inches  of  water  (0.6  psia).  Therefore,  the 
tube  size  and  tube  arrangement  must  be  carefully  selected.  To  increase  the  flow 
area  and  heat  transfer  area,  the  use  of  a  suitable  diffuser  to  connect  the  waste 
heat  boiler  with  the  gas  turbine  exhaust  box  becomes  necessary.  The  candidate 
diffuser  identified  during  Phase-I  study  will  serve  this  purpose. 

Pinch  Point:  It  is  difficult  to  assess  practical  limits  on  the  degree  of  heat 
recovery  without  considering  the  cost  of  the  equipment,  and  one  of  the  most 
important  parameter  in  sizing  the  waste  heat  boiler  is  the  pinch  point  tempera¬ 
ture.  Figure  1.1  shows  the  profile  of  the  turbine  exhaust  gas  and  the  water/ 
steam  temperature  for  a  typical  unfired  waste-heat  steam  generator.  The  pinch 
point  generally  occurs  where  the  liquid  reaches  its  saturated  state.  The 
selection  of  the  pinch-point  temperature  not  only  effects  the  liquid-side  flow 
condition  (flow  rate  and  pressure),  but  also  the  boiler  size.  As  indicated  in 
Ref,  1.11,  waste-heat  boilers  with  pinch-point  less  than  50*F  are  normally  not 
considered  to  be  economical. 

Temperature  Differential:  Unlike  that  in  a  conventional  oil  or  coal 
fired  boiler,  the  temperature  differential  between  the  two  working  fluids  in  a 
gas  turbine  waste-heat  boiler  is  low.  Accommodating  this  low  temperature 
differential  requires  a  special  design  in  terms  of  tube  arrangement  and  material 
selection.  Finned  tubes  with  high  thermal  conductivity  can  be  considered  as 
long  as  the  sum  of  the  material  cost  and  manufacturing  cost  do  not  exceed  the 
economic  limit. 
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1.1.2  Tube  Design  Considerations 


After  the  flow  conditions,  tube  material,  tube  size,  and  the  tube  arrangement 
have  been  selected  and  defined,  the  waste-heat  steam  generator  performance  can 
be  estimated;  the  selection  of  the  tube  size  and  tube  arrangement  has  a  significant 
effect  on  boiler  performance  and  size.  Use  of  small  diameter  tubes  yields  a  high 
heat  transfer  coefficient  on  both  sides  of  the  working  fluids  and  results  in  small 
boiler.  The  advantage  can  be  taken  of  using  small  tubes  only  when  working  with 
organic  fluids  or  extremely  high  quality  water  so  hardness  or  fouling  problems  are 
eliminated.  However,  from  a  practical  standpoint,  the  tube  size  selected  should 
be  sufficiently  large  to  accommodate  a  pneumatic  tube  reamer.  This  precaution  is 
taken  so  that  if  untreated  water  is  used  for  an  emergency  condition,  or  if  the 
cooling  water  (river  water  or  sea  water)  leaks  into  the  condensate,  the  tubes  can 
be  cleaned  mechanically  if  chemical  cleaning  is  impossible  or  if  the  tubes  become 
plugged  to  the  point  where  chemicals  cannot  be  introduced.  Therefore,  from  a  prac¬ 
tical  viewpoint,  tubes  smaller  than  3/4  inch  in  diameter  should  not  be  considered. 

Heat  exchanger  tubes  should  be  arranged  in  such  a  manner  that  thermal 
stress  concentration  can  be  avoided;  both  U-shape  and  coil  arrangements  are 
good  candidates  in  this  respect.  However  these  arrangements  are  not  generally 
regarded  as  being  compact  and  their  accessibility  for  maintenance  and  replacement 
of  parts  is  generally  poor.  A  modular  design,  similar  to  the  evaporator  of  the 
automobile  air  conditioning  unit  (with  finned  straight  tubes  used  as  the  heat 
transfer  core  and  with  the  ends  of  the  tubes  welded  to  a  U-shape  tube  joints 
which  are  located  outside  of  the  tube  sheet  as  shown  in  Fig.  II. 7),  may  be  a 
better  choice  for  marine  gas  turbine  waste  heat  applications. 

The  baffles  needed  to  act  as  tube  support  plates  and  flow  guide  vanes  must 
be  located  so  that  the  maximum  tube  length  between  support  plates,  or  between  a 
tube  sheet  and  a  supporting  plate,  does  not  exceed  36  inches.  Holes  for 
tubes  in  baffles,  baffles  clearances,  and  tie  rod  standards  must  be  designed  in 
accordance  with  the  latest  standards  of  the  ASME  boiler  design  code  (Ref.  1.16). 

The  selection  of  tube  material  affects  not  only  on  the  heat  transfer 
performance  and  the  initial  cost,  but  also  the  boiler  reliability  and  its 
operation.  For  landbased  waste-heat  recovery  systems,  carbon  steel  or  low 
alloy  are  commonly  used.  However,  for  naval  ship  propulsion  system  applications, 
high-temperature  stainless  steel  (such  as  304  or  316)  or  Incoloy  800  may  be 
used  to  cope  with  the  possible  dry-running  conditions. 

1.1.3  Performance  Degradation  C  ,s iderat ions 

In  designing  heat  transfer  equipment,  the  possible  performance  deterioration 
due  to  flow  leakage,  nonuniform  flow  distribution,  and  fouling  must  be  considered. 
Leakage  is  one  of  the  most  exasperating  problems  in  heat  exchanger  fabrication  and 
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maintenance,  for  it  not  only  effects  the  heat  exchanger  performance,  but  also 
requires  flow  make-up  and  clean-up  equipment.  Therefore,  use  of  all-welded 
tubes  may  be  considered  to  ease  the  leakage  problem,  Honuniform  flow  distribu¬ 
tion  has  some  effects  on  heat  exchanger  performance  (Refs.  1.14  and  1.15).  In 
Ref,  1.15  i'_  is  indicated  that  poor  flow  distribution  through  the  cores  of  a 
typical  counterflow  exchanger  can  cause  degradation  in  excess  of  30  percent  in 
the  operating  effectiveness  as  compared  with  the  values  predicted  for  the  ideal 
use  of  uniform  flow  distribution.  Therefore,  applicable  flow  distribution 
ivnlrol,  wherever  is  necessary,  must  be  incorporated  into  the  design  of  a  waste 
heat  boiler  to  avoid  any  unnecessary  performance  degradations. 

Fouling  has  also  been  a  problem  common  to  all  waste  heat  recovery  equipment. 
In  order  to  design  a  marine  waste  heat  boiler  capable  of  sustaining  its  design 
lapability  over  a  desired  period  of  operation  with  minimum  maintenance  and 
repair,  the  designer  must  give  serious  consideration  to  the  selection  of 
materials.  The  material  specified  must  be  able  to  offer  maximum  resistance  to 
corrosion,  and  to  the  fouling  characteristics  of  the  fluids  being  handled. 

As  an  added  design  burden,  consideration  must  be  given  to  the  varying 
degrees  of  inclination  encountered  in  sea  service.  In  naval  practice,  all  heat 
transfer  equipment  must  be  designed  to  perform  satisfactorily  under  conditions 
of  5  degrees  trim,  10  degrees  pitch,  15  degrees  list,  and  45  degrees  roll. 

1.1.4  System  Layout  Considerations 


The  physical  arrangement  of  the  gas  turbine  exhaust  relative  to  the 
location  of  the  heat  recovery  unit  has  considerable  effect  on  turbine  main¬ 
tenance  as  well  as  the  cost  of  the  overall  installed  recovery  system.  For 
industrial  applications,  the  horizontal  side-discharge  gas  turbine  exhaust 
(Fig.  1.2)  is  preferred.  This  arrangement  provides  good  access  for  turbine 
maintenance,  has  less  structural  support  for  the  waste  heat  recovery  components, 
and  provides  adequate  space  for  bypass  stack  and/or  supplementary  firing.  This 
arrangement  also  offers  the  opportunity  to  use  natural  circulation  (through 
vertical  tube  arrangement)  for  reliable  flow  circulation  and  uniform  heat 
distribution,  both  of  which  are  of  particular  importance  if  supplementary 
firing  were  required.  Because  of  space  and  weight  constraints  in  ship  propul¬ 
sion  system  applications,  a  veitical  top-discharge  gas  turbine  exhaust  (Fig, 

1.3)  is  more  desirable  (Refs.  1.9,  1,10,  and  1.12).  The  advantages  of  this 
arrangement  are  primarily  for  saving  in  cost  and  space  as  well  as  good  exhaust  gas 
distribution  across  the  boiler  heating  surface  provided  that  diffuser  is  properly 
designed.  Generally  speaking,  this  arrangement  does  not  create  any  special  gas 
turbine  maintenance  problems  because  the  engine  is  housed  in  its  own  enclosure  and 
can  be  removed  through  the  intake  for  major  services.  For  easy  installation, 
maintenance  and  replacement,  the  heat  exchanger  tube  elements  are  arranged  hori¬ 
zontally.  However,  this  arrangement  would  require  forced  circulation  of  liquid  to 
satisfy  such  operational  requirements  as  ease  of  control  and  dry-running.  In 
addition  to  having  the  characteristic  of  good  stability  and  reliability,  the 
forced  circulation  design  is  known  to  be  more  compact  and  lighter  in  weight  in 
comparison  with  the  natural  circulation  (vertical)  design. 
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The  physical  location  for  the  auxiliary  components,  such  as  the  pump,  the 
automatic  (pneumatic)  control  devices,  the  feed  vater  treatment  system,  and  the 
pipings  must  be  carefully  selected  so  that  accessibility  for  their  maintenance 
and  parts  replacement  is  adequate.  If  the  space  limitation  were  such  that  an 
integral  waste-heat  boiler  system  would  cause  problems  in  accessibility, 
dispersed  arrangement  of  some  secondary  auxiliary  component  must  be  considered. 

1.1.5  Structural  Rigidity  Considerations 

Due  to  the  rough  sea  service  condition,  all  components  of  the  waste- 
heat  boiler  must  be  provided  with  adequate  foundation  supports.  Additional 
allowances  must  be  made  in  the  design  of  the  heat  recovery  equipment  supports 
to  provide  for  expansion,  contraction  and  high-impact  shock.  Furthermore,  all 
design  features  must  conform  to  the  ASME  Boiler  and  Pressure  Vessel  Code  (Ref. 
1.16).  The  design  data  commonly  used  by  manufacturers  of  industrial  waste-heat 
boiler  and  waste-heat  economizer  (hot  water  heater)  and  the  critical  problem 
associated  with  marine  gas  turbine  waste  heat  boiler  design  are  discussed 
in  the  section  which  follows. 


1.2  Waste  Heat  Boiler  Design  Data  and  Critical  Problems 

A  comprehensive  survey  was  performed  to  identify  the  state  of  the  art  of 
waste  beat  boiler  design  (including  the  gas  and  liquid  flow  conditions,  the  unit 
capacity,  the  efficiency,  and  the  materials  used)  and  the  critical  technology 
problem  areas.  The  data  obtained  from  this  survey  are  shown  in  Table  1,2,  and 
the  critical  technology  problem  areas  are  summarized  in  Table  1.3. 

It  was  discovered  that  there  are  more  than  one  hundred  waste-heat-boiler/ 
economizer  manufacturers  worldwide  and  those  shown  in  Table  1.2  represent  only 
a  few  of  this  total.  The  design  approaches  used  by  these  firms  have  been 
varied  depending  on  the  heat  sources,  amount  of  heat  which  can  be  recovered 
economically,  the  end-use  of  the  recovered  heat,  and  space  and  economic  concerns. 
The  left  column  of  Table  1.2  shows  that  most  manufacturers  can  provide  custom- 
built  units  (shown  with  an  affixed  mark)  to  meet  a  specific  design  require¬ 
ment.  Therefore,  the  design  data  obtained  vary  over  a  wide  range.  For  example, 
on  the  fourth  line  SA  Babcock  Belgium  NV  can  provide  both  waste  heat  boiler  and 
waste  heat  economizer  (hot  water  heater)  for  gas  flow  rates  ranging  from  27  to 
333  cu.  meter/sec,  gas  temperatures  from  400  to  700*C,  and  gas-side  pressure 
losses  from  20  to  60  mm  W.G.  The  liquid  flow  rate,  liquid  temperature,  and 
unit  capacity  would  then  vary  according  to  the  design  requirements.  The  unit 
capacities  and  heat  recoverying  efficiencies  of  these  designs  vary  from  10  to 
200  MW  and  60  to  75  percent,  respectively.  The  design  data  for  other  manufac¬ 
turers  are  similar  in  nature,  but  different  in  level  of  absolute  values. 
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Among  the  design  data  presented,  the  gas-side  pressure  loss  information  is 
probably  the  most  useful  to  this  present  study.  It  was  found  that  (from  column 
No.  4  of  Table  1.2)  a  pressure  loss  between  100  to  200  mm  W.G.  would  be  a 
practical  value  for  marine  gas  turbine  waste  heat  boiler  design.  Other  informa¬ 
tion,  such  as  tube  materials  (shown  on  the  far  right  column  of  the  same  table) 
and  the  boiler  design  configuration  (not  shown  in  the  table),  offer  further 
insight  into  boiler  design.  For  industrial  applications,  finned  carbon  steel 
tubes  aie  the  most  commonly  used  although  stainless  steel  tubes  are  also  used 
in  some  designs.  The  boiler  configurations  are  mostly  once-through  designs. 

The  critical  technology  problem  areas  of  waste-heat  boilers  are  listed 
in  Table  1.3.  These  critical  problem  areas  are  generally  related  to  material 
selection,  mechanical  design,  or  operational  requirements.  From  available 
information,  it  appears  that  problems  with  materials  are  the  most  common  and 
serious  of  these  observed  in  the  steam  generator  equipment.  The  commonly 
encountered  material  problems  are  related  to  corrosion  damage  in  the  boiler 
tubes  including  those  of  denting,  pitting,  cracking,  and  erosion.  These 
problems  retults  from  the  attack  of  concentrated  aggressive  chemical  impurities 
on  the  tube  materials.  Laboratory  tests  made  on  samples  removed  from  dented 
steam  generator  tubes  indicate  that  denting  is  an  acid  chloride  reaction  (Ref. 
1.17).  Tube  wall  thinning  has  also  been  observed  in  the  region  near  the  tube 
sheet  when  phosphates  have  been  used  in  water  treatment.  Examination  of  tubes 
removed  from  a  once-through  steam  generator  has  also  revealed  that  stress 
corrosion  cracking  may  result  from  sulfuric  acid  attack.  In  view  of  these 
severe  material  problems.  Navy  initiated  a  material  study  program  in  FY79  at 
DTNSRDC/Anapol is  to  determine  the  best  material  for  use  in  boiler  tubes.  The 
results  of  the  study  indicate  that  Incloy  800  can  resist  to  oxygen  and  chloride 
stress  corrosion  as  well  as  that  from  sulfurous  and  sulfonic  acid  attack. 

The  cause  of  the  problem  of  tube  fretting  and  wear  can  be  traced  to 
vibration.  The  tube  vibration  can  be  induced  not  only  by  fluid  flow  perpen¬ 
dicular  to  the  tube  but  also  that  parallel  to  the  tubes.  Because  the  movement 
between  the  rubbing  surfaces  is  oscillatory  and  usually  small  in  amplitude,  the 
rubbing  process  taking  place  is  termed  "fretting".  It  is  well  known  that  the 
fretted  region  is  highly  susceptible  to  fatigue  cracks.  The  immediate  con¬ 
sequence  of  the  fatigue  cracks  is  the  leakage  of  working  fluid  and/or  cooling 
water.  The  leakage  in  the  boiler  affects  system  performance,  increases  feed- 
water  makeup  requirements,  and  demands  more  frequent  cleaning  (including 
deoxygenation)  operations.  Minimizing  the  flow-induced  vibration  is  a  critical 
task.  One  common  approach  has  been  to  use  an  all-tubular  boiler  with  no 
connections  other  than  the  water  inlet  and  steam  outlet  manifolds  to  reduce  the 
leakage.  However,  use  of  flow  distribution  control  to  achieve  a  more  uniform 
flow  distribution  is  even  more  essential,  and  in  the  long  run  may  prove  to  be 
the  most  beneficial  solution. 
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The  results  of  Phase-I  study  indicate  that  the  flow  distribution  within 
the  exit  diffuser  of  the  gas  turbine  exhaust  is  highly  nonuniform  in  the 
absence  of  flow  distribution  controls.  This  nonunifom  flow  distribution  not 
only  could  reduce  the  heat  transfer  performance  (Refs.  1.14  and  1.15),  but  also 
could  create  thermal  and  mechanical  stress  concentrations,  local  hot  spots,  and 
dryout  problems  (Ref.  1.18).  Therefore,  before  the  accurate  performance  can  be 
predicted,  both  analytical  and  experimental  programs  must  be  conducted  to 
investigate  the  actual  flow  distribution  pattern  inside  the  waste  heat  boiler 
predicted. 

Problems  related  to  transient  operation  of  the  marine  waste-heat  boiler 
must  also  be  addressed.  Because  of  the  self-cleaning  requirement,  the  waste-heat 
boiler  may  have  to  be  operated  under  dry  condition  for  a  period  of  15  to  30 
minutes  at  elevated  temperatures,  as  recommended  by  the  manufacturer.  Additional¬ 
ly,  the  boiler  has  to  be  operated  under  odd-design  condition  every  so  often  to 
meet  the  duty-cycle  requirements.  These  transient  and  off-design  operation  and 
the  routine  start-up  and  shut-down  procedures  will  undoubtedly  have  profound 
effect  on  the  boiler  reliability  and  life  expectancy.  However,  available  informa¬ 
tion  indicate  that  the  allowable  thermal  distortion  for  the  steam  turbine  would 
limit  the  rate  of  load  change  to  not  more  than  approximately  22  per  minute.  Even 
at  this  seemingly  slow  rate  of  load  change,  care  still  must  be  exercised  to 
control  the  boiler  and  its  auxiliary  system  so  that  the  pressure,  temperature,  and 
water  inventory  distributions  in  the  system  create  no  severe  conditions  throughout 
the  starting  period. 


X.3  Selection  of  Candidate  Steam  Generator  Configuration 

Based  on  the  results  of  Tasks  1.1  and  1,2  obtained  in  this  study,  and  the 
need  to  simplify  the  maintenance,  increase  the  reliability,  and  reduce  the  size 
and  weight  of  the  system,  a  once-through  cross-counterf low  type  boiler  was 
selected  as  the  steam  generator  configuration  for  this  analysis.  The  criteria 
used  in  this  selection  are  consistent  with  those  reported  in  Refs.  1.9  and 
I. 10,  and  therefore,  the  results  obtained  from  this  study  should  have  direct 
relevance  to  the  U.S.  Navy  RACER  Program.  A  sketch  of  the  conceptual  steam 
generator  configuration  is  shown  in  Fig.  1.4  and  a  summary  of  its  characteristics 
is  presented  in  Table  1.4. 

Since  the  objective  of  this  study  was  to  investigate  the  effect  of  gas  flow 
distribution  control  on  waste-heat  boiler  performance,  the  liquid  flow  has  been 
assumed  to  be  uniform.  As  shown  in  Fig.  1.4,  the  feedwater  would  be  supplied 
through  a  water  manifold  and  distributed  evenly  among  the  top  two  rows  of  the 
boiler  tubes;  the  superheated  steam  (or  hot  water)  would  be  discharged  from  the 
bottom  row  tubes  and  then  collected  in  a  steam  manifold.  The  gas  flow  would 
enter  the  boiler  from  its  bottom.  The  flow  distribution  and  flow  conditions  have 
been  based  on  the  results  of  the  Phase-I  study. 

In  the  section  which  follows,  an  analytical  model  of  flow  distribution 
control  formulated  for  the  candidate  steam  generator  is  presented,  and  results 
obtained  from  its  analysis  are  discussed. 
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TABLE  1.1 

CENERAL  DESIGN  CONSIDERATIONS  FOR  MARINE  GAS 
TURBINE  WASTE  HEAT  STEAM  GENERATORS 


.  Do  . ? gn  -Point  Performance  Considerations:  to  handle  the  large  amount  of 
gas  flow  at  low  pressure  and  high  temperature;  to  satisfy  the  gas-side  tow 
pressure  drop  requirement;  to  cope  with  low  temperature  differential  between 
the  two  working  fluids. 

,  Tube  Design  Considerations:  tube  size,  material,  and  arrangement, 

.  Performance  Degradation  Consideration:  flow  leakage,  nonuniform  flow  distri¬ 
bution,  and  fouling,  sea-service  condition. 

.  Economic  Considerations:  material,  sizing,  and  effectiveness. 

.  Operational  Requirement  Considerations:  thermal  load,  duty  cycle,  dry  running, 
emergency  operation,  and  control  devices. 

.  System  Layout  Considerations:  accessibility  for  maintenance  and  repair, 
reliability,  space  and  weight  limitations. 

.  Structural  Rigidity  Consideration:  shock  and  vibration,  structural  expansion 
and  contraction. 
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TABLE  1.3 

CRITICAL  PROBLEM  AREAS  OF  GAS  TURBINE  WASTE  HEAT  BOILERS 


.  Material  Problems:  Tube  denting,  pitting,  cracking,  erosion-corrosion 

.  Vibration  Problems:  Tube  fretting  and  wear,  high-cycle  fatigue,  stress 

corrosion 

Leakage  Problems:  Performance  degradation  and  feed-water  makeup  and 
cleanup  including  deoxygenation 

,  Flow  Maldistiibution  Problems:  Effective  flow  distribution  control, 

soot  formation  prevention  and  self-cleaning 
methods  , 

.  Transient  Behavior  Problems:  dry  cleaning  operation,  and  duty  cycle 

operation  of  gas  turbine,  regular  or  emergency 
shut-down  and  start-up. 
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TABLE  1.4 

CANDIDATE  STEAM  GENERATOR  CONFIGURATION 


.  Gas  Flow:  one-pass  flowing  upward  without  supplementary  firing 

.  Water  Flow:  once-through  forced  circulation,  counter  cross  to  gas  stream 

.  Tubes  and  Tube  Arrangement:  all-welded  finned  tubes  made  of  corrosion 

resistant  material  (IN  800);  placed  horizontally 
along  the  gas  turbine  centerline  direction 

.  Boiler  Geometry:  rectangular  corss-sect ion  of  10  ft  by  7  ft  (compatible  with 

the  diffuser  obtained  from  Phase-I  study),  height  be  less  than 
8  ft  (compatible  with  the  ship) 

.  Self-cleaning  on  the  gas  side  and  scaling  prevention  on  the  water  side 
.  Heat  Recovery  Capacity:  between  12,000  to  20,000  kw 
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SECTION  II 

PERFORMANCE  ANALYSIS  OF  FLOW  DISTRIBUTION 
CONTROT,  ON  WASTE-HEAT  STEAM  GENERATOR 

This  section  describes  the  formulation  of  an  analytical  model  and  the 
analysis  of  flow  distribution  control  on  marine  gas  turbine  waste-heat  boiler  (or 
steam  generator)  performance.  The  analytical  model  formulated  is  based  on  compact 
heat  exchanger  design  concept  (Ref.  2.1)  and  the  relaxation-approach  method 
(Ref.  2.2).  The  configuration  of  the  candidate  waste-heat  boiler  is  a  once- 
through,  counter  crossflow  heat  exchanger  selected  in  Section  1.3.  The  perfor¬ 
mance  characteristics  of  this  candidate  waste-heat  boiler  at  its  design  and 
otf-design  conditions  were  analyzed  for  three  different  flow  inlet  conditions 
obtained  from  Phase-I  study  (Ref.  2.3).  The  three  flow  inlet  conditions  are:  (1) 
a  uniform  flow  distribution;  (2)  a  nonuniform  flow  distribution  (based  on  actual 
flow  distribution  of  a  typical  marine  gas  turbine  exhaust)  without  flow  distribu¬ 
tion  control;  and  (3)  a  nonuniform  flow  distribution  with  flow  distribution 
cont  ro 1 . 

1 

II. 1  Formulation  of  Analytical  Model 

The  analytical  model  used  in  this  study  is  modified  version  of  the  distillate 
fuel  vaporization  model  originally  developed  by  Chiappetta  and  Szetela  (Ref.  2.2). 
Although  the  basic  assumptions  and  the  method  of  approach  for  the  fuel  vaporization 
model  are  adequately  described  in  Ref.  2.2,  additional  assumptions  as  well  as 
working  formulas  pertaining  to  this  present  study  are  presented  below  to  allow  a 
better  understanding  of  the  results  of  this  study. 

The  procedure  used  in  formulating  the  working  model  for  this  present  study 
consists  of  the  following  five  steps:  (1)  establishing  a  nodal  system  to  represent 
the  overall  waste  heat  boiler;  (2)  compiling  all  heat  transfer  and  pressure  loss 
working  formulas  and  empirical  constants;  (3)  developing  a  computer  program  for 
computing  the  thermal  and  physical  properties  of  water  substance  for  the  applicable 
flow  conditions;  (4)  establishing  a  numerical  computation  procedure  for  computa¬ 
tional  analysis  of  the  waste-heat  boiler;  and  (5)  developing  a  computer  program 
based  on  results  obtained  from  steps  (1)  through  (4)  to  facilitate  the  compu¬ 
tational  analysis  of  flow  distribution  control  in  marine  gas  turbine  waste-heat 
boilers.  A  detailed  discussion  of  each  step  follows. 
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1 1 .  1 . 1  Establishment  of  Heat  Exchanger  Nodal  System 

As  cited  in  Ref.  2.2,  in  order  to  use  a  nodal  system  for  easy  computation 
of  crossflow  heat  exchanger  performance,  several  assumptions  must  be  made. 

Two  basic  assumptions  made  are:  (i)  the  heat  exchanger  must  be  rectangular  and  of 
uniform  thickness;  and  (2)  the  working  fluid  on  the  shell  side  must  be  a  gas 
although  the  working  fluid  on  the  tube  side  can  be  either  a  gas  or  a  liquid. 

Based  on  these  two  assumptions,  the  overall  waste  heat  boiler  can  be  subdivided 
into  several  nodes  as  shown  in  Fig.  III.  1 .  Each  node  represents  a  rectangular 
parallelpiped  of  same  length,  i  ,  but  of  a  different  width,  Ax^  and  height,  Ay^. 

The  overall  waste-heat  boiler  then  can  be  described  by  a  two-dimensional  nodal 
array  in  i,  and  j,  where  i-1,2,3...  imax,  and  j=l,2,3...  jmax.  Nodes  are 
connected  to  form  several  groups  which  represent  the  flow  paths  and  flow  direction. 
For  example,  Fig.  II.  1  shows  a  5  by  10  array.  Assuming  that  the  gaseous  flow  for 
this  5  by  10  array  is  divided  into  five  paths  and  is  flowing  upwards,  then  the 
first  gas  path  can  be  represented  by  (1,1),  (1,2),  (1 , 3) , . . . (1 , 10) ,  and  the  second 
gas  path  can  be  represented  by  (2,1),  (2,2),  (2 ,3) , . . . (2 , 10) ,  etc.  If  the  water 
were  flowing  in  a  single  path  counter  cross  to  the  gas  stream,  then  the  nodal 
connection  for  the  water  path  can  be  represented  by  (5,10),  (4,10),  (3,10), 
...(1,10),  (1,9),  (2 ,9) . . . (5 , 1) .  The  flow  conditions,  including  the  mass  flow 
rate,  the  temperature,  and  the  pressure,  are  specified  at  the  inlet  of  each 
flow  path  to  provide  a  starting  point  for  numerical  computation. 

It  should  be  noted  that  each  node  is  treated  as  a  miniature  heat 
exchanger.  Because  of  the  nodal  arrangement,  the  exit  flow  condition  of  each 
node  will  automatically  become  the  inlet  flow  condition  of  the  subsequent  node. 

In  order  to  preserve  the  mass  flow  rate  for  each  gas  path,  a  plate-fin-type 
heat  exchanger  or  a  baffle  plate  tube  supporting  structure  must  be  specified  in 
the  construction  of  the  flow  path  for  the  waste  heat  boiler  if  experimental  and 
analytical  results  are  to  be  compared. 

The  heat  transfer  coefficient  for  the  "miniature"  heat  exchanger  is  computed 
based  on  the  averaged  local  flow  conditions.  Since  the  number  of  nodes  in  each 
path  and  the  number  of  paths  for  each  fluid  are  arbitrarily  selected,  proper 
arrangement  of  the  nodal  connection  can  be  made  to  simulate  any  two-dimensional 
heat  exchanger  geometry. 
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11.1.2  Heat  Transfer  and  Pressure  Loss  Working  Formula 

The  heat  transfer  coefficient  (convective  film  coefficient)  and  the  pressure 
loss  character ist ics  are  functions  of  the  surface  geometry,  fluid  properties,  and 
flow  conditions.  It  is  not  the  objective  of  the  present  study  to  develop  these 
functional  relationships,  but  rather  to  compile  the  existing  working  formulas 
from  the  open  literature  and  to  use  them  for  computing  the  heat  transfer  and 
pressuie  loss  characteristics  of  the  waste-heat  boiler  being  studied.  The  formula 
and/or  test  data  which  were  used  in  the  present  analytical  model  are  summarized  in 
the  following  sections. 

II. 1.2.1  Gas-Side_Working_Forraulas 

The  heat  transfer  and  pressure  loss  characteristics  for  the  gas-side  (shell- 
side)  working  fluid  are  based  on  the  data  presented  in  Ref.  2.1.  The  friction 
coefficient  (f)  and  the  Stanton  number  (St=h/C^  p(J)  were  expressed  in  terms  of 
Reynolds  number  for  various  surface  geometry  in  Chapter  IX  of  Ref.  2.1.  When 
extended  surfaces  were  considered,  the  fin-effect  formulas  presented  in  Chapter  II 
of  the  same  reference  were  also  used. 

II.  1.2.2  I^it^uid-Side^  Working  Formula 

Because  the  working  fluid  (water)  on  the  liquid  side  may  experience  phase 
changes  from  a  liquid  to  a  vapor  and  possibly  to  a  superheated  vapor,  it  was 
necessary  to  define  different  formulas  for  heat  transfer  and  pressure  loss  in  each 
phase.  In  addition,  the  working  fluid  might  be  in  laminar  flow,  turbulent  flow, 
or  supercritical  flow,  appropriate  working  formulas  had  to  be  used. 

For  liquid-phase  flow,  the  heat  transfer  coefficient  can  be  computed 
using  one  of  the  following  formulas: 


laminar  flow:  Nu 


turbulent  flow:  Nu 


supercritical  flow:  Nu 


C  /Rgb  ?rb\  C2  / \ C3 

l\  L/D  /  W/ 

C1(Reb)C2  (Prb)c3 


C1  (Rew>C2 


(la) 

(lb) 

(ic) 


where  C^,  C2,  C-j  and  C^  are  empirical  constants  presented  in  Table  II. 1. 

The  Reynolds  number  (Re),  Prandtl  number  (Pr),  viscosity  (ji)  and  flow  density  (p) 
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are  computed  based  on  either  the  bulk  temperature  or  the  wall  temperature  as 
designated  by  the  subscripts  of  b  or  w,  respectively.  The  pressure  drop,  Ap ,  can 
be  calculated  using  the  following  relationship  from  F.or.  2.1: 


(l-°2-Ke) 


P  . 
in 


P  „ 
ou  t 


(2) 


where  K£  and  Kg  are  entrance  and  exit  coefficients,  0  is  the  ratio  of  free¬ 
flow  area  to  frontal  area,  and  A^t  and  Ac  ^  are  the  heat  transfer  area  and  core 
flow  area,  respectively. 

For  boiling-phase  flow,  the  heat  transfer  coefficients  were  calculated 
using  the  correlations  developed  by  Chen  (Ref.  2.4),  who  assumed  that 
the  overall  boiling  heat  transfer  coefficient  consists  of  two  additive  basic 
mechanisms:  an  ordinary  macro-convective  mechanism  and  a  micro-convective 

mechanism  associated  with  bubble  nucleation  and  growth.  He  defined  these  two 
convective  heat  transfer  coefficients  as: 


hmac  ‘  °*023  (Rei)°-8  (Prt)°-8(^  F 


(At)0 .24  (4p)0 • 75  s  (3b) 


where  F  and  S  are  called  the  effective  two-p'nase  Reynolds  number  function  and 
the  bubble  growth  suppression  function,  respectively.  Both  F  and  S  are  determined 
from  empirical  correlations  of  heat  transfer  data  and  the  momentum-analogy  analysis 
and  are  presented  in  graphic  form  in  Ref.  2.4.  The  terms,  0  and  *  in  Eq.  (3b)  are 
the  vapor-liquid  surface  tension  and  latent  heat  of  vaporization,  respectively. 

The  difference  between  the  wall  temperature  and  the  saturation  temperatures,  Ax, 
and  Ap  is  the  difference  in  vapor  pressure  corresponding  to  this  At.  The  subscripts 
t  and  v  refer  to  liquid  and  vapor,  respectively. 

If  the  two-phase  boiling  process  were  isothermal,  the  pressure  loss  can  be 
estimated  using  the  model  developed  by  Lockhart  and  Martinelli  (Ref.  2.5).  In 
this  model,  the  overall  pressure  loss  consists  of  three  additive  components:  a 
gravitational  loss  (Ap^.),  a  momentum  loss  ( >  *nd  a  frictional  loss  (APp), 
which  are  defined  as: 


hmic  -  0.00122 


Kt°.79  Cpi0-*5  p^O.49  g0.25\ 


/Kt°-79  Cpt° 

\  <j0.5  WjO. 


29  *0.24  p  0.24 
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&PG  »■  (l-“)  Pj,  + 


inS  1 


sin 


(4a) 


a  P*U  A 

iPM  =  - T-  4 

M  1? 


%  = 


1  ♦  y2/n 


whore  a  and  Y  are  defined  as: 


a  = 


1  + 


1-X 


M  i 

X2  ' 

(4b) 

l\‘-7  °t 

+  “V 

'N  N  (&  p 

j  sv 

(4c) 

3/2| 

-1 

(5a) 

Y  = 


Revm  Ct  Pv 
Ret"  Cv  Pi 


(r) 


(5b) 


whore  m,  n,  Cp,  Cv,  N  are  constants  presented  in  Table  II. 1.  The  term,  x»  is  the 
steam  quality,  and  is  the  nodal  width  defined  in  Section  II. 1.1. 


For  the  supercritical  flow,  the  analytical  and  experimental  results  presented 
in  Chapters  V  and  VI  of  Ref.  2.1  were  used  to  calculate  the  heat  transfer  and  the 
pressure  losses  .  In  this  procedure,  the  friction  coefficient  and  Stanton  number 
were  expressed  in  terms  of  the  Reynolds  number  in  graphical  form,  which  were 
tabulated  as  input  data  to  the  computer  program  which  is  discussed  in 
Appendix  A. 


In  order  to  increase  the  flow  depth  to  satisfy  specific  design  re¬ 
quirements,  each  flow  path  may  be  reversed  alternately  to  form  several  passes. 
The  pressure  drop  associated  with  the  turns  were  calculated  using  the 
averaged  dynamic  head  evaluated  at  the  nodes  before  and  after  the  turn. 
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11. 1.3  Thermodynamic  and  Transport  Properties  of  Water  Substance 

Because  the  original  analytical  model  (Ref.  2.2)  was  developed  for  distillate 
fuel,  vaporization  application,  the  development  of  a  computer  program  which 
could  estimate  the  thermodynamic  and  transport  properties  of  water  at  any  flow 
condition  became  an  essential  part  of  the  analytical  model  for  the  present 
study.  In  this  model  the  thermodynamic  and  physical  properties  needed  are  the 
temperature,  pressure,  specific  volume,  enthalpy,  specific  heat,  viscosity,  and 
thermal  conductivity.  The  numerical  values  for  the  thermodynamic  properties 
(the  first  six  items)  can  be  obtained  from  a  fundamental  equation,  called  the 
Helmholtz  free  energy  equation,  which  is  described  in  Ref.  2.6.  The  advantage 
of  using  this  fundamental  equation  is  that  all  thermodynamic  properties  can  be 
obtained  from  its  derivatives.  Because  differentiation,  unlike  integration, 
produces  no  undetermined  functions  or  constants,  the  information  yield  is 
complete  and  unambiguous. 

To  calculate  the  thermal-conductivity  and  the  viscosity,  two  well-known 
relationships  contained  in  Refs.  2.7  and  2.8  were  used.  These  two  working  relation¬ 
ships,  along  with  the  derivatives  of  the  Helmholtz  free  energy  equations,  were 
then  incorporated  into  computer  programs  for  use  in  this  program. 

1 1 . 1 . 4  Numerical  Computation  Procedure 

The  numerical  computation  procedure  used  in  this  study  is  the  same  as  thqt 
presented  in  Ref.  2.2.  To  provide  a  better  understanding  of  this  analytical 
model,  a  brief  discussion  of  its  computational  procedure  is  given  as  follows. 

As  mentioned  in  Section  II. 1.1  that  each  subdivision  (called  node)  will  be 
treated  as  a  miniature  heat  exchanger.  The  performance  of  each  miniature  heat 
exchanger  will  be  calculated  based  on  the  averaged  temperature  for  each  fluid 
and  for  the  walls  in  each  node  during  the  previous  iteration  (to  be  described 
below).  Since  the  heat  transfer  areas  on  the  shell-side  and  tube-side  for  a 
finned-tube  bundle  are  not  equal,  the  overall  heat  transfer  coefficient  is 
conventionally  referenced  to  the  shell-side.  The  steady-state  heat  transfer 
rate  for  the  node  can  be  expressed  as: 


Qk  -  (U)K  (Ag)K  |(Tg)K  -  (Tt)KJ  <6) 

where  the  subscripts  g  and  £  denote  the  gas  and  liquid  sides,  respectively. 
The  term,  U,  represents  the  overall  heat  transfer  coefficient  determined  from 
the  film  coefficient  of  the  working  fluids  on  both  sides  of  the  tubes  and  the 
thermal  conductivity  of  the  tube  material. 
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The  averaged  fluid  temperatures  of  the  gas  and  liquid  (T^  and 
Tt )  are  functions  of  the  heat  transfer  rate  between  these  fluids.  For 
example  the  averaged  gas  temperature  is 


and 
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By  definition,  the  inlet  temperature  for  this  node  is  the  outlet  temperature  of 
the  previous  node  (or  Che  gas  supply  temperature  if  this  is  the  initial  node  for 
the  path  containing  this  node).  Similarly,  for  the  liquid  side: 


and 
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Substituting  Eqs.  (7)  through  (10)  into  Eq.  (6),  one  obtains 


(ig)m  cue*)  -  Tl)ra  d+Bj) 
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where  0  *  ^P^in/^P^out’  an^  t*,e  subscript  K  has  been  omitted  in  Eq.  (11) 
to  facilitate  typing. 


The  temperature  distributions  in  gas,  liquid  and  the  tube  walls  of  the 
waste  heat  boiler  were  calculated  using  a  relaxation  method.  This  method  is 
described  as  follows:  The  averaged  temperature  of  each  fluid  and  of  the  walls 
at  each  node  estimated  during  the  previous  iteration  was  used  to  calculate  the 
thermal  properties  and  the  heat  transfer  coefficients  for  each  fluid.  At 
the  completion  of  each  iteration,  the  outlet  temperatures  at  the  last  node  in  each 
path  on  each  side  of  the  working  fluid  were  compared  with  those  calculated  during 
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the  previous  iteration.  If  all  outlet  temperatures  were  within  a  specified 
tolerance,  the  iteration  was  then  terminated.  Otherwise,  the  new  averaged 
fluid  temperatures,  together  with  the  new  wall  temperatures,  were  used  as  bases 
for  the  next  iteration.  It  should  be  mentioned  that,  for  the  steady-state 
operation,  the  inlet  condition  were  used  as  the  initial  guesses  of  the  average 
temperature  distribution  for  each  path.  This  simple  iteration  algorithm  was 
shown  to  be  quite  stable  for  most  cases  analyzed  by  the  authors.  Typically 
convergence  occurs  within  fifteen  iterations  with  a  convergent  tolerance  of  5 
degree  F. 


II.  1.5  Development  of  Computer  Program 


In  order  to  minimize  the  computer  program  development  effort,  the  existing 
fuel  vaporization  heat  exchanger  program  was  updated  to  satisfy  the  study 
objective.  This  modification  was  made  with  permission  and  assistance  from 
Messrs.  Chiappetta  and  Szetela.  The  modification  included  allowance  for:  (1) 
isothermal  vaporization;  (2)  variable  nodal  size,  (3)  use  of  water  as  the 
working  fluid;  and  (4)  inclusion  of  overall  size  and  cost  estimates.  Descrip¬ 
tion  of  this  computer  program,  which  includes  the  program  capability  and 
input/out  format,  is  given  in  Appendix  A. 
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11. 2  Design  and  Off-Design  Flow  Condition 

The  shell-side  flow  distribution  and  flow  condition  at  the  inlet  to  the 
waste  heat  boiler  are  obtained  from  the  Phase-I  study  (Ref.  2.3).  In  Figs. 

II. 2  and  II. 3,  two  typical  diffuser  flow  distributions  are  shown  for  a  gas 
turbine  operating  at  50?i  power  with  and  without  flow  distribution  control.  The 
flow  distribution  and  flow  conditions  at  the  exit  of  the  diffuser  were  calculated 
for  eaih  of  the  several  sections  which  extend  the  entire  width  of  the  flow 
paths,  as  shown  in  Fig,  II. 1.  The  resulting  flow  conditions  are  tabulated  in 
Figs.  II. 4  and  II. 5  for  gas  turbine  operated  at  100-  and  50-percent  power, 
respectively.  Three  different  sets  of  flow  distribution  data  are  presented  for 
each  operating  condition,  namely,  a  uniform  flow  distribution,  an  actual  flow 
distribution  without  control,  and  an  actual  flow  distribution  with  control. 

For  the  first  two  data  sets,  the  flow  was  divided  into  five  paths  equally 
spaced  while  for  the  third  data  set,  the  flow  was  divided  into  six  paths 
because  a  pressure  difference  exists  between  the  two  sides  of  the  flow  guide 
vane.  The  averaged  flow  distribution  data  were  used  as  input  conditions  to  the 
parametric  performance  analyses  which  are  discussed  as  follows. 

11. 3  Parametric  Performance  Analyses 

Parametric  performance  analyses  of  the  waste  heat  steam  generator  were 
made  using  the  computer  program  discussed  in  Section  11.1.5.  The  parameters, 
which  were  varied,  included:  (1)  the  tube  arrangement,  (2)  the  effective  tube 
length,  and  (3)  the  water  flow  rate.  The  uniform  gas  flow  distribution  data 
presented  in  Figs,  II. 4  and  II. 5  were  used  as  a  reference  for  comparison  with 
the  results  of  the  nonuniforta  flow  cases.  The  range  of  the  parameters  used  were 
based  on  the  following  rationales. 

As  noted  in  Section  1.3,  the  candidate  steam  generator  selected  was  a 
once-through ,  counter  crossflow  design.  For  this  design,  optimal  heat  transfer 
surface  geometry  and  the  tube  arrangement  still  had  to  be  determined  in  order 
to  define  its  most  efficient  performance.  Four  different  tube  design  configura¬ 
tions,  depicted  in  Fig.  II. 6,  were  obtained  from  Ref.  2.1.  Each  configuration 
is  comprised  of  circular  finned  tubes  with  tube  diameter  equal  either  to  1.024 
inch  or  to  0.774  inch.  For  each  tube  diameter,  two  different  longitudinal  and 
transverse  pitches  and  fin  geometries  were  considered.  The  effect  of  the  tube 
designs  on  the  performance  of  the  candidate  waste-heat  boiler  were  then  inves¬ 
tigated. 

The  next  parameter  determined  was  the  effective  tube  length  which 
is  related  to  the  following  factors:  (1)  ease  of  installation  and  removal  in 
the  ship  (i.e.  maximum  height  of  each  heat  transfer  module  should  be  less  than 
8  ft);  (2)  acceptable  thermal  gradient  within  each  node  (assumed  to  be  equal  to 
or  less  than  25*F  for  numerical  stability);  and  (3)  design  and  manufacture  of 
boiler  tubes  according  to  the  ASME  Boiler  Codes,  With  those  concerns  in  mind, 
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the  maximum  nodal  width  would  be  approximately  24  inches,  and  the 

maximum  nodal  height,  (Ayj)  would  be  approximately  4  inches.  The  best  tube 
arrangement  which  would  meet  these  constraints  is  shown  in  Fig.  II. 7.  Based  on 
this  tube  arrangement,  the  maximum  effective  tube  length  was  estimated  to  be 
approximately  200  ft. 

The  maximum  water  flow  was  calculated  from  an  energy  balance.  From  the 
Phase-I  study,  it  was  determined  that  the  flow  rates  of  the  gas  turbine  engine 
exhaust  were  approximately  160  lb/sec  and  100  lb/sec  for  full-load  and  half¬ 
load  operation,  and  that  the  corresponding  temperatures  were  856*F  and  796°F, 
respectively.  If  the  gas  exit  temperature  from  the  waste  heat  boiler  were  kept 
above  300F  to  avoid  sulfuric-acid  condensation  (or  corrosion)  problems,  the 
maximum  amount  of  heat  that  could  be  recovered  would  be  approximately  21  MW  and 
12  MW  for  engines  operated  at  100-percent  and  50-percent  power,  respectively. 
Assuming  that  the  feed  water  enters  the  waste  heat  boiler  at  115. 7 0 F  (correspond¬ 
ing  to  a  condenser  pressure  of  3  inch  Hg. )  and  300  psia,  and  leaves  the  boiler 
as  a  superheated  steam,  the  maximum  water  flow  rates  for  a  gas  turbine  operated 
at  100-percent  and  50-percent  power  would  be  approximately  18  lb/sec  and  11 
lb/sec,  respectively.  In  parametric  performance  analyses,  the  water  flow  rates 
were  varied  in  increment  of  0.5  lb/sec  per  step  to  investigate  the  effect  of 
this  factor  on  the  performance  of  the  waste-heat  boiler. 

After  the  values  of  all  parameters  were  defined,  the  parametric  performance 
analyses  were  conducted.  The  results  of  these  analyses  (which  are  presented  in 
terms  of  steam  temperature,  gas  exit  temperature,  gas  side  pressure  loss,  and 
overall  heat  transfer  rate  as  function  of  water  flow  rate)  for  all  four  tube- 
design  configurations  and  for  gas  turbine  operated  at  full  power  are  shown  in 
Figs.  II. 8  through  II.  11.  The  water-side  pressure  loss  is  not  shown  because 
the  pumping  power  required  for  pressurizing  the  water  has  essentially  no  effect 
on  the  cycle  efficiency.  Results  shown  in  Figs.  II. 8  through  II. 11  indicate 
that  the  a-rangements  of  tube  design  indeed  have  a  significant  effect  on  the 
performance  of  the  waste-heat  boiler.  Generally  speaking,  Configuration  1 
would  perform  better  than  Configurations  3,  2,  and  4  in  that  order  if  boiler 
efficiency  were  the  only  concern. 

Figure  II. 8  shows  that,  for  a  given  tube  design  configuration,  the  water 
flow  rate  can  be  regulated  to  yield  a  wide  range  of  steam  temperatures  desired. 
However  if  the  water  flow  rate  were  greater  than  15.6  Ib/sec  for  Configuration 
No.  3  (Fig.  II. 6)  a  wet  steam  would  be  produced.  If  the  water  flow  rate  were 
too  low,  the  boiler  would  operate  below  its  attainable  efficiency  (see  Fig. 

II.  11).  The  experience  gained  at  UTRC  from  studies  of  waste-heat  recovery 
systems  indicates  that  a  700-degree  Fahrenheit  steam  at  approximately  300  psia 
woulc-  be  a  practical  design  for  a  Rankine-cyc le  power  conversion  system  applica¬ 
tion.  Therefore,  this  steam  condition  was  selected  as  reference  in  the  final 
selection  of  a  baseline  design  of  waste  heat  boiler  for  naval  applications. 
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In  order  to  examine  the  sulfuric-acid  condensation  problem  on  the  cold 
end  of  the  waste  heat  boiler,  the  gas  exit  temperatures  were  plotted  as  a 
function  of  the  water  flow  rate  for  all  four  tube  design  configurations;  these 
data  are  shown  in  Fig.  II-9.  The  constant  temperature  lines  (shown  in  dash¬ 
line)  were  obtained  from  Fig.  11.8.  It  is  seen  that  for  an  effective  tube 
length  of  200  ft,  the  gas  exit  temperature  would  be  between  370°F  and  480*?,  or 
*1  i*'  ve  the  sulfuric  corrosion  formation  temperature  of  300°F,  Another 
implication  of  these  data  is  that  it  is  possible  to  improve  the  waste-heat 
boiler  performance  by  increasing  the  overall  heat  transfer  area  if  there  is  no 
space  limitation  and  if  the  gas-side  pressure  loss  can  be  tolerated. 

The  effects  of  tube  design  conf igurat ion  and  water  flow  rate  on  the 
,’,i.=  -side  pressure  loss  and  on  the  overall  heat  transfer  rate  are  shown  in  Figs. 

II.  10  and  II.  11,  respectively.  It  can  be  seen  from  these  figures  that  he 
pressure  loss  varied  between  0.7  psia  and  2.0  psia  for  the  parameter  ranges 
considered.  These  pressure  loss  values  were  used  in  identifying  the  correction 
facto  for  gas  turbine  output  power  which  will  be  discussed  later  (Section 
11.5). 

II. 4  Off-Design  Performance  Analyst s 

In  order  to  compare  the  performance  characteristics  of  a  given  waste  heat 
boiler  design  at  different  operating  conditions,  the  uniform  flow  distribution 
data  presented  in  Fig.  II.. 5  for  gas  turbines  operated  at  50-percent  power  were 
also  considered.  The  results  of  this  parametric  performance  analysis  were 
compared  with  those  obtained  previously  based  on  the  design-point  conditions. 

In  Figs.  11.12  through  11.15,  these  comparisons  are  shown  in  terms  of  steam 
temperature,  overall  heat  transfer  rate,  gas  exit  temperature,  and  gas-side 
pressure  loss,  respectively.  The  solid  lines  are  results  for  gas  turbine 
operated  at  100-percent  power  and  the  dotted  lines  are  for  50-percent  power 
cases.  The  relationships  of  performance  to  tube  design  configuration  is 
presented  in  detail  in  the  following  Section  (II. 5). 

Figure  11.12  shows  that  if  these  candidate  waste-heat  boilers  (Configurations 
l  through  4)  were  integrated  with  a  given  gas  turbine  which  was  operated  at 
half-load,  the  water  flow  rate  must  be  reduced  significantly  in  order  to 
generate  the  same  quality  steam  (i.e.  with  the  same  steam  temperature).  For 
example,  'f  the  steam  temperature  required  is  700*F,  the  water  flow  rate 
would  have  to  be  between  6.8  and  8.4  lb/sec  for  gas  turbine  operated  at  half 
power,  and  between  11.0  and  14.5  lb/sec  for  gas  turbine  operated  at  full 
power.  The  overall  heat  transfer  characteristics  which  correspond  to  these 
operating  conditions  are  shown  in  Fig.  11.13.  There,  it  can  be  seen  that  the 
maximum  heat  transfer  rate  attainable  by  these  candidate  boilers  would  range 
between  9600  and  11500  Btu/sec  at  half  power  and  between  15400  and  19400 
Btu/sec  at  full  power  of  the  gas  turbine  engine  considered  (LM  2500  or  similar 
model). 
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The  gas  exit  temperature  and  the  gas-side  pressure  loss  characteristics 
for  the  waste-heat  boilers  are  shown  in  Figs.  11.14  aid  11.15.  It  is  seen  that 
the  lowest  gas  temperature  shown  in  Fig.  11.14  still  exceeds  300*F,  which 
implies  that  the  sulfuric  corrosion  would  not  occur.  Figure  11.15  shows  that 
when  these  candidates  waste-heat  boilers  were  integrated  with  gas  turbine  at 
50-percent  power,  the  gas  turbine  back  pressure  would  be  between  0.3  and  0.8 
psia,  indicating  that  those  tube  design  configurations  are  technically  acceptable 
from  the  turbomachinery  performance  view  point. 

II. 5  Baseline  Waste  Heat  Boiler 

The  background  information  which  was  used  to  select  a  baseline  design 
configuration  for  waste-heat  boiler  in  marine  propulsion  applications  are 
presented  in  Figs.  11.16  through  11.19.  Figure  11.16  summarizes  the  performance 
characteristics  of  four  candidate  tube  design  configurations  at  steam  temperature 
of  700*F  for  a  gas  turbine  operated  at  100-percent  power.  Similar  performance 
data  for  a  gas  turbine  operated  at  50-percent  power  are  presented  in  Fig.  11.17. 

Based  on  the  heat  transfer  rate  shown  on  the  far  right  frame  of 
Figs.  11.16  and  11.17,  Configuration  No.  1  would  yield  better  performance 
than  Configuration  No.  3,  No.  2  and  No.  4  in  that  order.  However,  the  gas-side 
pressure  losses  for  these  boiler  configurations  also  decreases  in  the  same 
order.  Unfortunately  the  higher  the  gas-side  pressure  loss,  the  higher  is  the 
back  pressure  to  the  gas  turbine,  and  according  to  the  correction  factor  for 
exhaust  pressure  loss  shown  in  Fig.  11.18  the  greater  will  be  the  loss  in 
turbine  power  output  . 

It  is  known  that  a  most  desirable  waste-heat  boiler  design  should  be  one 
which  can  provide  the  greatest  net  gain  in  power  output  when  coupled  with  a 
Rankine  cycle  power  conversion  system.  In  order  to  evaluate  the  net  gain 
in  power  output,  the  cycle  efficiency  of  the  Rankine  cycle  power  conversion 
system  must  be  identified.  From  the  results  of  a  waste-heat  recovery  system 
study  conducted  at  UTRC  (Ref.  2.9),  it  was  determined  that,  for  steam  condition 
of  700*F  and  3000  psia  and  a  condenser  pressure  of  3  inch  Hg,  the  cycle  efficiency 
of  a  typical  steam  Rankine  system  is  approximately  22  percent.  Therefore,  the 
net  gain  in  overall  power  system  is  equal  to  the  difference  between  the  Rankine 
cycle  power  output  and  the  loss  in  gas  turbine  output  power  due  to  the  increased 
back  pressure.  The  results  of  this  comparison  is  shown  in  Fig.  11.19. 

The  left  frame  of  Fig.  11.19  shows  the  net-gain  power  for  a  gas  turbine 
operated  at  100-percent  power.  It  was  found  that  Configurations  No.  2  and  No. 

4  would  provide  almost  equal  value  of  net-gain  power  and  this  gain  is  sub¬ 
stantially  higher  than  that  estimated  for  the  other  two  configurations.  In 
contrast,  the  right  frame  of  the  same  figure  shows  that  at  a  50-percent  (gas 
turbine)  power  condition,  Configuration  No.  3  is  the  most  desirable  selection. 
Because  improvement  in  the  propulsion  system  efficiency  at  cruise  conditions  is 
of  primary  concerns  in  naval  ship  operation,  Configuration  No.  3  was  selected 
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as  baseline  design  configuration  for  the  marine  gas  turbine  waste-heat  steam 
generator  application.  The  design  conditions  of  the  baseline  waste-heat  boiler 
can  now  be  defined  as  those  which  correspond  to  gas  turbines  operated  at  50-percent 
power,  and  off-design  conditions  are  defined  as  those  where  gas  turbines 
operate  at  any  power  level  other  than  the  50-percent  point. 

The  temperature  distribution  inside  the  baseline  waste  heat  boiler  operated 
at  its  design  condition  (as  detined  immediately  above)  is  shown  in  Fig.  11,20, 
and  the  overall  heat  transfer  coefficients  (defined  as  UA  in  Equation  6)  are 
shown  in  Fig.  11.21.  These  results  were  obtained  from  performance  calculations 
using  the  miniature  heat  exchanger  approach  and  the  nodal  system.  Similar 
information  for  off-design  condition  (corresponding  to  gas  turbine  100-percent 
power)  are  shown  in  Figs.  11.22  and  11.23.  From  the  temperature  distribution 
maps  (Figs,  11.20  and  11.22),  it  was  determined  that  this  waste-heat  boiler 
can  be  divided  into  three  distinctive  regions;  a  liquid-phase  heating  region,  a 
boiling-phase  region,  and  a  superheat ed-vapor  region.  In  liquid-phase  region, 
which  occupies  about  one-third  of  the  upper  portion  of  the  boiler,  the  feedwater 
is  heated  from  1 1 5 “ F  to  its  saturation  temperature  (which  is  417*F  for  steam 
pressure  of  300  psia).  In  the  boiling-phase  region,  which  is  shown  by  va  shaded 
boundary  and  occupies  approximately  60  percent  of  the  boiler  volume,  the 
teedwater  would  go  through  an  isothermal  boiling  process.  After  boiling,  the 
saturated  steam  would  be  superheated  in  the  last  10  percent  of  the  boiler 
whereupon  it  would  be  discharged  at  700*F. 

From  these  two  temperature  distribution  maps  in  Figs.  11.20  and  11.21,  it 
was  determined  that  the  temperature  gradient  between  any  two  nodes  is  less  than 
20'F  and  relatively  uniform  under  the  two  steady-state  operations.  Therefore, 
the  thermal  stress  concentration  should  not  be  a  problem.  However,  based  on 
the  performance  and  design  requirements,  the  water  velocity  wa3  calculated  to 
be  approximately  0.5  to  1.2  ft/sec  and  the  gas  velocity  was  approximately 
between  80  to  130  ft/sec.  The  transient  response  of  heat  transfer  characteristics 
and  thermal  stress  concentration  could  become  severe  under  the  dry-running  or 
changing  load  operations  and  this  should  be  investigated  before  the  final 
design  is  undertaken. 

The  product  of  the  overall  heat  transfer  coefficient  (U)  and  the  gas  side 
heat  transfer  areas  (Ag)  for  all  nodes  (miniature  heat  exchangers)  are  shown  in 
Figs,  11.21  and  11.23  for  baseline  waste  heat  boiler  operated  at  design  (gas 
turbine  50-percent  power)  and  off-design  (gas  turbine  100-percent  power) 
condition,  respectively.  These  values  were  computed  using  the  averaged  tem¬ 
peratures  shown  in  Figs.  11.20  and  11,22.  From  these  results,  it  was  determined 
that  the  value  of  UA  varied  between  1720  and  2585  Btu/hr-F  for  the  design 
condition  and  between  2100  and  3000  Btu/hr-F  for  the  off-design  operation.  It 
is  believed  that  those  values  are  mainly  determined  by  the  gas-side  film 
coefficient  rather  than  by  the  water-side  film  coefficient,  and  are  generally 
within  the  range  of  current  design  practices. 
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II. 6  Effect  of  Flow  Distribution  Control  on  Baseline  Waste-Heat  Boiler 

The  actual  flow  distribution  data  with  and  without  flow  distribution  control 
(see  Figs.  11.4  and  II. 5)  were  used  as  input  to  the  analytical  model  for 
analyzing  the  effect  of  flow  distribution-control  on  the  performance  of  the 
baseline  waste  heat  boiler.  The  results  of  this  analysis  are  presented  in 
Figs.  11.24  to  11.31. 

The  effects  of  flow  distribution  control  on  steam  temperature  of  the  baseline 
waste  heat  boiler  operated  at  design  (gas  tur'viut  50-percent  p<-> wpt)  and  off-design 
(gas  turbine  100-percent  power)  conditions  are  shown  in  Figs.  II  24  and  11.25, 
respectively.  It  should  be  noticed  that  only  those  water  flow  rates  which  can 
provide  superheated  steam  were  considered.  The  steam  temperatures  attainable 
for  the  uniform  flow  distribution,  actual  flow  distribution  with  and  without 
flow  distribution  controls  cases  are  shown  in  solid,  dashed,  and  semi-broken 
lines,  respectively.  These  results  show  that  if  the  baseline  waste  heat 
boiler  were  designed  for  uniform  flow  distribution  and  for  generated  steam 
700°F,  and  if  it  were  operated  with  actual  flow  distribution  without  flow 
distribution,  the  steam  temperature  would  decrease  to  approximately  450*F. 

However,  the  steam  temperature  could  be  maintained  at  700*F  if  the  water  flow 
rate  were  reduced  from  . .9  lb/sec  to  6.6  lb/sec.  If  this  occurred,  the  overall 
heat  transfer  rate  would  be  reduced  from  10167  to  8540  Btu/sec.  On  the  other 
hand,  from  Fig.  11.24,  it  can  be  seen  that  if  flow  distribution  control  were 
employed  and  if  the  water  flow  rate  were  maintained  at  its  design  condition 
(7.9  Ib/sec),  the  steam  temperature  would  increase  by  approximately  25*F. 
Alternately,  if  the  steam  temperature  were  to  be  maintained  at  the  design 
condition  (700*F),  the  water  flow  would  increase  to  8.1  lb/sec.  The  factors 
attributing  to  this  improvement  in  boiler  efficiency  are  believed  to  be  partly 
attributable  to  more  uniform  flow  distribution  and  partly  to  increased  gas  flow 
rate  in  boundary  layer  separation  control  (see  Fig.  11.3). 

The  heat  transfer  performance  characteristics  of  the  baseline  waste  heat 
boiler  operated  at  design  and  off-design  condition  for  cases  with  and  without 
flow  distribution  controls  are  shown  in  Fig.  11.26.  Again  the  solid  lines 
represent  the  results  of  the  assumed  uniform  flow  condition,  and  the  dash  lines 
and  semi-broken  lines  are  for  cases  with  and  without  flow  distribution  controls. 
The  asterisk  represents  the  conditions  where  the  700*F  steam  would  be  generated. 
One  can  readily  see  that  significant  improvement  in  boiler  efficiency  would  be 
expected  if  flow  distribution  controls  were  employed.  The  percentage  of  the 
boiler  efficiency  improvement  is  shown  in  Fig.  11.27  for  the  baseline  waste- 
heat  boiler  operated  at  its  design  condition. 

In  Fig.  11.27  which  illustrates  the  effect  of  flow  distribution  control  on 
the  performance  of  a  waste-heat  boiler,  the  overall  heat  transfer  rate  for  the 
assumed  uniform  flow  distribution  was  used  as  reference.  The  overall  heat 
transfer  rates  for  constant  water  (or  steam)  flow  rate  and  for  constant  steam 
temperature  were  obtained  from  Fig.  11.26.  In  Fig.  11.27  it  can  be  seen  that, 
without  flow  distribution  control  (blank  bars),  if  the  water  flow  rate  were  held 
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constant  at  7.9  lb/sec,  the  overall  heat  transfer  rate  would  be  reduced  by 
approximately  10  percent  which,  in  turn,  would  lower  the  steam  temperature  by 
approximately  250*F.  On  the  other  hand,  if  the  water  flow  rate  were  regulated 
so  as  to  maintain  a  constant  steam  temperature  of  700°F,  the  overall  heat 
transfer  rate  would  be  reduced  by  approximately  16  percent.  Although  this  loss 
is  greater  than  that  for  the  constant  flow  rate  case,  a  constant-steara-terapera- 
ture  operation  would  probably  be  more  suitable  for  Rankine  cycle  power  conver¬ 
sion  system  applications.  Similarly,  the  shaded  bars  shown  in  the  same  figure 
are  for  the  results  with  flow  distribution  control.  It  was  found  that  approx¬ 
imately  a  20  percent  improvement  in  boiler  efficiency  would  be  expected  for  the 
baseline  waste  heat  boiler  if  flow  distribution  control  were  considered. 

The  predicted  temperature  distribution  for  the  baseline  waste-heat  boiler 
operating  at  design  conditions  and  based  on  the  actual  velocity  data  with  and 
without  flow  distribution  controls  are  shown  in  Figs.  11.28  and  11,29,  respectively; 
the  corresponding  overall  heat  transfer  coefficients  for  these  two  flow  cases  are 
shown  in  Figs.  11.30  and  11.31.  From  Figs.  11.28  and  11.29,  it  is  seen  that  the 
gas  temperature  near  the  exit  becomes  highly  nonuniform,  and  that  this  nonuniformity 
is  more  apparent  in  the  absence  of  flow  distribution  controls.  This  temperature 
distribution  nonuniformity  is  certain  to  cause  uneven  thermal  expansion  in  the 
heat-exchanger  tubes,  which  in  turn  may  result  in  thermal  stress  concentration 
problems.  It  is  also  seen  that  due  to  the  no-flow  condition  occurring  in  the 
far-right  column  (uncontrol  case),  the  gas  temperatures  and  the  liquid  temperature 
may  be  in  equilibrium  (see  Fig.  11.28),  indicating  a  region  of  "no-heat-transfer" 
within  the  waste-heat  boiler  (see  Fig.  11.30),  The  performance  degradation  of  the 
waste-heat  boiler  can  then  be  expected  when  no  proper  flow  distribution  control 
methods  are  implemented. 

The  effects  of  gas  inlet  temperature  on  the  baseline  waste  heat  boiler 
performance  were  also  investigated  and  the  results  obtained  are  shown  in  Fig. 

11.32.  From  these  results,  it  can  be  seen  that  change  in  gas  inlet  temperature 
would  have  more  profound  effect  on  the  baseline  waste  heat  boiler  performance 
under  constant  steam  temperature  operation  than  under  constant  steam  flow  rate 
operation. 
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TABLE  II. 1 

(A)  Constant  for  Equation  (la),  (lb),  and  (lc) 


Laminar  Flew  ; 

C1 

-  0.595 

c2 

■  0.498  C3  ■  0.140 

Turbulent  Flow  ; 

C1 

*  0.0046 

C2 

-  0.927  C3  -  0.628 

Supercritical  Flow; 

C1 

-  0.003354 

c2 

-  0.951  .  C3  -  0.435 

C4  -  0.38 

IB) 

Constant  for  Equations 

(4a),  (4b), 

(4c) , 

,  (5a),  and  (5b) 

Turbulent-Turbulent 
NRef>2000,  NRe  >2000 

— - — - - - P -  - 

Viscous-Turbulent 

NRef<1000,  NRe  >200 

-  ■  -  " _ _ _ _ 

Turbulent -Viscous 
NRef>2000,  NReR<1000  NRe 

Viscous-Viscous 
f<1000,  NRe  <1000 

m 

0.25 

0.25 

1.0 

1.0 

n 

0.25 

1.0 

0.25 

1.0 

0.079 

16.0 

0.079 

16.0 

c. 

0.079 

0.079 

16.0 

16.0 

N 

4.0 

3.50 

3.50 

2.75* 
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FIG.  H. 


ACTUAL  FLOW  DISTRIBUTION  WITH  FLOW  DISTRIBUTION  CONTROL 

•  GAS  TURBINE  50%  POWER 


•1- 
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FIG.  U4 


GAS  FLOW  INLET  CONDITIONS  FOR  MARINE  GAS  TURBINE 
WASTE  HEAT  BOILER  -  GAS  TURBINE  100%  POWER 

AX  =  WIDTH  OF  THE  GAS  FLOW  PATH,  (in  > 

M  =  AVERAGED  FLOW  RATE,  (ib/sec) 

P  =  AVERAGED  STATIC  PRESSURE,  (psial 
T  =  AVERAGED  STATIC  TEMPERATURE.  <°F> 


AVERAGED  FLOW  PARAMETERS 
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24 

24 

24 

24 

24 

— 

• 

M 

32 

32 

32 

32 

32 

— 

T 

856 

856 

856 

856 

856 

P 

14.94 

14.94 

14.94 

14.94 

14.94 

— 

actual  flow 

DISTRIBUTION 

WITHOUT 

CONTROL 

AX 

24 

24 

24 

24 

24 

— 

• 

M 

80.3 

50.6 

15.5 

13.6 
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856 

P 

14,96 
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14,89 

14,89 

14.89 
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FIG  E5 


GAS  FLOW  INLET  CONDITIONS  FOR  MARINE  GAS  TURBINE 
WASTE  HEAT  BOILER  -  GAS  TURBINE  50%  POWER 

AX  -  WIDTH  OF  THE  GAS  FLOW  path  (in  I 
M  •  AVERAGED  FLOW  RATE,  Hb'sec) 

P  i  AVERAGED  STATIC  PRESSURE.  (psi;il 
T  -  AVERAGED  STATIC  TEMPERATURE  I°F| 
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FIG.  E6 


CONFIGURATIONS  OF  CIRCULAR  FINNED  TUBES  FOR  MARINE  WASTE-HEAT  BOILER 

APPLICATION 
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ARRANGEMENT  FOR  MARINE  GAS  TURBINE  WASTE  HEAT  STEAM  GENERATOR  APPLICATIONS 


TUBE  O  D  =  1  024  m  .  9  FINS/m 

FIN  LENGTH  =  0  356  in  FIN  THICKNESS  =  0  01 2  in 


100%  POWER 


SATURATION  TEMP  =  417*F 


TRANSFER  RATE.  Btu/sec 
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FIG.  1.13 


LOSS,  psi 
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FIG.  H15 


COMPARISON  OF  GAS-SIDE  PRESSURE  LOSS  FOR  WASTE  HEAT  BOILERS  OPERATED  AT 

DESIGN  AND  OFF-DESIGN  CONDITIONS 

•  UNIFORM  FLOW  AND  EFFECTIVE  TUBE  LENGTH  =  200  It 


WATER  FLOW  RATE,  Ib/sec 


PERFORMANCE  CHARACTERISTICS  OF  WASTE  HEAT  BOILER  AT  CONSTANT  STEAM  TEMPERATURE 


WATER  FLOW  RATE  GAS  EXIT  TEMP.  GAS-SIDE  PRESS.  LOSS  HEAT  TRANSF  RATE 

(Ib/sec)  (°F)  (psi)  (1000  Btu/sec) 
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FIG.IT.21 


DISTRIBUTION  OK  OVERALL  HEAT  TRANSFER  COEFFICIENT  FOR  BASELINE  WASTE  HEAT 
BOILER  OPERATED  AT  DESIGN  CONDITION  (GAS  TURBINE  50%  POWER,  UNIFORM  FLOW) 

•  HEAT  TRANSFER  COEFFICIENT  IN  Btu/ni-F 
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R82-955750-4  FIG  U  23 

DISTRIBUTION  OF  OVERALL  HEAT  TRANSFER  COEFFICIENT  FOR  BASELINE  WASTE  HEAT 
BOILER  OPERATED  AT  OFF-DESIGN  CONDITION  (GAS  TURBINE  100%  POWER, 

UNIFORM  FLOW) 

•  HEAT  TRANSFER  COEFFICIENT  IN  Bln  hr-F 
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STEAM  TEMPERATURE. 
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FIG.  2  25 


EFFECT  OF  FLOW  DISTRIBUTION  CONTROL  ON  STEAM  TEMPERATURE  OF  BASELINE 
WASTE  HEAT  BOILER  AT  OFF-DESIGN  CONDITION 


- - ASSUMF  D  UNIFORM  f  LOW 
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_ actual  flow  with  control 
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OFF-DESIGN  PERFORMANCE  CHARACTERISTICS  OF  GAS  TURBINE 
WASTE-HEAT  STEAM  GENERATOR 

*  STEAM  TEMPERATURE  -  700  »F 

— _  ASSUMED  UNIFORM  FLOW  DISTRIBUTION 
«___  ___  ACTUAL  FLOW  DISIHI8UTION  WITHOUT  CONTROL 
ACTUAL  FLOW  DISIMIBUTION  WITH  CONTROL 
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RATIO  OF  OVERALL  HEAT  TRANSFER  RATE.  Q/Q,. 
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FIG.  E27 


EFFECT  OF  FLOW  DISTRIBUTION  CONTROL  ON  WASTE  HEAT  RECOVERY 
STEAM  GENERATOR  PERFORMANCE 
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FIG  11.29 


DISTRIBUTION  OF  AVERAGED  TEMPERATURE  FOR  BASELINE  WASTE  HEAT  BOILER  OPERATED 
AT  DESIGN  CONDITION  WITH  FLOW  DISTRIBUTION  CONTROL 


T,  =  LIQUID  TEMPERATURE.  *F 
Tg  =  GAS  TEMPERATURE.  *F 
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FIG  11.31 


DISTRIBUTION  OF  OVERALL  HEAT  TRANSFER  COEFFICIENT  FOR  BASELINE  WASTE  HEAT 
BOILER  OPERATED  AT  DESIGN  CONDITION  WITH  FLOW  DISTRIBUTION  CONTROL 


•  HEAT  TRANSFER  COEFFICIENT  (UA1  IN  Blu/hr-F 
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SECTION  III 

FORMU1.AT ION  OF  EXPERIMENTAL  PROORAM  AND  SCHEDULE 


The  effect  of  flow  distribution  control  on  a  narine  gas  turbine  waste- 
heat  boiler  has  been  analyzed,  and  the  results  of  this  analysis  are  presented 
in  Section  II.  It  is  obvious  that  a  complex  system  like  a  waste-heat  boiler 
cannot  be  designed  and  operated  successfully  without  first  conducting  a  carefully 
designed  experimental  program  to  generate  sufficient  technical  information  and 
operational  experience.  In  this  section,  an  experimental  program  formulated 
for  the  candidate  waste-heat  boiler  is  presented.  The  technical  information 
and  operational  experience  desired  from  this  experimental  program,  the  experiment 
program  plan,  and  the  overall  program  schedule  and  effort  are  discussed  in 
detail. 


III. I  Technical  Information  and  Operational  Experience 
Desired  from  the  Experimental  Program 

As  shown  in  Table  III.l  the  objectives  of  this  suggested  experimental 
program  can  be  divided  into  two  categories:  (1)  that  which  would  provide 
sufficient  technical  information  so  comparisons  with  the  analytical  results  can 
be  made,  and  (2)  that  which  w?ulu  allow  operational  experience  to  be  gained 
with  the  use  of  gas  turbine  waste-heat  recovery  propulsion  systems.  In  order 
to  obtain  this  technical  information  and  operational  experience,  an  experimental 
program  must  first  provide  for  the  preparation  of  components  and  instrumentation 
necessary  for  testing.  Such  preparation  would  include  design  specification  of 
the  component  and  instrument  designs,  then  fabrication,  preliminary  demonstra¬ 
tions,  and  testing  apparatus  check-out. 

A  demonstration  of  steady  state  operation  of  an  experiment  model  would  seem 
the  first  step  in  the  experiment  program  following  the  preparation  tasks. 
However,  before  the  experimental  model  can  be  operated  at  its  steady  state 
condition,  the  scart-up  and  shut-down  procedures  would  have  to  be  specified  and 
evaluated.  Then  che  technical  information  and  operational  experience  expected 
to  be  gained  from  this  steady  state  operation  can  he  obtained  from  flow  visual¬ 
ization,  and  temperature  and  pressure  measurements.  Once  Che  steady  state 
demonstration  i3  completed,  the  transient  characteristics  of  the  model  can  be 
demonstrated,  particularly  as  they  relate  to  naval  ship  propulsion  applications. 
The  technical  information  and  operational  experience  expected  to  be  gained  from 
this  transient  (dynamic)  operation  demonstration  should  include  those  related 
to  flow  and  temperatue  stability,  thermal  performance  response,  and  possibly 
thermal  stress  concentration  problems. 
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One  major  consideration  in  the  design  of  a  waste-heat  steam  generator  for 
naval  ship  propulsion  application  is  the  characteristics  and  limitation  of 
the  additional  shaft  power  obtainable  from  waste-heat  recovery  at  different 
part-power  operating  condition  of  the  gas  turbine  engine  in  order  to  match  the 
duty-cycle  operational  (i,e.  speed  and  time)  requirements.  The  suggested 
experiment  program  should  identify,  at  least  qualitatively,  the  nature  of  these 
cnaracterist ics  and  limitations,  if  not  specifically  by  estimating  their 
magnitude.  The  effect  of  flow  distribution  control  on  power  output  limitation 
must  also  be  estimated  experimentally  so  that  the  results  obtained  can  be 
compared  with  the  analytical  results  of  Section  II  in  this  report. 

To  assess  control  methodology,  the  experimental  program  should  also  provide 
sufficient  information  relating  to  the  use  of  either  pneumatic  or  electronic 
controls  to  regulate  the  performance  of  the  waste-heat  boiler  system.  This 
control  system  should  be  able  to  regulate  the  boiler  pressure  and  feedwater 
flow  rate  so  the  temperature  and  flow  rate  of  the  gas  turbine  exhaust  can  be 
matched  with  the  specified  steam  outlet  temperatures.  Since  the  waste-heat  boiler 
will  be  required  to  operate  under  dry-running  condition  for  self-cleaning  purposes, 
the  control  system  required  to  cope  with  the  lost-of-coolant  problem  may  not  be 
critical,  however  appropriate  devices  to  prevent  such  operating  condition  will  be 
required.  ' 

The  procedure  to  assess  the  dry-running  operation  of  a  marine  gas  turbine 
waste-heat  boiler  has  yet  to  be  established,  although  many  of  the  routine 
operating  procedures  for  marine  waste-heat  boilers  may  be  similar  to  those  for 
landbased  combined  cycle  system  operation.  However,  for  the  dry-running 
operation,  special  consideration  must  be  given  to  the  rate  at  which  the  feed 
water  is  drained  and  recharged  to  avoid  undesirable  damage  or  deterioration  of 
the  boiler  lubes.  This  experimental  program  should  produce  valuable  information 
to  assist  in  establishing  the  dry-running  procedure. 

The  U.S.  Navy  is  expected  to  be  quite  interested  in  identifying  the  manpower 
requirements  for  a  marine  waste-heat  recovery  propulsion  system.  Again,  the 
experimental  program  should  clarify  this  question  through  an  assessment  of 
system  maintainability  and  reliability.  Areas  of  manpower  needs  should  be 
quite  similar  to  those  of  landbased  waste-heat  recovery  plants,  while  the 
number  of  men  needed  in  a  marine  propulsion  application  could  be  reduced  if  the 
waste-heat  boiler  were  designed  with  less  complexity  and  higher  reliability. 

The  importance  of  gaining  as  much  data  and  operating  experience  as  possible 
cannot  be  overemphasized  whan  dealing  with  the  maintainability,  reliability  and 
safety  characteristics  which  are  common  concerns  for  any  new  system. 
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III. 2  Experiment  Program  Plan 

An  experimental  program  plan  for  the  flow  distribution  control  study  of 
marine  gas  turbine  waste  heat  boiler  has  been  prepared.  This  program  consists 
of  four  major  tasks  as  shown  in  Table  III. 2,  The  first  task  involves  the 
design  and  fabrication  of  the  experiment  model,  while  the  second  task  is 
devoted  to  the  set-up  of  experimental  apparatus  including  the  acquisition  of 
necessary  auxiliary  components.  Tasks  3  and  4  are  directed  coward  conducting 
the  actual  experiment,  including  data  recording  and  post  test  evaluation. 

1 1 1 . 2  . 1  Design  and  Fabrication  of  Experiment  Model 

The  first  step  in  the  design  and  fabrication  of  the  experimental  model 
is  to  determine  the  model  size.  Because  the  cost  of  building  a  full-scale  test 
model  as  well  as  the  heat  source  required  for  the  experiment  would  be  enormous, 
a  one-fifth  scale  model  is  suggested.  It  is  believed  that  this  scaled  model 
can  be  designed  and  fabricated  in  a  reasonable  time  frame  and  at  an  acceptable 
cost  that  would  provide  these  desired  information  described  in  Table  III.  1 . 

The  thermal  condition  (flow  lace,  temperature,  and  pressures)  of  the  working 
fluid  are  usually  determined  from  the  availability  of  the  test  (auxiliary) 
equipment  and  by  using  the  principles  of  similitude.  In  principle,  the  test 
model  and  the  full-scale  unit  should  have  the  same  Nusselt,  Prandtl,  Reynolds, 
and  Mach  numbers.  The  flow  passage  in  the  test  model  and  that  of  the  full-scale 
unit  should  also  be  geometrical lv  similar.  Because  the  detailed  temperature  and 
flow  distributions  are  the  primary  concerns  in  the  present  study,  it  is  more 
desirable  to  use  larger  flow  passages  with  fewer  number  of  tubes  in  the  model. 

As  long  as  the  flow  conditions  are  based  on  the  •  rinciples  of  similitude,  the 
pressure  loss  and  heat  transfer  characteristics  in  the  test  model  should  differ 
little  from  the  full-scale  unit. 

The  instrumentation  needed  for  the  present  study  must  be  capable  of  measuring 
temperature,  pressure,  and  flow  conditions.  In  a  conventional  heat  exchanger 
experiment,  as  few  as  four  temperature  measurements  might  suffice  (that  is,  the 
inlet  and  outlet  temperatures  for  the  hot  and  cold  fluids).  However,  in  the 
suggested  program  where  «.he  boiler  is  a  single-pass  counter  crossflow  heat 
exchanger,  the  outlet  temperatures  will  not  be  uniform.  Accordingly,  at  least 
fifteen  temperature  measurements  would  be  required  to  identify  the  inlet  and 
exit  flow  conditions.  If  internal  temperature  distributions  are  also  to  be 
measured  (as  shown  in  Figs.  11.20  and  11.22),  additional  temperature  probes 
would  be  needed. 

The  magnitude  cf  pressure  drop  across  the  heat-transfer  matrix  particularly 
on  the  gaa  side  of  the  waste-heat  boiler  is  as  important  as  the  heat-transfer 
performance.  The  test  rig  may  be  designed  so  the  duct  cross  section  is  the 
same  as  that  of  the  inlet  face  of  the  heat-transfer  matrix  under  test,  in  which 
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case  simple  scatic  pressures  in  the  duct  may  be  satisfactory.  If  this  is  not 
practical,  allowances  should  be  made  for  differences  in  the  kinetic  pressure 
head  which  changes  with  flow  passage  size.  It  is  important  that  at  least  ten 
diameters  of  straight  duct  preceed  the  heat-transfer  matrix  to  assume  a  uniform 
velocity  distribution  across  the  face  of  the  duct.  Because  pressure  drop  data 
are  important,  it  is  more  desirable  to  use  the  piezometer  ring. 

The  simplest  and  most  accurate  means  of  measuring  the  gas  flow  for  this 
experimental  program  would  be  to  use  a  flow  nozzle  mounted  at  the  air  inlet  and 
a  draft  fan  be  mounted  on  the  outlet  side  of  the  heat-transfer  matrix.  This 
would  preclude  errors  stemming  from  turbulence  and  poor  velocity  distribution 
in  the  flow-rate  measurements.  To  measure  the  flow  distribution  in  each  gas 
path  Isee  Fig.  11.20),  pitot  static  tubes  can  be  used.  At  least  two  pitot 
static  tubes  must  be  used  for  each  gas  path  and  each  tube  should  be  installed 
downstream  of  the  heat  transfer  matrix  in  order  to  avoid  disturbing  the  flow 
fteld  in  the  test  section.  Finally,  care  must  be  exercised  to  minimize  the 
t  low  leakage,  heat  loss,  and  boundary  effect  during  the  design  and  fabrication 
ot  the  experiment  model  as  severe  flow  leakage  and  heat  loss  could  cause 
difficulties  in  the  analysis  of  test  results. 

It l . 2 . 2  Auxiliary  Equipment  Set-Up 

The  auxiliary  equipment  needed  for  this  experiment  will  include:  a  hot-air 
lor  gas)  supply  and  discharge  system;  a  pressurised  feedwater  and  stean  handling 
system;  a  control  device  to  regulate  the  flow  rate,  flow  distribution,  temperature, 
and  pressure  for  both  working  fluids;  and  data  acquisition  and  recording 
devices.  One  possible  arrangement  of  the  experimental  apparatus  is  shown  in 
Fig.  lll.l. 

The  nol -air /hot -gas  supply  and  discharge  system  would  require  a  draft  fan  which 
should  be  mounted  at  the  downstream  of  the  test  model,  and  a  combustor  which 
burns  either  natural  gas  or  propane  for  generating  hot  gas  needed  to  simulate 
the  gas  turbine  exhaust.  Both  draft  tan  and  combustor  could  be  controlled  from 
a  central  control  box  to  obtain  desirable  gas  flow  rates  and  gas  temperature. 

The  pressurized  feedwater  aupply  system  would  consist  of  a  pressurized  feedwater 
manifold,  a  flow  rate  regulator,  a  sterna  manifold,  a  radiator,  a  condensate 
tank,  and  a  pump.  Both  the  flow  rate  regulator  and  the  pump  would  be  controlled 
by  the  central  control  device.  A  baffle  plate  should  be  installed  in  the 
diffuser  to  regulate  the  desirable  gas  flow  distribution.  All  the  test  data 
including  temperature,  pressure,  flow  rates  would  be  recorded  by  means  of  an 
automatic  data  acquisition  and  recording  device  for  later  evaluation  and 
snalysis. 
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1 1 . 1 . 3  Test  Procedure 

The  test  procedures  to  be  conducted  in  this  suggested  exper iment a’  program 
must  include,  as  a  minimum,  all  items  shown  in  category  C  of  Table  111.2.  The 
first  step  in  conducting  the  experiment  would  be  to  assemble  and  check-out  the 
major  components  according  to  the  layout  drawing.  Shakedown  testing  would  then 
follow  to  demonstrate  the  functional  capability  and  strctural  integrity  of  the 
experiment  apparatus.  Some  minor  modi f icat ions  and  adjustment  might  be  necessary 
in  the  early  stage  of  the  experiment  before  substantive  test  program  can  be 
commenced.  Piping  and  wiring  details  would  be  determined  prior  to  the  assembly 
of  rhe  experimental  unit,  and  instrumentation  and  controls  of  the  flow  rate, 
temperature,  and  pressure  would  be  installed  and  calibrated  before  the  actual 
experiments  were  performed. 

An  extremely  important  aspect  of  the  suggested  program  is  flow  visualization 
test  in  the  flow  distribution  control  study  of  gas  turbine  waste-heat  recovery 
steam  generator.  A*  cited  in  Phase-I  study  (Section  11.2  of  Ref.  11.3)  the 
actual  flow  distribution  at  the  exit  of  the  gas  turbine  exhaust  is  highly 
irregular  and  nonuniform.  The  major  portion  of  the  flow  was  found  to  be  near 
the  rear  section  of  the  elbow  and  some  reversed  flow  as  observed  in  the  regions 
near  the  front  section.  To  complicate  matters,  these  flow  distributions  are 
actually  three-dimensional.  In  order  to  gain  an  insight  into  flow  distribution 
nonuniformity  in  the  waste-heat  boiler  performance,  two-dimensional  flow 
distributions  were  assumed  in  the  analytical  study.  However,  the  results  of 
such  an  analytical  study  can  only  be  compared  with  those  of  the  experimental 
study  on  the  same  basis,  i.e.  of  a  two-dimensional  flow  experiment.  Therefore, 
for  naval  applicators,  the  experiments  must  also  include  three-dimensional  flow 
distributions,  since  the  results  obtained  from  three-dimensional  flow  testing 
would  be  beneficial  in  any  modifications  of  the  analytical  model,  which  is 
deemed  necessary. 

Flow  visualization  for  the  hot  combustion  gss  can  b _■  conducted  by  attaching 
tufts  of  thread  or  yarn  to  the  passage  walls,  or  by  attaching  these  tufts  to  a 
wire  probe  that  can  he  moved  about  in  the  flow  field.  Smoke  can  alco  be 
employed,  but  its  use  is  usually  not  very  satisfactory  because  the  smoke 
filament  tends  to  be  dispersed  so  rapidly  by  turbulence  that  the  technique  is 
applicable  only  lor  relatively  low  Reynolds  numbers  and  simple  geometries.  It 
should  be  obvious  that  flow  visualization  tests  can  be  conducted  more  conveniently 
if  the  models  are  made  of  a  transparent  plastic,  such  as  Lucite. 

When  preparing  for  heat  transfer  performance  testing,  particularly  with  nonunifonn 
flow  distributions,  consideration  must  be  given  to  the  flow  stebility  to  aesure 
that  the  teat  data  are  coneietent  and  repeatable.  During  the  literature  eurvey 
conducted  in  Phase-1  study,  it  wae  noticed  that  the  flow  in  the  transitory 
stall  region  of  a  two-dimensional  diffuser  is  inherently  unstable,  and  any 
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disturbance  could  shift  the  stall  region  from  one  wall  to  the  other.  Therefore, 
a  flow  stability  test  must  be  conducted  in  conjunction  with  flow  visualization., 

A  simple  way  of  assuring  tnat  fluid  flow  is  stable  is  to  take  periodic  flow 
measurements  (five  to  ten-minutes  intervals  are  suggested)  at  each  fixed 
operating  condition.  This  procedure  should  be  continued  until  three  successive 
readings  of  flow  measurement  show  negligible  change. 

The  last  item  in  the  suggested  test  procedure  of  Table  111.2  is  the  heat 
transfer  performance  test  which  should  consist  of  a  stead-state  operation  test 
and  a  transient  (dynamic)  operation  test.  In  the  steady-state  operation,  the 
thermal  output  characteristics  and  limitation  of  the  waste  heat  boiler  under 
various  flow  conditions  should  be  determined.  The  effect  of  variations  in  flow 
rate  and/or  pressure  of  feedwater  must  be  assessed,  and  the  magnitude  of  such 
parameters  as  critical  temperature,  pressures,  stresses,  control  feedbacks  must 
be  determined  and  examined  as  well. 

The  transient  operation  test  would  investigate  the  dynamic  characteristics 
of  the  marine  gas  turbine  waste  heat  boiler  during  the  dry  running  and  off-design 
operations.  The  first  dynamic  test  should  assess  the  effect  of  the  heat  input 
which  varies  according  to  the  duty  cycle  operations.  There  are  three  possible 
operation  modes:  (1)  variable  gas  temperature  with  constant  flow  rate;  (2) 
variable  gas  flow  rate  with  constant  gas  temperature;  and  (3)  both  variable  gas 
temperature  and  flow  rate.  The  response  time  and  stability  of  heat  exchanger 
performance  under  these  dynamic  teste  should  be  determined  and  examined.  If 
the  steam  temperature  were  maintained  constant,  the  control  procedure  and 
control  requirements  would  have  to  be  identified.  The  «o3t  important  information 
to  be  acquired  from  a  dry  running  cest  would  be  the  rate  of  draining  and 
recharging  the  feedwater.  In  addition,  following  dry  running,  the  metal 
temperature  would  be  approximately  900  to  1000*F,  and  therefore  the  requirement 
of  cool-down  process  must  be  determined.  Finally,  the  possibility  of  wef. 
running  with  static  flow  condition  should  also  be  explored  during  the  transient. 

III. 2.4  Foot  Teat  Evaluation 

The  last  task  of  the  suggested  experimental  program  plan  should  be  the  post 
experiment  inspection  and  evaluation  of  the  test  results.  This  task  would  seek 
to  assess  the  results  of  the  experiment  in  terms  of  establishing  requirements 
for  future  modifications  to  the  waste  heat  boiler  design  and  analysis. 

The  first  phase  of  this  last  task  would  be  routine  inspection  and  examination 
of  the  test  model  to  determine  its  general  condition  and  to  estimate  whether 
any  degradation  in  conditions  may  have  affected  the  test  results.  Evaluation 
of  the  test  results  would  then  be  made  with  particular  emphasis  being  directed 
toward  the  technical  information  and  operational  experiences  desired  (refer  to 
the  discussion  in  Section  111.1).  Recommendations  regarding  to  the  future 
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design  of  marine  gas  turbine  waste  heat  recovery  propulsion  systems  should  be 
included  as  part  of  this  last  task;  methods  directed  at  removing  operational 
limitations  should  be  explored;  and  reliability  and  maintenance  requirements 
should  also  be  assessed. 

III. 3  Overall  Program  Schedule  and  Effort 

The  overall  program  schedule  formulated  and  man-hour  effort  estimated  to 
conduct  this  suggested  experimental  study  are  shown  in  Fig.  III. 2.  Although 
overlaps  in  program  schedule  for  certain  activities  are  necessary  because  of 
the  nature  of  a  particular  test  or  because  of  the  need  to  shorten  the  performance 
period,  it  can  be  seen  that  the  activities  described  are  generally  consistent 
with  the  program  plan  discussed  in  Section  III. 2. 

The  longest  period  required  for  this  suggested  experimental  study  would  be 
those  for  model  preparation  (activities  No.  1  and  No.  2)  and  heat  transfer 
performance  testing  (activity  No.  8);  these  would  require  approximately  five 
and  four  months,  respectively.  The  experiment  apparatus  setup  time  (including 
that  for  the  acquisition  of  control  devices,  data  recording  system,  pumps, 
blower,  burner  test  site  preparation,  utility  hook-up,  piping,  and  wiring)  as 
well  as  the  shakedown  would  require  approximately  three  months.  Collectively, 
the  flow  visualization  and  flow  stability  tests  would  require  approximately 
two  months,  and  finally,  two  months  would  be  needed  for  post  test  evaluations, 
analysis  of  test  results,  recommendation,  and  reports  preparations. 

The  last  column  of  Fig.  111.2  shows  the  man  hours  estimated  to  completed 
each  task.  It  can  be  seen  that  most  of  the  engineering  effort  would  be  spent 
in  model  design  definition  and  post  test  evaluation  while  the  technician's 
time  would  be  directed  toward  model  fabrication,  test  rig  setup,  testing,  and 
data  recording.  It  should  also  be  emphasized  that  a  joint  effort  from  engineer 
and  technician  is  also  necessary  during  each  task  to  assure  the  success  of  the 
experimental  program.  Therefore,  for  the  entire  suggested  program  the  total 
engineering  time  is  estimated  to  be  approximately  1000  hours  and  that  for 
technician,  approximately  2050  hour#. 
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TABLE  III. 1 

Technical  Information  and  Operational  Experience 
Wanted  From  Experimental  Program 


A,  Related  to  Flow  Distribution  Control  Study 

,  Design  Specification  of  Test  Apparatus 
Steady-State  Operation  Demonstration 
.  Transient  (Dynamic)  Operation  Demonstration 
.  Output  Characteristics  and  Limitation  Identification 

B.  Related  to  Naval  Propulsion  System  Applications 

.  Control  Characteristics  Assessment 

Dry  Running  and  Duty  Cycle  Operating  Procedure  Assessment 
,  Maintainability,  Reliability,  and  Safety  Assessment 
.  Operational  Manpower  Requirement  Assessment 
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TABLE  II I. 2 

Experiment  Program  Plan 

A.  Design  and  Fabrication  of  Experiment  Model 

.  Size  of  experiment  model 

Test  condition  and  consideration 
Pr  inc  iples.  ot  similitude 
Adequate  and  potent  instrumentation 
.  Leakage,  heat  loss,  and  boundary  effects 

H  Auxiliary  Equipment  Setup 

.  Hot  air/gas  supply  and  discharge  equipments 
.  Water  supply  and  steam  handling  equipments 
.  Control  devices 

.  Data  acquisition  and  recording  devices 
0.  Test  Procedure 

.  Shake-down  flow  and  structural  tests 
.  Calibration  of  instrumentation 
.  Flow  visualization  test 
.  Flow  stability  test 
,  Heat  transfer  performance  test 

l).  Post  Test  Evaluation 
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POTENTIAL  LAYOUT  OF  EXPERIMENT  APPARATUS  FOR  FLOW  DISTRIBUTION 
CONTROL  STUDY  OF  MARINE  WASTE  HEAT  STEAM  GENERATOR 
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APPENDIX  A 

DESCRIPTIONS  OF  WASTE-HEAT  BOILER  COMPUTER  PROGRAM 


The  heat  exchanger  computer  program  which  was  used  to  study  the  effect  of 
flow  distribution  control  on  marine  waste-heat  steam  generator  performance  was 
developed  from  the  analytical  model  presented  in  Section  II  of  this  report. 

This  program  can  be  used  to  predict  the  overall  performance,  size,  and  manu¬ 
facturing  cost  for  many  types  of  crossflow  heat  exchangers.  It  is  applicable 
to  nearly  any  kind  of  gas  and  liquid  provided  that  their  heat  transfer  and 
pressure  loss  correlations  are  expressed  in  the  form  shown  in  the  book  "Compact 
Heat  Exchangers"  by  Kays  and  London,  The  surfaces  of  the  heat  exchanger  core 
can  be  of  plate-fin,  finned-tube,  or  screen-matrices  geometries.  Although  the 
program  can  be  easily  extended  to  other  heat  exchanger  applications,  it  is 
currently  limited  to  the  cases  where  the  gases  are  flowing  on  the  shell  side 
and  the  liquid  is  on  the  tube  side.  Tne  liquid  may  undergo  phase  changes 
(from  liquid  phase  to  boiling  phase,  and  then  to  superheated  vapor,  but  not  in 
reverse  process)  depending  on  the  design  requirement  specified. 


Program  Structure 

The  program  was  organized  in  hierarchical  structure  as  shown  in  Fig.  A.l. 

The  main  program  which  is  called  HXMAIN  is  the  commanding  portion  used  to  call  the 
three  subroutines  (HXINPUT,  HXCALC,  and  KXOUTPUT)  which  perform  the  specific  tasks 
as  indicated  in  their  respective  boxes  in  the  figure.  The  subroutine  HXINPUT 
reads  the  input  data,  interprets  and  initializes  these  data  prior  to  the  heat 
exchanger  performance  calculations,  and  finally  stores  all  of  the  relevant  informa 
tion  in  the  common  blocks.  The  input  data  consists  of  a  job  title,  job  control 
parameters,  the  inlet  flow  conditions,  the  interconnection  of  flow  paths,  flow 
properties,  and  heat  transfer  and  pressure  loss  correlations,  ail  of  which  are 
explained  in  the  next  section. 

The  subroutine  HXCALC  is  the  calculation  section  of  the  heat  exchanger  progr* 
which  is  based  on  the  flow  equations  and  the  computation  process  described  in 
Section  II  of  this  report.  This  program  takes  the  input  data  from  the  common 
blocks  as  needed  during  the  computation  process  and  also  stores  the  computed 
results  in  the  common  blocks.  The  function  of  the  HXOUTPUT  routine  is  to 
translate  the  results  of  the  heat  exchanger  performance  computed  by  the  HXCALC 
routine  into  practical  engineering  units  and  print  them  as  hard-copy.  The 
output  results  consist  of  convergence  information,  summary  results,  temperature 
distributions  for  gss,  liquid,  snd  tube  walls,  and  heat  transfer  coefficients 
throughout  the  entire  hest  exchanger  core. 
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The  program  which  is  written  in  Fortran  language  is  implemented  on  a  UNIVAC 
1110/80  computer  system  and  requires  approximately  50K  core  storage.  The 
computational  time  required  for  a  typical  study  varies  according  to  the  tolerance 
of  convergence  and  the  number  of  nodes  specified  for  the  heat  exchanger  core. 

For  the  cases  studied,  each  of  which  consist  of  five  gas  paths  and  one  liquid 
patii,  one  hundred  nodes  in  heat  exchanger  core,  and  a  f ive-degree-fahrenheit 
tolerance  on  the  exit  temperatures  for  both  working  fluids,  the  computational 
time  was  approximately  30  to  40  seconds. 

In  the  following  section,  the  input  parameters  are  discussed,  the  results 
of  a  example  run  are  presented,  and  the  listing  of  the  Fortran  statements 
for  the  main  program  and  three  major  subroutines  are  given. 


Description  of  Input  Data 

The  input  data  for  the  Waste-Heat  Boiler  Computer  Program  are  listed  in 
Tables  A.l  and  A. 2.  Table  A.l  contains  four  different  types  of  input  data;  the 
job  title,  the  job  control  parameters,  the  inlet  flow  condit ions ,  and  the  nodal 
connection  method  for  each  flow  path. 

The  job  title,  which  may  be  comprised  of  up  to  72  character ist ics ,  must  be 
punched  on  a  BCD  card.  The  job  control  parameters  (there  are  twenty  two  of 
inem)  are  defined  as  follows: 


NI 

NJ 

NPTHA 

NPTHB 

NPRNT 

NDUMP 

KOMPLX 

NITER 

YL.EN 

XLEN 

ZA-ZB 

IWWAL 

TOTITR 

TclKNLA 

TUKNLB 

NCOST 

NTYPE 

OTCORE 

OTSHKL 

FACTF 

FACTE 


No.  of  nodes  in  I-direction  (<■  30) 

No.  of  nodes  in  J-direction  (£  30) 

No.  of  paths  for  gas  side  (£  10) 

No.  of  paths  for  liquid  side  (£  10) 

option  for  printing  the  intermediate  iteriation  results  (“0  or  1) 

option  for  dumping  the  detail  calculation  for  each  iteration  («0  or  1) 

option  for  using  or  not  using  the  boiling  heat  transfer  model  (“0  or  1) 

maximum  No.  of  iterations  (default  value  *  23) 

overall  core  height  (■  £N  j A j  ^  inches) 

overall  core  width  (■  iN^AXj  inches! 

overall  core  depth  (inches) 

tube  wall  thickness  (inches) 

convergence  tolerance  for  iteration  of  the  exit  temperature 
of  the  fluids  (*F) 

factor  for  turn  loss  on  the  gas-side 
tactor  for  turn  loss  on  the  liquid  side 
option  for  cost  estimate  ("0  or  1) 
types  of  heat  exchanger  (1  to  3) 
types  of  tubes  material  ( 1  to  8) 
types  of  shell  material  (1  to  8) 
fabrication  complexity  factor 
escalation  factor  from  Hid  '70  dollar  value 
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In  addition  to  the  job  control  parameters,  there  are  several  sets  of  input 
data  which  were  used  to  describe  the  inlet  flow  conditions  and  the  nodal 
connection  method  for  each  flow  path  (see  line  6  to  line  28,  or  line  29  to 
51,  etc.  of  Table  A.l).  The  number  of  these  data  sets  is  equal  to  the  number  of 
gas  paths  (NPTHA)  plus  the  liquid  paths  (NPTHB).  The  first  two  data  cards  for 
each  data  set  contain  ten  parameters  which  are  defined  as  follows: 

WDOT  *  gas  flow  rate  (lb/sec) 

PZRO  *  gas  inlet  pressure  (psia) 

TZRO  ■  gas  inlet  temperature  (*F) 

DHYD  *  hydraulic  diameter  (inches) 

DELTAX  ■  nodal  width  (inches) 

FAOFA  ■  flow  area/frontal  area 
SAOV  *  surface  area/volume 
FINTHK  *  fin  thickness  (inches) 

FINLEN  ■  fin  length  (inches) 

FINSRF  *  fin  area/surface  area 

The  remainder  of  the  input  data  are  for  flow  direction,  number  of  nodes,  and 
nodal  connection  sequence.  The  liquid-side  flow  conditions  and  nodal  connections 
(lines  122  to  223)  are  similar  to  those  for  the  gas  flow  except  that  the  last 
five  parameters  are  replaced  by  a  NTUBES  parameter  which  is  used  to  specify  the 
number  of  tubes  for  that  path. 

The  development  of  this  heat  exchanger  computer  code  was  developed  to 
be  independent  of  the  working  fluids,  and  therefore  the  user  has  complete 
freedom  of  choosing  a  working  fluid  to  meet  a  specific  need.  Consequently,  the 
thermal  and  physical  properties  of  working  fluids  must  be  specified  as  part  of 
the  input  data  for  the  program.  For  application  in  the  present  study,  the 
properties  of  air  and  water  are  tabulated  in  a  special  format  as  shown  in  Table 
A. 2.  It  should  be  noted  that  each  data  set  is  preceeded  by  an  integer  number 
which  specifies  the  number  of  entries  to  be  read. 

The  thermal  and  physical  properties  for  the  liquid-side  working  fluid  are 
given  in  lines  No.  1  through  209  of  Table  A. 2,  Lines  No.  1  to  No.  24  (which 
were  not  used  in  the  present  study)  are  the  coke  (scale)  properties,  the  coke 
thickness  as  function  of  temperature,  and  the  coke  formation  hiatory  as  function 
of  time.  Lines  26  to  30  tabulate  the  aaturation  pressures  (psia)  and  temperatures 
(*F).  The  formats  for  these  entries  are  E1Q.5.  Lines  31  to  33  are  the  empirical 
constants  for  computing  the  convective  heat  transfer  coefficients  (see  Eqs,  la  to 
lc  in  Section  II  of  this  report)  for  laainsr-,  turbulent-,  or  supercritical-flow. 
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Lines  35  to  56  provide  the  heat  transfer  correlations  in  terms  of  Reynolds 
2/3 

number  and  StPr  for  the  vapor  phase.  Lines  58  to  66  tabulate  the  liquid 
properties  including  its  the  temperature  (*F),  viscosity  ( lbm/ f t-sec ) ,  thermal 
conductivity  (Blu/ft-F-Hr) ,  specific  heat  at  constant  pressure  (Btu/lbmF)  and 
density  (lbm/ft^).  The  input  format  is  also  in  E10.5.  Line  67  consists  of 
three  parameters  which  represent  the  critical  pressure  (psia),  critical  temperature 
(*F).  and  molecular  weight. 

Lines  68  to  119  are  the  tabulations  of  vapor  properties  which  consist  of  the 
pressure  (psia),  temperature  (*F),  density  (lb/ft^),  viscosity 
( lbm/ f t-sec ) ,  thermal  conductivity  (Btu/ft-F-Hr)  and  specific  heat  at 
constant  pressure  (Btu/lbm~F).  Each  parameter  will  be  read  in  E10.5  format. 

The  number  of  the  pressure  entries  is  specified  by  the  first  parameter  on  line  68 
and  the  numbers  of  entries  for  other  parameters  at  a  given  pressure  are  specified 
on  line  69.  The  second  parameter  shown  on  line  68  represent  the  type  of  the 
working  fluid:  1  for  distillate  and  2  for  pure  substance. 

Lines  121  through  125  tabulate  the  saturation  pressure  (psia)  and  the  heat 
capacity  (Btu/lbm-F)  for  the  boiling  mixture  if  the  boiling  heat  transfer  model 
is  not  used  (i.e.  KOhtLX»0).  Lines  128  to  139  are  the  pressure  (psia),  the 
temperature  (*F)  and  the  density  (lbm/ft^)  above  the  critical  point.  These 
data  are  also  read  in  E10.5  format.  There  are  three  sets  of  pressure  data  (as 
shown  on  line  126)  and  at  each  pressure  value,  there  are  four  sets  of  data  for 
the  temperature  end  density.  Lines  141  to  1  45  give  five  saturation  pressures 
(psia)  and  their  corresponding  values  of  heat  of  vaporization  (Btu/lb).  They 
all  have  the  same  input  format  of  E10.5. 

Lines  147  and  148  presents  the  data  for  temperature  (*F)  and  surface  tension 
(dyne/cm).  The  F-function  and  the  S-function  required  for  boiling  phase  heat 
transfer  computations  (see  Eqs.  3a  and  3b  in  Section  II  of  this  report)  are  presented 
in  lines  150  to  170  and  lines  172  to  186,  respectively.  Finally,  the  friction 
coefficient  (which  was  defined  as  Ap  **  4fpV^L/(2gO)  as  function  of  Reynolds  number 
are  presented  in  lines  188  to  line  209. 

The  gas-side  thermal  and  physical  properties  are  tabulated  in  lines  210  to 
265.  Line  210  specifies  the  molecular  weight  of  the  gas  while  the  temperature 
(*F),  molecular  viscosity  (lb  /ft-sec),  thermal  conductivity  (Btu/ft-F-Hr), 
and  specific  heat  at  constant  pressure  (Btu/lb^-F)  are  presented  in  lines  212 
to  230.  The  friction  coefficient  and  the  Stanton  numbers  as  functions  of 
Reynolds  number  are  given  in  lines  232  to  246  and  lines  248  to  262,  respectively. 
Finally  the  wall  (tube)  thermal  conductivity  (Btu/ft-F-Hr)  as  function  of 
temperature  (*F)  arc  tabulated  in  lines  264  and  265.  All  these  entries  are 
also  read  in  E10.5  format. 
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Sample  Results 

Based  on  the  input  data  shown  in  Table  A,1  and  A. 2,  the  computed  results 
for  heat  exchanger  performance  are  presented  in  Tables  A, 3  through  A. 7. 

Table  A, 3  presents  the  convergence  informal  ion;  i.e.,  number  of  iterations 
and  the  tolerated  errors  in  exit  temperature  of  the  working  fluids  between  the  last 
two  iterations.  The  summary  results,  which  include  the  flow  conditions,  the 
heat  exchanger  size,  and  the  manufacturing  cost  estimate  (not  shown  in  this 
example),  are  shown  in  Table  A. 4. 

Table  A. 5  shows  the  temperature  distribution  for  the  entire  heat  exchanger 
core,  including  the  inlet,  and  the  outlet  as  well  as  the  mean  temperatures  for  both 
working  fluids  and  the  wall  temperatures  on  each  side  of  the  tubes  for  each 
node.  The  distributions  of  the  convective  heat  transfer  coefficient,  the 
Reyno>ds  number,  and  the  overall  heat  transfer  coefficient  are  shown  in  Table 
A. 6.  The  liquid-side  pressure  loss  characteristics,  the  steam  quality,  and  the 
boiling  heat  transfer  coefficients  are  shown  in  Table  A. 7. 


List  of  Computer  Programs 

The  listing  for  Fortran  statement  for  the  four  major  computer  programs 
(HXMA1N,  HXINPUT,  HXCALC,  HXOUTPUT)  and  the  common  block  allocations  are  shown 
in  Tables  A. 8  to  A. 12. 
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FIG  A  1 


STRUCTURE  OF  WASTE-HEAT  BOILER  COMPILE*  PROGRAM 
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TABLE  A. 1  INPUT  DATA  FOR  MARINE  WASTE-HEAT  BOILER  STUDY 
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?  TWUNCON  Nir?r,f;J  =  5,NPTHA:S,NPTHt  =  l ,  NPRNT  =C  ,  NDUMP=0  ,6  0HPL*=  1 , 

3  YLEN270.':,  XLEN  =  120.0,2  A  =  84.  L,  ZB  =84,0,  THKW AL=" .^95 , 

«  T  CL  I  TRrb.C,TUfiNLA:C.r:,TuRNLB:C.u,NC0ST  =  ufNTYPL  =  4f 
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TABLE  A. 1  Cont'd 
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TABLE  A. 2  THERMAL  AND  PHYSICAL  PROPERTIES  FOR  WATER  AND  AIR 
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TABLE  A. 8  LISTS  OF  HEAT  EXCHANGER  MAIN  PROGRAM 


INCLUDE.  PPOCO,  LIST 
DIMENSION  T  B  A  o  (  ?,MAXY,MAXX  ) 
-SET  LATA 

DATA  LCTTA/2P244/ 

CATA  N C>\ S T R , NCK S AV / 1H , r  / 
DATA  NSTOWE/9/ 

I  SDYNr 1 

CALL  !\cLT(5irO,  i2DC > 
CONTINUE 

NU*RtC  r  2»  M  *  NJ 
MJ  M  w  W  0  =  10 
NSAVE  :  8 

CALL  ER TkAN ( o  ,  • 3A S& ,  T  S  .  • 
CALL  E  R  TR  A  N  (  6  ,  '  3  A  SQ  ,  T  S  .  • 


DEFINE  E a l E  NSAVl  (NUMRECiNUMWRD,U,LREC  ) 

CALL  CALC  1  L  0  T  T  A  )  ' 

1  UOu  continue 

DO  lir?C  1  =  1 .  U 
DO  line  jii.r.j 

11CL  CONTINUE^’'  r,5#(T^LU1'I'JI+Tu^IM,JII 
IF(  ISPYN.EJ.UGO  TO  1C50 

URITE(kw,6S0)  gnet,suhcp,dtau,tau 

b5l\ /?PvAIr  ’FP0*  RESPONSE  MODEL...'/ 

X/l  CX,'0-NET  =  ’  ,  E  1 2  •  5  *  1  x  ,  'BTu/SEC*  • 

X/10X,'SUUM*CP*DT)  =',E12.5,lx!'ETu', 

w}rJ’!?ALTA  T4U  -**|I2.5, IX, 'SECONDS', 
X/1CX, 'T  Aj  =  '  ,E  12.  S,  IX,  'SEC ONES' , 

1JSU  CONTINUE 

ISDYNZI ?Lyn*1 
CALL  lNpLT($3’“,J«l4ur  J 
ODD  CONTINUE 

CALL  CALC(LOTTA) 

do  I2nr  1 ; i , m 

DC  1200  j;l,Nj 


L20C 


CONTIN'U  J”J,:0'*5*<T‘,'4LL<  1,I,J)*TUALLC2,I,J)  1 

I - MASS  PER  NODE 

DMASS=ELMASS/(NI*NJ) 

SUMCP  =  C  .C 
DO  1300  1=1, SI 
DO  1  30C>  J=  1 , NJ 

TEE  =  ^.5*(ToAR(l,Itj)*TijAR(2,I,J)) 

CALL  LOCK  !Tcpvet.cpmet,ncpmet,tee,cpee,kk) 

13QD  CONTINUE^*  rM  A^S>>cPe£*  1TB  AR  (2,1  ,J)-T&AR(  1,1,  jl  ) 

OTAU  =  AB  S i 2.C*SUMCP/QNET ) 

TAU  =  TAU*OTAU 
GO  TO  l^T 

400  CONTINUE 

CALL  ERTriAM  b.'SFREE  9  .  •) 

WRITE (K  W, 600  INSAVE 

GOO  FORMAT <// 1UX , 'UNIT  ',15, IX, 'HAS  NOT  BEEN  FREED'//) 
200  CONTINUE 


NCK  S  A  V » NE • 2 ) GO  TO  210 
2 1 r 


IFINCKSAv.NE.l  .AND. 

WRITE  <K  k, 600  )NCKSTR 
IF<NCKSAv.NE.l)GQ  TO 
NDUm:ma XX*MAXY 

call  NT  ran (NCKSTR  ,22, 10,22,1 ,12,TITLE.LL,22. 

X  1,NDUMITMKCK,LL,22,10,22) 

WRITE«KW,610)NCKSTR 

610  FORMAT! /10X. 'COKE  THICKNESS  DATA  HAVE  BEEN  STORED' 
X  1 X  » *  ON  UNIT’, 15/) 

210  CONTINUE 

STOP 

END 
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TABLE  A. 9  LISTS  OF  HEAT  EXCHANGER  INPUT  PROGRAM 


SUBROUTINE  INPUT!*,*) 

If.CLuOE  PROCO,  LIST 
DIMENSION  N TEE (  10)  *  TYTLEI12) 

DtTA  KR , k  w/ 5 , 6  / 

DATA  P  I,GC/3.  1  4  1  592  ,  32.17**/ 

DATA  ELHASS/C.?/  .  „  „  ,  , 

DATA  NCP*ET ,NTYPt,MTCQRE,MTSHEL,FACTF,FACTF/n,  1  , 1 ,  1  ,  1. 0,7.3/ 


u«  I  B  V  ^  ^  U  I  I  •  I  »  r  t  «  V*  V  ^  L  f  *  1  *  w  N-  f  f  *  L 

NAMELIST  /INPUT/  NCPMET f CPMET ,TCPMET .ELMASS 
NAMELIST  /RUNCON/  NI  ,NJ  ,NPThA ,NPTHB ,NPRNT , NC 

.mao*  >  vtiKP  i.  *  a  *  n  u  t  i  ,  r  .  i-  »  (•  t  r  r  t 


p.homplx, niter, 

L  » 

TURNLA .TURNLfc 


READ  INPuTS  -  CHECK  FOR  ERRORS 


•LUMMY  PARAMETERS  SET  FOR  USE  IN  WRITE  STATEMENTS 
M  A  X  J  T  M  A  X  X 
MAXI  Z  K  A  X  Y 
M AXP  Z  h A  xP  TH 
MAXT  -  M  A  X  T  A  " 

maxn  z  maxnuD 


NERRzn 

IF!  ISDYN.Eu. 1 )G0  TO  1C010 
READIKR,  INPUT  ,ERR  =  lCOOO,ENO=incno ) 

ISDYnX=ISDYN-T 

WRITE I  K  L , 10035) I SCYNX  . 

10305  FORMAT  t  1h  1/ 1  Dx  f  «  >x>»*'>>x*m*m*ST  ART  ING  DYNAMIC  RESPONSES 
X  lx,  ’CODE  -  STEP  S  15.  IX  ,  ‘  ***.***»>.>**  •  /  ) 

go  to  irozu 

1CG1G  CONTINUE 
C - "“•CARD  1 

HE  AD ( K  R  ,  5HC )  (TITLE (I), 1  =  1, 12) 

SOU  FORMAT (  12A6 ,18) 

5  1  U  FORMAT (9110) 

52U  FORMAT ( SLlU.5) 

REALMS,  RUNCON) 

1GJ2C  CONTINUE 

C - CARDS  1C  -  19 

DO  1  DO  N=1,NPTHA 

READIKR  ,  5  1 L' )  1ST  .NODE  ,LSTEP 

1STARTI  1,N)=  1ST 

NODES ( 1  ,N )  =  NODE 

K  ANSTPI  N)  =(, 

IF (LSTEP.GT,C)KANSTP(N)=i  „  „ ,  , 

RE  AO ( KR  ,520 )  WDOT (1,N) , PZRO ( 1 , N ) ,TZR0I1»N) ,OHYQ( 1,N) , 

1  DE  L  T  A  X  i  to.) 

READIKR, 520)  FAOFAIN)  ,SAOVIN)  ,FINTHK  IN)  ,  F  I N  L  E  N  I  N  )  ,F  INSRF  IN) 
DO  11?  L=  1 , NODE 

READIKR, 511)  I.J,  TIN!  1 , 1 , J ) ,  HASIDEII.J)  ■ 
ill  FORMAT! 2I1C, 2E1D.5) 

512  FORMAT! 3E10»5,I1G,2E10.5) 

1ARAY I  1 ,N,L ) =  I 
JARAYI 1  ,N,L)=  J 
110  CONTINUE 
100  CONTINUE 
C . CAROS  21-29 


DO  2U0  N=1.NPTHB 
READIKR, 5 1G)  1ST,  NODE 
ISTARTI2.N):  1ST 
NODE  S I E  ,N  )  =  NODE 

READIKR, 512)  WDOT  < 2 ,N ) ,PZ*0 ( 2 , N ) ,TZRO 1 2 , N ) , 
X  NTUBESIN)  ,DH¥D(2,N)  ,OELTAYIM 

DO  210  L=  1  , NODE 

READIKR, 510)  I AP A Y I  2 , N , L ) ,  JAPAY(2,N,L) 

210  CONTINUE 
20U  CONTINUE 

IFIISDYN.NE.HGO  TO  10030 
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TABLE  A. 9  Cont’d 


C - CARD  72  COKE.  DATA 

READIER  ,  b  2C  )  XK  coke 
READ(KR,51C>  NCOKE 
DC  3ir  L=1,NCCKE 

RE  AO < KP  ,52b >  TCTAB(L),  THKCTIL) 

31u  CONTI  S’ Uf 

c 

RfAD(KP,52U)  H  ft  S  v  4  X 
KEAD(Kft,3li»  NHOuRS ,NCK SA V 
DC  311  L=1,NhCURS 
RMD(KP(d2H!  H  C  U  ft  S  (  L  t 
311  CONTINJf 
C 

C . CAftL  32 

READ(Kft  ,5!.  J  NS  A  T 
DC  3c?  L:l,NSAT 

READ ( Kt  ,52C  )  PS  ATT P (L >  ♦  T S ATTB ( L  ) 

3  2U  CONTINUE 
C 

C - CARDS  37-39 

C 

llSCOCr  1 

fiEAD(K»,s2C,EPft:721l  ALAM.BLAM, CLAM, CLAM 
READIER  ,  5  2C  >  A  T  U  ft  b  ,BTUP6,CTUftb 
READIKR, S2j>  ASUP,PSUP,CSUP,DSUP 
GO  TO  3Z2 

321  CONTINUE 

C - LIQUID  ^  FACTOR  USED  INSTEAD  OF  CORRELATIONS 

L  IOCOP: 2 

RCAL' (  KR  ,B1G  )  NLJAY 

DO  32T  LA  1, NLJAY 

READIKRorOftENLICCD.STNLIOa) 

3 23  CONTINUE 
r 

322  CONTINUE 

READ ! kR  ,5  ID  )  NSTNTB 

DO  3  3n  LA1.NSTNT3 

RE  A  0  (  K  R  ,  5  ?  D  )  REN STB  <L ) , STNTB (L  ) 

33 J  CONTINUE 
C 
C 

C - CARDS  4  E-49 

C 

READ  <  KR  ,510)  NT A3B 
DO  42C  L=1,NTABB 

READ(KR,52C)  TEM3  T  (  L  )  t  VI SBT ( L J ,  XK8T < L 1 , CPB T ( L ) » RHGBT ( L ) 
420  CONTINUE 

READ(KR,520)  PCRITB.TCRITB,  A“UB 

c - FUEL  VAPOR  PROPERTIES 

READiKR  ,510)  NP.ISPURE 
NPR0PS=5 

NTRY  ;  NPWOPS*  1 
READ(KR,5lO)(NTEE<L)fL:l,NP) 

T  AB  V  A  P  <  1 ) =NP 
T  A  B  V  A  P (  2) ANPRCP  S 
DO  43^  L=1 ,NP 
TABVAP(?*L)=  NTEE(L) 

430  CONTINUE 

LAST:  2 *2*NP 
00  43?  N=1,NP 
NTrNTEE  <N ) 

DO  4  34  K: 1 j  NT 

READ ( KR  ,520  >  PEE,  ( T A B V AP ( L A  ST *L ) , L  =  1  ,NTR V ) 

I E  <  K , E 0  •  1 )  TABVAP ( 2  *NP ♦ N >  APE  E 
LASTaLA ST^NTRV 
434  CONTINUE 
C 

432  CONTINUE 

C - MIXTURE  HEAT  CAPACITY 

READIKR  *  5  l  0 ) NM I  X 
DO  438  LA1.NMIX 

REAOIKR  ,520 )  PM  I  X TB ( L 1 , CPM I XB < L > 
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436  CONTINUE 


-----DENSITY  ABOVE  CRITICAL  CONDITIONS 
REACIKR  ,61C)NP 

NPROPS;iblC1  ,NTmL,*L:l»*p> 

NTRY  r  NPHOPS+1 
TABCRTI 1 > :  NP 
TABCRT l 2) R  NPPOPS 
DO  440  1=1,., P 
T  ABORT ( 2  +  L  )  R  N  T  £  E I L i 
44^  CONTINUE 

LASTR2*2*NP 

DC  440  NR  1 , NF 
NTrNTEE IN ) 

DC  444  k:  1  ,NT 

SfISIeS  :!?°T4ECst ,L-1 ,KTBV , 
L  A  S T  RL  A  ST  +NTS  Y  - 
444  CONTINUE 

4  4  CONTINUE 


- HEAT  OE  VAPORIZATION 

READ! KR  ,51U )  N  V  A  P  T  P 
DO  46C  L=l,NVAPTb 

4 8 U  CONTlSS#520*  rL4MTB“-J*HVAPTB(L» 


•SURFACE  TENSION 


READ  I KR  ,5 1L I  S  I  u  *  A 
00  4QI  L~1,NSIGUA 

481  C^rl^fS2C>  ISIC-MAILS.fSIfiTABlL) 


- F  FUNCTION  F C c  BOILING 

REAOIKR  ,51C)  NOV  XF 
00  4  B  2  LR  1  i  NOVXF 

com{55|“c’  L..ro,. MB, L) 


- S  FUNCTION  FOR  BOILING 

READIER  ,sil>  nstab 

DO  4  6  3  L  =  1  ,NSTap 

MB  CONUNuf521'1  S*Ett3ILI.ST«6,ll 

•— 5?iEII^5i8rif8iEM 

DO  446  L=l,NrR3 

»M  Coft9ifiS£S2!-'  e“F8*Ll.FBMB,L, 


• - CAROS  5  T- 5 9 

REAOIKR  *520  AHUA 
gEA0(KR,510)  NT  A  B  A 
00  450  L=1»NTABA 

45u  CONTINu£S?U>  T^M*T(L  >  i  V  1 5  A  T (L ) » XK  A T (L ! ,  CP  AT (L ) 

REAOIKR,  51C)  NF  P I C 

DO  460  LRl.NFPII 

REAOIKR, 52D)  RE NF  ( L  > , F T AB I L 1 
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•♦60  COST  I  NU  r 

REAC(K»,S1.)  \  S  T  A  N  T 

DC  472  L=1,NSTANT 

R  E  A  0  (  X  R  ,  b  2  L  >  PESSTIL  )  .STSTABIL  ) 

4  7  L)  CONTINUE 

- CARL'S  6  0-69 

READ! KR  , SID  )  NWALK 
DC  6  2 C  L 1 1 , N « A L  K 

»rAD(KR,S?r»  T  »  T  A  h  I  L  I  »  XKWTAB(L) 
t >20  CONTINUE 

END  or  INPuT 


WRITE  INPUTS 


IFINPPNT  .EQ.  D  >  CO  TO  74  71 
WRI  TL  <K  W, 6. rw ) 

6230  F  C P M A  T «  1M  ) 

*RITE(Kw,6Clw)  « T  I  T  L  E  I  I  )  ,  1  =  1  , 12 I 
6  u  1 C  FCR4ATI  /10a, 'INPuT  SPEC  IF TCAT  IONS  OF  ' , 1 2  Ab / > 

wpite(kw,6l2li  m,nj,nptha,npthr,nprnt,noump,kcmplx 

6  i 2.  FORMAT ( 

X/ILX.'NC.  OF  NODES  IN  Y-CIR.  ',I1L, 

X/1CX,'N0.  OF  NODES  IN  X-DIR.  '.lie, 

X/IDX.’NO.  OF  A  -  PATHS  ',110, 

X  /  1  :j  X  ,  'NO.  OF  P  -  PATHS  *,I1L, 

.110, 

tile, 

.110. 


X/ILX.'NC.  OF  NODES  IN  Y-CIP.  ',I1L, 

X  /  1 C  X  ,  ’NO.  OF  NODES  IN  X-DIR.  '.lie, 

X/IDX.’NC.  OF  A  -  PATHS  ',110, 

X/l'jX.'NC.  OF  P  -  PATHS  ’.IlL, 

x / 1  c x ,  vjPkNT  ’.iio, 

X/1DX,  *  N  p  u  **  P  MU, 

X/1CX, 'K CMPUX  ’,110, 

X/  ) 

.RITEUW.6L3L)  YLEN,XEEN,2Af 2PSTHK*AE,T0LITR, 
X  T LRNL A  ,  TuPNCb 

6,30  F0RMA7( 

X/12X , 'Y -SIDE  LENGTH  » , E U . 5 , 1 X , • INC 

X/1DX , ’ x -SIDE  LENGTH  *,£12.5, lx, 'INC 

X/1GX , *A -SIDE  DEPTH  ’  ,£12*5, lx, 'INC 

X/1GX, '3-SIOE  CEPTH  *,£12.5, IX, 'INC 

x/Icx,  *•  all  thickness  •  »Ei2.5,Ix,  mnc 

x/!ux, • I TtRATION  TOLERANCE  •  ,  E 1 2 . 5 , 1 X , • OE G 

X/  1 0  X ,  '  S  IDE  A  TURN  LOSS  FACTOR  *,E12.5, 
X/1CX,*SIDE  B  TURN  LOSS  FACTOR  *,£12.5, 

X«RITt <K w, 6031 )  S.'EEP 


*  ,£12*5, IX, 'DEGRtES*/ I 


•  ,  1 12  , 

' ,112, 


6031  FORMAT ( 

x/iox,  'Sweep  angle  ',ei2.5,ix,'degri.es’/) 

WRITE (K  k, 6049 ) 

6049  FORMATJ  lnl ,/irX  ,  'A-SIDE  PATH  DESCRIPTIONS'/! 

C 

10030  CONTINUE 

CO  7uPD  N=1,NPTha 

WRITE  I K  W ,  6:.5w  IN, 1ST  ART  (  1 ,  N  I ,  NODES  <  1 ,  N  > ,  KANSTP(N)  ,  wDOT( 1 ,N> 
X  P  2R0  ( 1 ,  N I  ,  TZRO  ( 1 ,  N  I  ,  DHYDU.N) 

6050  FORMAT < 

X / 1  OX , 'PATH  NUMBER  ',112, 

X/ICX, 'START  INDICATOR  ',112, 

X/JOX, 'NLMPER  OF  NOOES  ',112, 

X/IOX, 'STEPPING  SWITCH  ',112, 

X/ 

X/lOx,  'FLOWRATE  »»E12*5,1X,'LBM/S£C'» 

X/1UX, 'INITIAL  PRCSSURE  •  ,E  1 2 . 5 , 1 X , • PS  I  A  ', 

X/IOX, 'INITIAL  TEMPERATURE  •  ,E  12 . 5 , IX , • OEG  F 

X/IOX, 'HYDRAULIC  DIAMETER  », £12. 5, IX, 'INCHES  *, 

X/  I 

WRITE  <  K  W,  606C  I  F  AOF  A  ( N  )  ,  SAOVIN),  FJNTHK  (  N  )  ,  FINLENIM, 

X  FINSWF(N) 

6060  FORMAT! 

X  1  OX , 'FLOW  ARfA/FRONTAL  AREA  »,E12.5, 

X/IOX, 'SURFACE  AREA/VOLUME  *,E12.5,  IX,'FT*6-1  ', 

X/13X, 'F  IN  THICKNESS  ',£12. 5,  IX, 'INCHES  ', 


»,E12.5,lX,'LBM/SEC*t 
',E12.5,1X,'PSIA  ', 

', £12. 5, IX, 'OEG  F 
», £12. 5, IX, 'INCHES  *, 


I X , *FT**» 1 
IX, 'INCHES 
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LENGTH 

ARE  A/SURFACE 


X/iCX  ,  '  F  IN 
X/1CX ,  *F  IN 
X/  ) 

NODE  =  NODES! 1 »  N  ) 
IF  (  KAN'S  TP  <  N  )  .C-T  .  . 
WRITE (K  W»  6L7C  ) 
b  C  7  L  FORMAT!  / 1  Q  X  •  'NOPE 
/  (  /  IPX  ,  Ig ( I  3  , 

GO  TO  7  COO 
7  u  i  J  CONTI NU  F 

WRITE (  K  W ,  6  715 
6371  FORMAT (  /1RX.  * 

X  /IP*,  8  X  , 

DO  7L2C  L  =  l,NODF. 

I  =  IAPAY(1,N,L ) 

J  =  JARAYI  1  *  N  ,  L  ) 
WRITE <K W, 6u7E  ) 
t-  372  F  ORMA  T  l  1-X , 21 4 , 

7  3  2-j  CONTINUE 

wRITt (KW.6D73 ) 

6  j?i  FORMAT (  ///  ) 

7jQU  CONTINUE 


AREA 


’  *  E  1  2  »  5  , 


i*.  *  ' 


> GO  TO  73  lr 

(IAKAY(l,N,L),JAPAY(lfN#L),L: 
CO-ORDINATES  <  I  ,  J  I  =  <  Y  ,  X  )  -  5  , 
•,•,13,3X1)  ) 


1 ,NO-E  ) 


I 


TEMPERATURE 
DEG  F 


H  * 

btu/hr-soft-f • ) ) 


IARAYU,N,L),JARAY(J,\',L),TIN<1,1,J),HASIP£(I,J) 

F 1 2 • 3 ,  E  1  4 , 5  ) 


WRITE  IK  k, 6130 ) 

6 13U  FORMAT (  lh 1 , / IPX 


- .  / ii  *  ,  ■ 

DO  7 i 00  N=l,\PTwj 
WRITE <K W.6L53  )  N.ISTARTI2 

X  P2RO<2,N)  ,  - - 

6353  FORMAT! 
x/ 1GX  , 

X/ f uX  , 

x  /  ’  :  x , 

X/10X, 

X  / 

'  X/ICX, 

x/iox. 


X  /  1 

N0DE  =  N0rLS(2,M 
WRITE  (K  w,  6-7-  1  (IARAY(2,N,L) 
71P0  COM  I NU  r 

IF ( ISDYN.NE . 1 >G0  TO  1PC40 


•B-SIDE  PATH  DESCRIPTIONS*/) 


ISTART  (  2  •  N  )  ,  NODES  (  2  ,  N  1  ,NTUBES(N),W’CCT(, 
T  Z  PO ( 2 , N ) ,  6hyD(2,N) 


,  N  » 


•PATH  NU“BE  R 
•STAPT  INDICATOR 

•NLMGLR  Of  NODES 

•  ,112 

•NLM9ER  OF  TLBES 

*  *  1 1 2 

•FLOWRATE 

•  *E12 

•INITIAL  pressure 

•  ,E12 

•INITIAL  1EMPEPATUPE 

•  ,E12 

•hY-RAUlIC  DIAMETER 

•  ,E12 

•TLGE  DENSITY 

•  ,E12 

: ,  5 , 1  X  , 
.  S  *  1  X  , 
>  5 , 1  X , 
.5, IX, 


•PSIA 
•DEG  F 
• INCHuS 

•IN**-2 


JARAY (2,N,L ) , L= 1 , NODE  ) 


WRITE  UW,611C»  XKCOKE 

6110  FORMAT  (// IL'X  ,  ’COKE  THERMAL  COND.  *,E12.5,1X, 

X  *Ei  Tu/LBM-F  T-HR*  /  ) 

WPITE  CK  W, 612C  ) 

612-  FORMAT! /IPX, *  N  TEMPERATURE  COKE  THICK*, 

X/15X, ’  DEG  F  INCHES’/) 

00  7110  N=1,NC0KE 

WRITE !K W,613C  )  N,  TCTAb (N  > ,  THKCT(N) 

613-  FORMAT!  1-X, 13, 2X,  2EI2.5  ) 

7 1 1U  CONTINUF 

C 

c 

WRITE(KW,6131)HRSMAX 

b 1 3 1  FORMAT! /1DX.*C0KE  CURVE  REPRESENTS  CONDITIONS  A T • , E 1 2 . 5 ,  1 X , 
X  ’HOURS’/) 

IMNCKSAV  .EO.  1)  WRITE  CKW, 6121  INCKSTR 
IF ( NCKSAV  « E 0  •  ?) WRITE (KW ,6122 INCKSTR 

6121  FORMAT! /10X, ’SAVE  CCkE  THICKNESS  DATA  ON  UNlT’,15/1 

6122  FORMAT! /10X, ’READ  COKE  THICKNESS  DATA  FROM  UNIT’,15/1 

WRITE (KW,6132) 

6132  FORMAT!  //1QX, •  N  T I ME  * / 1  OX , JX , 2X , 7X , *  HOUR S * / ) 

DO  7111  N=I,NHOURS 
WRITE  !K k, 61  30  )  N,HOURS(N) 

7111  CONTINUE 
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PSI A  TEG 

P  S  A  T  T  9  !  N  )  ,  T  S  A  T  T  B  (  N  ) 


/  ) 


W  R  I  T  L  !  K  W  ,  6  1  4  0  ) 

6 1 4  j  FORMAT (  /10x , ‘SATURATION  TEMPERATURE  TABLE’, 

X  /IDX,’  N  SAT, -PRESS  SAT. -TEMP’, 

X  /I  OX  ,  ’ 

DC  71 ?P  N-l.NSAT 
•  RITE (K  w,613D)  N, 

712u  CONTINUE 

GC  TO  (7121,71221,  LlQCOk 

7121  CCNTINUF. 

WRITE (K *,62001  AL AM ,BL AM ,CL AH  ,DL AM  ,  AT URfa , BT UR b , CT UR B  , 

X  A  SuP.BSl'P.CSUP.DSUP 

6232  FORMAT (  //IDX, 'HEAT  TRANSFER  CORRELATION  CO -E F F I C I F NT S ’ , / 

X/IGX.’FLOW  A  i‘  C  O’ 

X/1DX,  '  L  A  K I N  4  n  ’  ,  4E 1,  .5, 

X/1CX, ’ TIRBOLE’ T  ’,  3E1C.5, 

X/ 1CX  , ’SLPERCR'IT  .  ’,  4E10.S, 

X/  ) 

GO  Tu  7129 

7122  COM  I  ML*  F 
»PITE  (*  k,  620  j  ) 

6225  FORMAT (  ZiOx , ’5-S  ICE  STANTON  NUMBER  TABLE  (LIQUID)’/ 

X  / 1  c  X  ,  ’  \  ’,  •  REYN  NO.  $TvPR**2/3V  1 


DO  7123  Ml.NLJAY 
WRITE CKW, 631 3)  N. 
7123  CONTINUE 

’/12V  CONTINUE 


REYN  NO 
RENLIO(N),  STNLIQ(N) 


STANTON  NUMBER  TAILE  (VAPOR)’/ 
REYN  NO.  ST*PR**2/3 ’/ ) 

?R  r  N  ST  3  (  N  )  ,  S  T  N  T  B  (  N  ) 


WRITE  «  « W, 62  1L ) 

621.  FORMAT (  /i^x, ’r-SIDC 

X  /  1 C  X  ,  ’  \  • ,  • 

DO  713C  NRl.NSTNTB 
»  R  I  T  L  (  K  «  ,  k  3  1  l  )  N 
7  1  3  w  CC.NTINUF 

WRITE  Uk,  6  3(10) 

6303  FORMAT  (  lHl, /MX, ’LIQUID  THERMAL 
X//1CX ,  ’  N  TEMPERATURE 

X  1  SPEC  HEAT  DENSITY’  . 

X/15X,  7  X, *  DE  0  F*,'  LBM/F  T-5EC  BTU/FT-F-HR  BTU/LGM-F • , 

X  ’  LbM/CU.FT*/) 

CC  7  3C0  N;1,NTAFB 

WRITE (wk, 6313) N,  TE“PT(N),VIS9T(N),XKBT(N),CPBT(N),PH0BT(N) 
6313  FORMAT!  1 0  x  ,  I  3 , 2  X  ,  S  £  1 2 . 5  > 

73Du  CONTINUE 


PROPEPTY  DATA  ’, 
VISCOSITY  THERM 


COND 


WRITL(Kk,632u) 
632.  FORMAT! 

X/1GX,  ’LIQUID 
X/  1CX ,  * 

X/1UX,  • 

X/  ) 


PCRITb.TCRITB, AMUB 

PRESSURE  !  , E 1 2 • 


CRI  T  . 

CRIT.  _  _ 

MOLECULAR  WEIGHT 


TEMPERATURE  ’  ,Ei2tS,Ix - • 


, E 12 . 5 , 


5, IX, 


PSIA’  . 
DEG  F* 


wRITE(KkM33  :  JISPUEE 

6  3  3C  FORMAT!  1H1./1CX, ’e-SIOE  VAPOR  THERMAL  PROPERTIES  USPLRE  =  ', 

X/1DX,’  NP  '  NT  PRESSURE  TEMPERATURE 

-  THERM  COND  SPEC.  HEAT 


X/i:x,12X,  bX  , ’PS  I A •  ,  7X , *  DEG  F’, 
X  ’  LBM/FT-SEC  BTU/FT-F-HR 

A.  Am  V  A  A  II  A  A  I  1  I 


6  34  C 

7  32  . 
7  313 


NP:  TADVAPI1) 

NTRY:  TAtjVAP(2)*l 
LAST:  2»2*NP 
DO  7 3 1C  N-  1 , NP 
NT:  TABVAP(2*N) 

PEE:TA3  VAP! 2*NP*N ) 

DO  7  320  K :  1  ,  N T 
wPITE (KW.6340)  N ,K ,PEE  • ( 
FORMAT (  13x,2I5,2x,7E12.5) 
LAST :LA ST*NTRY 
CONTINUE 
CONTINUE 


BTU/LBH-F 


DENSITY*  , 
QUALITY’ , 

LEM/CU.FT' , 
PERCENT’/) 


PEE . (TAB VAP( LASTED ,L=1 , NTRY) 
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WRITE(KW,6140) 

614.  FORMATl  /inx, ‘SATURATION  TEMPERATURE  TABLE’, 

X  /10X,’  N  SAT. -PRESS  SAT.-TEMP’,  . 

X  / 1  0  X  ,  *  PSIA  TEG  F  *  n 

DC  7  N:1 ,NSA  T 

«RITE(KW,613L)  N»  PSATTB(N)  ,  T  S  A  T  T  B  (  N  ) 

TICw  C  ONT  I  NU  r 

GC  TO  (  Tin, 7122)  ,  LIQCOR 
CONTINUE 

WRITE (k  .,620C)  A L AM , BL A M , CL  AM  , DL AM  ,  ATURb , BT URb , CTURB  , 

X  ASuP,BSl'P,CSUP,OSUP 

6202  FORMAT (// 1CX , *HE AT  TRANSFER  CORRELATION  CO - E F F I C  1  F N T S * , / 

X/  1GX  ,  'FLOW  A  ?  C  O', 

X  /  1  0  X  ,  *L  AKlNAR  ’,  4E1,.S, 

X/1CX, *TLRDULENT  ’,  3E1C.5, 

X./1CX,  ’SLPERCRIT  .  ’,  4E10.S, 

X/  ) 

GO  Tu  7 129 

7122  COMINl'F 
wPIU  (Kw,62as  » 

620b  FORMAT  (  /iOX  ,  ’B-SIDE  STANTON  NUN.  BE  P  TABLE  (LIQUID!  *  / 
x/icx,*  *,  ’,  '  REYN  NO.  STvPR**2/3’/  ) 

DO  7123  Nrl.NLJAY 

WRITE (kw,631u)  N,  RENLIQ(N),  STNLIQ(N) 

7123  CONTINUE 
C 

7129  CONTINUE 


WRITE f» W, 621L) 

621.  F  OR  M A  T (  /  irx, ’  P -S I  DC  STANTON  NUMBER  TA2LE  (VAPOR)’/ 

X/1CX,’  f,  •  REYN  NO.  ST  *PR**2/3  ’  / ) 

DO  713C  Nrl.NSTNTB 

WRIU  (*  «,fc3ic  )  N,  RFNGTS(N)  ,STNTB(N) 

713^  CCNTINuF 

WRITE (K  W,h3f)Q> 

6  30u  FOT?MAT(  lHl,/lrX, ’LIQUID  THERMAL  PROPERTY  DATA  •  , 

X//1CX.  ’  N  TEMPERATURE  VISCOSITY  THERM  CO.ND’, 

X  1  SPEC  HEAT  DENSITY’, 

X / 1 5  X ,  7  X, ’DEG  F*,’  LBM/FT-SEC  BTU/FT-F-HR  BTU/LUM-F  *  , 

X  •  LbM/CU.FT'/) 

DO  7  3C0  N=1,NT*BD 

WRITE (KW, 6310)  N,  TEMBT(N),VIS9T(N),XKBT(N),CP8T(N),PHOBT(N) 

6  3 1 u  FORMATl .  1JX, I3.2X , 5E12.5) 

7  3DG  CONTINUE 
C 

WRITE (KW.632D)  P CR I T B , T CR  I TB , A MUS 
632-  FORMAT! 

X/ 1DX , *L IQUID  CRIT.  PRESSURE  ’,E12.5,lX,’  PSIA’, 

X/1CX,’  CRM,  TEMPERATURE  ’  ,E12. 5,  IX  ,  ’  DEG  F‘, 

X/10X,'  MOLECULAR  WEIGHT  ’,E12.5, 

X/  ) 

C 

wRITE(KW,623  :)ISPURE 

6  3  3  C  FORMATl  1H1 , /1CX , • B-SIOE  VAPOR  THERMAL  PROPERTIES  (ISPLRE:*, 
X  12  IX  •  )  1  / 

x/iux,’  N P  ’  NT  PRESSURE  TEMPERATURE  DENSITY’, 

X  •  VISCOSITY  THERM  COND  SPEC.  HEAT  QUALITY', 


X/1DX.12X,  bX, ’PSIA’ ,7X, ’DEG  F' ,  *  LbM/CU.FT', 

X  ’  LRH/F  T-SEC  BTU/FT-F-HR  BTU/LSM-F  PERCENT’/) 

NP=  TIB  VAP (  1  ) 

N  TR  Y  r  TAdVAP(2)*l 
LAST:  2  *2*NP 
DO  7 i 1C  N=1,NP 
NT:  TABVAP(2*N) 

PEE=TA3 VAP(2*NP*N> 

DO  7320  H :  1 , N T 

WRITE (KW .634Q)  N,K,PEE,(TABVAPILAST*L),L:1|NTRY) 

6340  FORMAT!  l.jX,2I5,2x,7El2»b) 

LAST:LAST*NTRY 
7  32  1  CONTINUE 
7310  CONTINUE 
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«  R  !  T  c.  (  k  «  ,  6  3  5  u  ) 

o  3  5  0  FORMAT! //ICX.’MJxTuPE  H  E  A  T  CAPACITY  TABLE’/ 
x/iDX,*  *,  •  PRESSUPE  SPEC.  HEAT*, 

x/1Ux,  fx,  •  PSIA  BTU/LBM-F*/) 

DC  7  3  3  C  Kz 1 ,  r.  *  I  x 

*»I  TE  (  K  .,6310  N , PHI  XTB (N  >  ,  CPMI XE  <  N  > 

7  3  3  u  C  0  -  *  T  I  \  u  r 

.RITLIw.,6360 

C36-  f  n  R  M  A  T  !  //*ux,*pf\sity  above  CRITICAL  point*/ 

X/1CX,*  *.p  NT  PRESSURE  TEMPERATURE  DENSITY*, 

X/1DX,  l'x,  *  PSIA  CEG  F  LBM/Cw.FT*/) 

NP: T ABC  RT  t  1  ) 

N  T  R  Y  :  T  A  d  C  R  T  (  2  )  ♦  1 
lAST-0* 

u  0  7  3  4  \  z  i ,  \  P 

NT:  TAB  CRT  12**.) 

PEERT  ACjRRT  I  2  ♦  N  P  ♦  N  ) 

DO  7  3  4  5  •»  -  1  ♦  T 

«RITL(Kw,f34u)N,<\,PEE,(TA9CRT(LAST*L),L:i,NTRY) 

LASTiLA ST *NTRv 
7J4S  CONTINUE 
734D  CONTINUE 
C 

WRIU  (.  fe3SC  ) 

6  36  J  FCRmaTI  //luX,  ’HEAT  OF  VAPORIZATION  TABLE*,/ 

X/lOX,’  N  ’,  *  PRESSURE  H-VAP*, 

X/1CX,  Ex,  *  PSIA  PTU/L9M'/) 

DC  7  3  BO  N=i,\VAPT3 

.RITE (k «,631C )  N,  PLAMTB(N) .HVAPTBTN) 

7383  CONTINUr 
C 


6  38  1 

7  39  1 


.  R  I  T  L  (  k  k  ,  6  3  8  1  ) 

FORMAT!//!  X, 'SURFACE  TENSION 
X/1DX,*  \  *,  •  TEMPERATURE 

X  /IDX ,  Ex ,  *  DEG  F 

DO  7  3  B 1  N:i,NSIG“A 


table*/ 

SIGMA*  , 
DYNES/O’/  > 


.PITE!*V,6310  N,  TSIGm  At  N  >  ,  SIGTAB  (N  ) 

C  0  N  T  I  N  U  r 


wRITl(kw,6392> 

6  38  2  FORMATI  //1LX, ‘F-FLNCTICN  FOP  ROILING*/ 

X/1CX,*  f,  *,  *  1/9IGX  F*/) 

DO  7382  Nri.NCVXF 

WRITE(KL,631G)N,  XOVFTO«N»,FOVFTB(N> 

7382  CONTINUE 


WRITE  CK .,  628  3  ) 

6333  FORMAT!  //1UX,  'S-FUNCTION  FOR  BOILING’/ 

X  / 1 G  X ,  *  (,  ',  •  FE*F#*1,2E  S  *  / » 

DC  73P4  n=1.NSTA'3 

WPITC(KW,63I3)  N,  SRELTB(N) , STAB tN> 

7384  CONTINUE 
C 


WRITE  <K  W, 6370  ) 

6  370  FORMATI  // 1LX, *B-SIOr  FRICTION  FACTOR  TABLE*/ 

X/1CX,*  N  *,  ’  REYN  NO.  F-F  ACTOR  * / > 

OG  7350  N -  1 , NFR B 

WRITE (W W, 63iw >  N,  R  E  NF  B  1  N  I  »  F5TABIN) 

7  35C  CONTINUE 


WRITE ( K k,64Ct 1AMUA 
6  400  FORMAT!  1H1./1CX, ‘A-SIDE  THERMAL  PROPERTY  DATA  *  . 

X/1CX  ,  ’MOLECULAR  (.EIGHT  SE12.5/ 

X/lOX,*  N  *.  ’  TEMPERATURE  VISCOSITY  THERM  COND*. 

X  *  SPEC  HEAT*  , 

X/15X,*  DEG  F  LBM/FT-SEC  BTU/FT-F-HR  BTU/LBM-F • / ) 

GO  7400  f.rl.NTABA 

WRITE(KW,631C>  N,  TEM A T ( N  )  ,  VI  S A T ( N  ) , XK AT < N » , CP AT < N ) 

740C  CONTINUE 
C 
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* 

k' 


£,4  50  FCR-AT(  1hi!/1PX,**-SIDE  FRICTION  FACTOR  TABLE*/ » 
x/l-x,*  N  *,  *  REYN  NO.  F-FACTOR*/) 

D  C  *  74  5C  N-IjNFRIC 

bR  I  TE ( K W , 6 3 10  )  N,  RENFINl,  FTAB(N) 

7450  CONTINUE 

WRITE (Kb, 6460) 

64fe.  FORMAT!  //1GX,  'A-SIDE  STANTON  NUMBER  TABLE*, 

//  1GX,*  N  *  REYN  NO,  ST  *PR  **2/ 3  * /  ) 

DO  7460  N-l  ,NST  ANT 

bPI  TL ( K k, 6il U  1  N  ,  RENST(N),  STNTAB(N) 

7460  CONTINUF 

k  R  I  T  E  (  k  w ,  6  4  7  ;  ) 

647C  FORMAT! //lit, ’WALL  THERMAL  CONDUCTIVITY  TABLE/, 

x/iox,*  *,  '  temperature  therm  cond*, 

X / 1 5  X , *  DEG  F  BTU/FT-F-hR*/ J 

DO  7470  U  =  1 , N W A L H 

WRITE  (k  ■»,  6311  )N  ,  TwTABINI,  XKWTAb(N) 

7  4  7  j  CONTINUF 

7472  FORMAHIHI,*  P-SIDE  FLOW  RATES  FOR  EACH  PATH  ARE; 
1  E  1 2 . 5  >» 


IF(ISPuRE.LT.l)ISPURE=l 
IF (  ISPuRt  •  GT  .  2 )  ISPUBE=2 
IF ( I SPURE .GE . ?  .AND.  NPRNT.EO.D) 

I  IWD0T (2,N),N:1,NPTHB) 

CONVlRT  TO  A  CONSISTENT  SET  OF  UNITS 

•  I  POUNDS 

T  T ME  -  SECONDS 

TEMPERATURE  -  RANkINE 


bRITE(KW,7472) 


1D04C 


XLEN 
YLEN 

ZA 

ZB 

THKwAL 

XKCOKE  = 
PCRITB  : 
TCRITB  = 
CONTINUF 
DO  auDD  n: 
P  ZRO  t  1,M 
PZRO ( 2, N) 
TZRO  l l, M 
T  ZRO (  2,  M 

OHYD ( 1 , M 
DH YD  12. A  I 
OELTAX ( N ) 
DELTAYt  N) 

F INTHK I M 
F INLEN! N) 


;  XLEN  /K.% 

;  ylen  / 12.0 
;  za  /  1 2  •  0 
=  29  /12.0 

:  THKWAL/12.0 

XKCOKE  / 36CC.D 
PCRITB*144  .C 
TCRITB>4tO.O 

:  1  , M  AX  P 

:  P  Z  R  0  <  l ,  N  )  * 
P  2  R  0  I  *  ,  N  )  * 


KO I t ,N1  * 

T  J  R  C (1 ,N)*46w,- 
TZRO ( 2 .N J  *460,0 


144.0 

i:8*c 


DHYD ( 1 ,N  » 
DHYO (2,N ) 
DELT AX (N ) 
TELTAY IN) 
FINTHK ! N ) 
FINLEN  <N )/ 


/  1  2 .  D 
/ 1  2.C 
/l?.r 
/  1 2 . 0 
/ 1  2 . 0 
12.0 


9  gOQ  CONTINUF 

IF ! ISOYN.NE .  1  S G 0  TO  1P0SC 
DO  8100  NRl.MAXT 
TCTAB(N)  =  T C T AS  )  N ) ♦ 460 .0 
THKCT(N)  =  ThKCT(N»/12.0 


PSATTB ( N ) 
T  S  A  T  T  B  (  N ) 

TEMPT (N  ) 
XKB  T l N ) 


PSATTB(N)*144.0 

TSATTB(N)*46u.C 

TEMBT (N  **460.0 
XKBT (N! / 3600 . 0 


// 1C ( lx  , 
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C 

C 

C 

C 

C- 

C 


P  m  I  X  T  <7  (  N  )  =  P‘'IXTB(N)*l44.a 

tekat  «n  )  ♦**6C.r 
x  k  A  T <N>  /360C.C 


T  E  *  A  T  ( N  > 

X  K  A  T  (  \  ) 

T-TAd  in  ) 

X  K  T  A  0  (  K  > 

P  L  A  «  T  B  <  r. ) 


;  TwTAri(N)*4bC.n 

=  x^wTAei v'/36cn.r 

:  PlAHTP (N)*144.C 


TSIGMAIN)  -  TSIGMA(N)  ♦  460. C 
•DYNE/C*  70  Lt>F/FT  „  „  ,  Q 

SIGTAfc(N)  =  S 1 0  7  A  ;•  J  N  )  *  2.248E-6*  3C.48 


8100  CONTINUE 


8210 

820C 


N p ;  7  Ac  VA  P  (  1  > 

NTRY  =  7APVAPI2I*1 
LAS7  =  2*  2HP 
DO  2  2**t  N  =  1  *  N  p 
NT  =  7AbVAP(2*M 
7  A  b  V  A  P (  J*NP*N>: 

DC  6  210  K  =  1  »  M 
7  A B  V  A  p  (  LA  S  7  ♦  1  ) 

7ABVAPILAS7*  4) 

7 A 6VAP<  LA $!♦  £  ) 
LAj7^LA  E  T  ♦  N  T  ;<  Y 
C  CN  T  I NU  r 
CCS  7 iNUf 


8  2  3u 
6220 
c 

10050 


T  Ab  ViP  t  2*NP*N  i  »  144  ,C 


:  7AtVAPILAS7*l> 
z  7APVAP(LAS7*4) 
:  7|EVAp(LAS7*6J 


460.0 

3bCr.O 

p.n 


N=>  =  7  Ab  CRT  (  1  ) 

MR  Y  =  7  AbCR  7  (  “  J  » 
L  AST:?*  >*P 
DO  8 k 20 

NT:  7  A  b  C  -  T  (  2  *  ►.•  > 
7ABC«T( 2‘NP*N)= 
DC  6  2  3P  K  =  l,r.7 
7AbC»7(  LAStM  1  = 
LAS7=LAS7*N7AY 
CONTINUE 
C0N7 INUf 


T  AbCPT  I  2*‘lP*N>  »14tf.O 
TAbC&TILASTM  )*460.0 


8241 

8240 


C0N71NUr  . 

DO  3240  Ml  ,  NP7HA 

ir ( K6NS YP < N >  .LE  .C »G0  70  8240 

NODE  -  N  COE  S (  1  .N  ) 

LO  8241  uAl.NCOr 

1=  I  A  0  A  Y  (  1 ,N,L  ) 

j:  JARA  Y< 1 ,N.L  )  ,  , _  _ 

TIN(l»I.U>=  TIN*  1*1 *J»  ♦  4*>0.C  n 
HASID^IJ.J)  =  HASlDEUtJ)  /  3600. 0 
CONTINUE 
COM  INUF 

IF< ISDYA.NL. 1 )G0  70  1006C 


INITIAL  12  A  7  I  ON 

00  5000  1=1. M 
DO  5100  J  =  1 » N  J 
DO  5200  Ml. 2 
IF ( K .EC  .2  )  tlN(K 
TOUT  IK ,  I, J)  = 
TWALLtK.I.J)  =  ^.C 
T “E AN ( 4  ,1, J»  =  0.0 
CPHEANJ  4,  1  ,J»  =  C.C 
5200  CONTINUE  ,  n 

QDOT i  I . J »  =  C.C 

TMNCKlltJ)  =  r-o 

ABEEt!. Jt  =  t.r 


tI»J)=u.U 

C.C 
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T  E  R  u  S  ( 1 »I ,U>  =  0.C 

TERMS  < 2  , 1  .  J  )  =  C.C 

TERMS! 3,1, J)  =  C.C 
5  10 J  CONTINUE 

5  JQO  continue 

c 

10360  CONTINUE 
C 

c  •  •*  ‘"Initial  tempera  turf  assignments 

C 

DO  6 5 DC  N=1,NPTNA 
NODE-  N  CCE  S ( i,N) 

DO  6  510  L=1,NCDE 
irlARAY  (1  ,N,L  ) 
jrjARAY  !  1  ,N,L  ) 

IF(KANSTPIN)  .RT.OGO  TO  6520 
TIN(  1,1  ,J>=  T  2  R  0  (  1  , N  ) 

TOUT  t  1 ,  I, J  )  =  TZRO!  1  ,N ) 

TMCAN ( 1  ,1 , J  )  =  T  ZR  C  ( 1 , N ) 

GO  TC  6510 

6  52  C  CONTINUE 

TOUT  (  1 ,  I,  J»  :  TIMl.I.JI 
THE  AN  (  1  ,1 ,  J)  =  T  I N  ( 1 , 1  ,  J'/ 

6  5 1 3  CONTINUE 
65CC  CONTINUE 

IF (  ISDYN.NE .  1  ) G  0  TO  10070 
CO  6600  N=l,NPTHb 
NOCE  =  NCDESiZ.N) 

UO  6o 1 C  L-  1 » NODE 
I  =  IARAY(2,N,L) 

J  =  JARAy ( 2 , N , L ) 

T I N (  2 ,  1,01  :  TZR0(2,N) 

TOUT  12,  I , J )  =  TZD0 I  2 ,N ) 

TwAll(1,I,J>  :  TZ  D0 I  2 , N ) 

T  W  A  L  L  (  2  , 1 , J  >  =  TZRO  <  2 , N ) 

TCOKE ( I  , J )  R  T2»0«2,N) 

Tm£  AN  (  2  , 1 ,  U  )  =  TZR0I2*M 

ThkCk ( I  ,j  )  =  O.C 

6  fa  1  U  CONTINUE 
660C  C  ON  T I NU  r 
C 

IF ( N C K S  A V .NE • 2  >  G C  TC  6 7 CO 
C 

C  RECALL  COKE  THICKNESS  DATA 

C 

nhum;ma  xx a x  y 

CALL  NTRAN(NCKSTR,2?,10,22,?«12,TYTLE,LL»22, 

X  2 ,NDUM , THKCK ,LL , 22 , 10,22  ) 

WRITE  (KW,671UMTYTLE(I),  1=1, 1  2  > 

6710  FORMAT ( /IPX , ’COKE  THICKNESS  DATA  HAVE  BEEN  RECALLED 
X / 1 0 X ,  ’FROM  CASE ’ ,  IX, 12A6/ 1 
6700  CONTINUE 
C 

C  END  OF  INITIALIZATION 


IFINEPR  . L£  .  u  I  RETURN  1 
10070  CONTINUE. 

IF ( IS0YN.GT.2 1RETURN  1 
C 

C- - RESET  PRINT  OPTION 

C 

NPPNTRO 

C 

c- - COKE  THICKNESS  CANNOT  CHANGE  DURING  TRANSIENT 

C 

HOURS (  1  )  =  HOURS ( NHOURS ) 

NHOURS:  1 
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WRITE.  <K  W,  1U07D)  ELMASS 

10072  FORMAT!  //K  X,  ’ELMASS  =  • ,  E 12 .  S ,  IX  ,  ’LBM  •  /  > 
WRITE (K  w,  10074  ) 

10D7w  FORMAT ( //10X , ’ME TAL  HEAT  CAPACITY  TABLE’/ 

X  /  1  U  x  , *  N',  *  TEMPERATURE  CP’, 

X/ luX , 2x  ,  •  DEG  F  BTU/LBH-F • /  ) 

DO  1UH7E  N=1,NCPmET 

WRITE  «K  k, 1CG79 ) N,  TCP  HE  T (N  > ,CPMET  «N  » 

1CG78  FORMAT l  1LX , 12 ,2E12. E > 

TCPMET  (  M  t TCP *ET  (M  ♦460.0 
10076  CONTINUE 

WRITE  <KW, lCOoD) 

1.'  ZB'  FORMAT  (  ///) 

RETURN  1 
1G00C  CONTINUE 

ERROR  SUMMARY 

RETURN  " 

END 
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SUBROUTINE  C  A  LC  <  L  OT  T  A  1 
INCLUDE  PROCO.LIST 

DIMENSION  Z  T  A  P  (  1C)iT£RR0R(2.MAXPTH) ,QPATH(2,MAXPTK) 

DIMENSION  FILDAT (  10) .INVRSAIMAXY ,  M  A  X  X )  , INVRSB1 MAXY,MAXX  )  , 

X  OSAT(MAXPTH),UA!MAXYfMAXX),MATEPL(in),TRATt(MAXY,MAXX) 

DATA  (MATERL(N) ,N=1,1D1/  6HCARB0N,  6H  STEEL,  6HM0LYBD,  6HENUM 

X  6HSTNLS  ,  6HSTEEL  . .  '  - -  - - 

DATA  R  ?E  R 0  /  154  5,0/ 

REAL  NEWMOM 

COSB=COS( S»EEP*PI/18C. ) 


6HM0NEL 


6  H 


,  DHHUUTBU,  cmtniun  , 

,  6HSUPER  ,  6HALL0YS/ 


•SET  UP  INVERSE  ARRAY 


C 

C 

c 

c 

c 

c 

c 

c 

C 

c 


C 

c 

c 

c 


DO  t'u 
NODE=NODL$<  1  , N > 
DC  L  -1  ,  NODE 
I  =  IARAY  (1  ,N,L  > 
JZJARAY  (1,N,L1 
I NVRSA (  T  ,  J  >  =N 

20  CONTINUE 

00  21  N - 1 , NP  THb 
NODE=NOCES(2,N> 
DO  21  L  = 1 , NODE 
I=IARAY ( 2 , N • L ) 
JUARAY  U,N,L  > 
INVRbBl I, J)=N 

21  CONTINUE 


OUTER  L  COP  CONTROLS  AGING  OF  HEAT  EXCHANGER 

do  spncr  ISH0UR=1,NH0URS 

TIME 

THEHR  =  HOURS ( I SHOUR ) 

SET  COKE  THICKNESS 

IF(ISH0UR.LE.1)G0  TO  5QIDD 
DHCURS  =  ThEHR  -  HOURS < ISHOUR-1 ) 

IFIISH0UR.GT.2JG0  TO  50C50 

GET  RATE  OF  COKE  THICKNESS  FORMATION  BASED 
UPON  FIRST  SOLUTION 

DO  5U01C  1  =  1, NI 
DO  5QC1C  j: 1 , N J 

CALL  LOCK ( TCTAB ,THKCT ,NCOKE,TWALL (2. I» J) .THICK ,KK ) 

TRATE (  I  , JIRTHICK/HRSMAX 

CONTINUE 

CONTINUE 


50010 

c 

50050 

c 

c 


UPDATE  COKE  THICKNESSES 
DO  5UD6C  1=1, NI 
DO  5C06C  J=1,NU 
THKCK(I,J)=  THKCK (I ,  J  >  ♦ 
SO 06 Q  CONTINUE 


TR ATE ( I , J  )  *  DHOURS 


50100 

C 


CONTINUE 


IFINITER  . LT •  1)  NITER=25 
QNET=0. C 
THTAxO. C 

DO  1000C  ITER=1, NITER 

SAVE  PREVIOUS  PASS'S  OUTLET  TEMPERATURES 
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CO  1JC  mi.nptha 
NODE:  N  CD  C  S (  1  ,N  ) 

ILAST:  IAPAY I  1 , N, NODE  ) 

JLAST:  JAPAY(I.N.NODE) 

TOuTSV(l.N):  TOU T ( 1  , 1 L A S T  ,  JL A S T  ) 

10b  CONTINUE 

DC  110  K=1,NPTHP 
NODE  :  NODE  S I  2 , N  > 

HAST  :  1  ARAY  (2  ,N  ,NCCE  ) 

JLAST  :  JA  R  A  Y (2  , N  .NODE  ) 

TOUTSV(",N)=  TOUT ( 2 , IL AST  .JLAST > 

110  CONTINUE 

A-PATH  CONDITIONS 

DO  2CD  N-1,NPTHA 
NODE  :  NjDtS(l.N) 

DO  2 1 n  Li 1 . NODE 
I  =  I AR  AY ( 1 ,N,L  ) 

J  =  JARAYll.N.L) 

NP  =  INVR  SB  I  I  ,  J  > 

EL  :  DELTAX(N) 

IT ( ISTA PT ( 1 ,N  )  ,f Q.l  )  EL  =  DELTAYJNB) 

IF(S»EEP.GT.l.1)  GC  TO  111 
F  R  0  N  A  P  ( N )  :  2  A*EL 

F  L  AR A ( N  )  :  F  A  OF  A  (  N )  *  FRONAR(N) 

SRF  A  A  A ( N ) :  $AOV IN  >*ZA*DELT AX  IN  )*DELTAY (NB  > 

FINARIN)  :  FINSPF(N)*  SRFARA(N) 

GO  TO  112 

111  CONTINUE 

FAOFA(N)  =  l.C 
F  RON  A  R  (  f, )  :  EL*  2A  *  COSB 
FLARA  <  N  )  :  F  °0N  A  R ( N  ) 

SRFARA(N)  :  DELTAX(N)*DELTAYINe  > 

SAOV(N)  :  *  •  0 / 2  A 
FINSRF(N)  :  G.C 
FINAR(N)  =  L.C 

112  CONTINUE 
etaf  =  o.  r 
ETAZRO:  l.C 
RENiU.G 

HA  =  HASIDl(I.J) 

IF (KANSTP(N) .GT . C )  GC  TO  205 
C - THERMAL  PROPERTIES 

CALL  LOOkITEMAT,  VISAT,  NTABA,  THE AN (  1 , I , J ) ,  VIS,  KK) 

CALL  LOCK I TEMAT  ,  XKAT,  NTABA,  T ME  AN ( 1 , I , J 1 ,  XK,  KK) 

CALL  LOOKITEMAT,  CPAT.  NTABA,  TMEAN ( 1 , 1 , J  )  ,  CPMLAN(l.I.J).KK) 

c - TRANSPORT  COEFFICIENTS 

GEE=  WDOTll.N)/  FLARAIN) 

RTN  :  GEE*  DHYOll.N)  /VIS 
PRN  :  VIS*  CPMEAN(1,I,J)/XK 

Call  look*  renst,  stntab,  nstant,  ren,  stant.kk) 

HA  :  STANT  /  PRN**C. 66667  *  GEE  *  CPME AN  1 1 , 1 ,  J ) 

C - FIN  EFFECTIVENESS 

CALL  LOOK ( TUT AB * XKUTAB  ,  NWALK,  T W ALL I 1 , I , J ) ,  XKUA.KK) 

EM=  S0RT(2.0*HA  /  I X  K  W  A  *  FINTHK(N))) 

ETAF  :  TANH(EM*FINLEN(N)  )  / ( EM*F INLEN  C  N ) I 
ETAZRO  =  1.0  -  FINAR(N)  /  SRFARA(N)  *(1.0-  ETAF) 

205  CONTINUE 

FILDATl  1)  :  HA 

F ILDA  T ( 2)  :  ETAF 
FILDATl 3)  :  E  T  A  ZR  0 
FILDATI4)  :  REN 

TERMS ( 1 .I.J)=1.0/(ETAZR0  *SRF ARA ( N )  *  HA) 

LREC  =(  1-1  )*NJ*J 

WPI TE I N SAVE’LPEC) (FILDATl  KKJ  , KK : 1 , NUM WR D ) 

IF(  ITER.GT.l  )  GO  TO  210 
ONET  =  ONE  T  *1 ,0/TEPMS ( 1 , 1 , J  )  •  ( TMEAN ( 1 , 1 , J )  - 
X  0.5*<T«ALL(1,I,J>+TWALL<2.I,J>>  > 

THTA=THTA  ♦  SRFARA(N) 

210  CONTINUE 


200  CONTINUE 
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6-PATH  CJNDITIONS 

CGAIf.r  '  ,i 

DC  3C.C  K=1,NPTH9 

NODE  :  N0DES(2»N> 

X  -  G  *  0 

X  C  K  -  u « C 

H9V=U.O 

Q  S  A  T  ( f  j  )  r  .  _ 


heat  addition  fob  constant  tcnpebatobe  vapobijatiot  ppcce 
o!.1vlapl=SSI(^T!5;5)!5,'ts’  »a-po.2,n,,  h«ap, 

DO  3D  I  - 1 , M 
DO  3  0  J - 1 »  N  J 
wUAL  (  I,  J>  =c  ,L 
C  C  NT  I  N  u  E 

-SATURATION  TEMPERATURE 

D0L3lfc°LrKNSJETD,TS4TTiJ’  NSM»  p^0«2.N),TSAT(N),  KK) 

I  =  I  AR  AY  (  2  ,  N  ,L  ) 

J  =„<J*R  4V«.  .N,L  ) 

NA  --  IKVPS4(I,J) 

EL=DELT  ax ( na  j 

jf JiSrEpIAf :SIi >  =  *le  TgL!?ELT'"' IN* 

p  D  AP  B ( N  )  ;  NTUBES(N)*PJ/4,D  *  f  H  Y  D  (  2  N)*#? 

Ji  •  i0M?g!!;S!:;5KN«u.EL 
GO  TO  irN  '  ‘,TuBES(M*  pI  *  DhY0(2,N)*EL 
CONTINUF 
EL  =  DLLT  A  a ( NA ) 

NTUBES?nI=1,N,*EJ*I>  ELrDELTAYJN) 
nhvn?5\‘,r  EL*ZB*COSB 

*  ^ I u  '  *  -  *  •  Q  *  EL*  ZB  *  C  0  S  B  /  (FL*rn^R*7a» 

•  ^J.DDjl*N)  =  DEL  T  A  X  I  N  A  )  *DEL  T  A  Y  (  N  ) 

A«ALU2,M  =  A  WALL  1 1 , N ) 

CONTINUE 

IF(P2R0  (2,N)  .GT.PORITBI  GO  TO  311 
■GET  QUALITY 


CONf?Nllt?il’:UJ’’iS'>U,!E 
distillate  substance 

CUAui?S)=P2TSlKAHPZP0,‘,N,*™EAN<2*r'J,*2TAe*K'<» 

GO  TO  333 
CONTINUE 

PURE  SUBSTANCE,  QUALITY  DETERMINED  FROM  PREVIOUS  NODE. 

IFIL.EQ.1)  GO  TO  323 
IP=IARA Y( 2»N, L-l ) 

JPRJAPA Y( 2  *  N , L -  1  ) 

8y£Vliil^,=oui,-,IP--JP» 

(I , J)  ,LE .C.O  )  GO  TO  311 

ihqual  u»j),cr,a.o  .and.  oualii.jj.lt.i.oj  go  to  312 
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C . FLOW  IS  1CGX  VAPOR 

C 

IFLAG(I.J)  =  T 

* . CALLML0ChSpiTA&VAP,P2R0(2,N),TkEAN(2,I,J),2TAB,kk) 

XK  -  Z  T  A t  <  3  > 

VIS  =  ?  T  Afa  i  ?) 

HtiiCB  =  Z  T  AB  t  1) 

CPMEANI  2,1  i  4  > 

C . TRANSPORT  COt^PICIENTS 

GrE  R  VDOT t 2 ,N ) /FL APB M > 

KEN  ^  ^c.  E  *  uHYQ<r,K')/VIS 

?5!;u:LO^:C^i:siS:IUlJB:NSTKNtB,Rt5.ST.M,K«, 

HB-  STANT*  GEE*  CPMEAN I  2  *  I , J >  /  PRN**  -  •  fc66t>  7 
HCONV  R  ME 
H  0  U  R  o  R  C  .  L' 

HP  RE  A  L  R  Lb 
GO  TO  3  1  .. 

311  CONTINUE 


C 

C 

C 


C - 


Flow  is  lcm  licuID 

IFL AG  < I  ,  J  1  R  i 
■thermal  properties 
call  LOCK ( TEMPT , 

CALL  LOCKITEMPT 
CALL  LOCK ( TE  MB  T 
CALL  LOCKITEMPT 
CALL  LOCK  (  TE MET  ,  • 

--TRANSPORT  COEFFICIENTS 
GEE  R  -CCTI2.N)/ 

RENR  GEE  *  uHYDI 


E  MP  T  \  1  ^  V  I SBT  ,  NTAPE,  T ME  AN l 2 » I  * U  )  ,  Vl^,  **{ 
rCoT  JjjvT.  N  T  A  0  P  •  T ME  AN ( 2 ,  I  « J  )  ,  XK,  KK) 
FMPT*  c  P  b  T  ,  N  T  A  p  fc  ,  TME AN  I  2 , 1 , J  )  »  CP"ijljl2i!t 
EMPt)  RHOB  t  ,  NT  ABB  ,  TMEAM?,i,J),  ^  } 

FMBT,  RHCBT,  NTABB,  TCOKEU.J),  RHOV ,  KK) 

r  '-.r  ffTCTf^U 


J) 


FLARB IN  ) 
■’,N>  /VIS 


p»n; 


C  PMC  A  N l 2,1 , j)  /  XK 


LICCOR 


C 

c- 


34  1 


345 


TO  340 
THAN  PC 


C- 

c- 


3?g 


34G 


GC  TO  I  74 1  ,  342  > , 

CONTINUE 

IF ( P2R0  12, N> »GT .  pCRIT3>G0 
-T  LCSS  Than  TSAT  /  P  LESS 
IFISEN.GT.23Cr.n)  GC  TO  32 U 

CALL^LCCkI TEKBT »  VI SE T , NT  A t B ,  T COK E 1 1 , J J *  V I SU ,  KK) 

ANURALA“*IREN*PPN/IX/0HYD(2,N)  )  )  *  *  t)L  AM 

-MINIMUM  ^lAMINAR^NUSSELT  NO.  IS  FOR  FULLY 

’iFUNU^LT  .3^66)  ANUR3.6E 

HE  R  AN  L*  XK  /  OhYOI2,N) 

HCONV  R  h8 
HGURG  R  u.C 
HBREALRHb 
GO  TO  3  To 
CONTINUE 

A  NU  =  ATUR3*REN«* 

HB  =  AND*  XK  /  0HYDI2.N) 

HCONV  :  ne 
HGURG  = 

HEREALRMu 
GO  TO  3  T 
CONTINUE 

■P  greater  THAN  PC  . 

IFIREN.LT. 23LC.0)G0  TO  345 
CALL  LOCK  I  TEMPT ,  VlSbT,NTABB, 

CALL  LOCkITEMPT,  XKBT,  1 

LOCKITEMPT,  CPBT,  NTAfaB, 

GE  E*Oh YD  I  2 , N 1/VISW 
VISW*  CPW/XKW 

AS  UP*  REN**BSUP*PRN**CSUP*(PHOW/RHOS)**DSUP 


BTURB  *  PRN**  CTURB 


CALL 
CALL 
REN  = 
PRNR 
XK  R 
ANU  = 


TCOKE  C I , J ) , V ISW  » 
TCOKt (I,J» .XKW, 
TCOKE (  X  *  J  )  tCPW , 


KK  ) 
KK  ) 
KK) 


,  KK  ) 
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H3  :  AN  U*  XK  /  DHYD(Z,N) 

HCONV  =  HB 
HGURG  :  0.0 
HBRE  AL- Ha 
GO  TO  37y 
J42  CONTINUE" 

C - LIQUID  STANTON  NUMBER 

CALL  LOCK  ( RE  N  L I Q  »  $TNLICi*NLJAY,REN,STANT  ,  KK  ) 

HB  =  STANT*  GEE*  CPME AN < 2 , I , J )  /  PRN  **  C. 66667 
HCONV  =  HB 
HGUPy  -  u.L 
HBPEAL-Hj 
GC.  TO  3  7J 
312  CONTINUE 

- FL0»  IS  A  MIXTURE 

CALL  L0CK(PMIXTE,CPMIXb,NMIXE,P2R0(2,N>,CPMEAN(?,I,J),NK) 
IF(K0MPLX  ,GT.  0)  GO  TO  375 
HB  =  1 »  Cc.  ♦  20 
HCONV  =  mR 
HC-URb  =  G.D 
HCREAL=r.D  , 

I  F  L  A  G ( I  ,J)=? 

GO  TO  370 


375  CONTINUE 

----TWO  PHASE  BOILING 

IF(QUAL(I|J)  .GT,  0.7)  GO  TO  371 
GEE  =  LOOT i 2  .  N ) /  FLARy(N) 

CALL  LOOK ( TEKBT  ,  VISPT,  NT  ABB  .  T ME  AN < 2 , I ,  J  >  ,  VISL,  KK) 

CALL  LOCnUEMBT  ,  •  XK*8!  ,  NT  ABB  *  THE  AN  (  2 , 1 ,  J  )  ,  XKL  »  KK) 

CALL ■ LOCK < TEHBT  ,  CPRT,  NT  ABB  ,  T M£ A N ( 2 , 1 , J > ,  CPL,  KK) 

CALL  LOCK) TEMBT  •  RHCBT,  NTABB,  THE  AN ( 2 ,  I  , J > ,  PHOL  ,  KKI 

CALL  LOC.v  (  PLAmte  ,  HVAPT3  »  NVAPTfa,  PZR0»2,N),  HVAP,  KK  ) 
call  lock j i sigma , sigtab,  NSIGMA,  THE  AN  <  2  « I . J  )  .  SIGMA, KK) 

CALL  LOCK) TSATTR.PSATTB.  NSAT,  TME AN ( 2 , I , J ) ,  PSTT8,  KK) 

CALL  LOCK ( TSATTB.PSATT6,  NSAT,  TCOKE(I,J),  PSTTW,  KK) 

CALL  LOCkUPI TABVAF,  PZR0(2,N),  TME AN ( 2 ,  I , J ) ,  ZTA&.KK) 

RHOV  =  Z  T  AS  l  1 ) 

V ISV  =  2  T  A  B ( 2) 

V  F  G -  1 . / RH 0  V  -  l./RHOL 

REN  =  GEE*( 1 .O-QUAl ( I ,J ) )  *0HYD<2,N>/  VISL 
PRN  :  VISL*  CPL/XKL 

XFUNCT  =  (CUAL(I,J)/(1.0-0UAL(I,J) ) ) **Q • 9 

X  *  (RHOL/RHOV  )**0.5  *  ( VISV/VISL  )**0.1 

CALL  LOCK ( XOVFTB ,FOVFTB,  NOVXF,  XFUNCT,  FVAL,  KK) 

HCONV  :  0,023*XKL  /  0HYD«2,NI  *  RCN**0.8  *PRN**0.4  *  FVAL 

SARG  =  REN*  FVAL**!, 25 

CALL  LOCMSRtLTB,  STAB,  NSTAB  *  SARG,  ESS,  KK ) 

HGUR a  -  C  •  w 
OT  :  0. C 

OP  =  o.n 

IF(TCOKF( I,J)  .LT.  T  S  A  T  <  N ) )  GO  TO  372 
OT=TCOKE< I . J )  -  TSAT(N) 

DP=778.16*0T*HVAP/(TSAT(N)*VFG) 

HGURG  =  0,00122  *  XKL**0.79  *CPL**0.45  *RH0L**0.49  *  GC**C.25 

X  / (SIGMA**0,25  *  V ISL**0»29  *  HVAP**n,24  *  RH0V**0,24) 

X  *  0T**0,24  *  0  P**0 » 7  5  *  ESS 

372  CONTINUE 
HMAC  =  HC  CNV 
HMIC=HGLRG 
HB  =  HC  CNV ♦HGURG 
HBREAL=Hd 
GO  TO  3  7 J 
371  CONTINUE 
C 

C . HIGH  QUALITY  MIXTURE 

C 
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IF  ( MB  V  .GT.  0.0  GO  TO  373 

CALL  LOCKUP ( T ABVAP,PZRO (2,N) ,TMEAN<2, I, J> ,2TAB,kk) 

XK:ZT  At (31 

VIS-ZTA912) 

RwOBpZT  At (  1  ) 

CPV  =  2TAP<  k  ) 

GEE=*DO  T  t  2 , N  > /F  L A&b  <N  ) 

R E N  =  oE E  *l'HYO  <  ?, N  )  /V  IS 
PBfc=tfIS*CPV/XK 

CALL  LG  CM  REN  STB  ,  ST N T 8 , NSTNTe , RE N , ST ANT  ,KK  ) 

HP  V  =  STA  M *GEl*CPV/PRN **0.66667 

37  3  H?0NV=H&V  ♦  I  <  1  .-OU  AL  1 1  ,  J  )) /C.  3)**0.S*«  HMAC-HBv  ) 
HGuPGRH^IC**  ( 1.0-CUAL  (I i J)  )/C  • 3 ) **Q. 5 
HB  =  hCONV  ♦  Hul.RG 
hBRGAL=  PL 
C 

37  ij  CONTINUE. 

FILOATI  1)  :  Hf 
FILOAT ( 2)  =  R  r  N 

FILOATI  U)  z  hCGNV 
FlLOAT(S)  =  hCUPG 
C 


C . COKE 

C 

T  hk  C  =  TnKCKd.J) 

C - TERMS  IN  EFFECTIVE  HE»T  TRANSFER  EQUATION 

TERMS(3,I,J)  =  ThkC  /(XkCGKE*  Awall(2,n>) 
AeEE«I.J)=NTuPES»N)*PI*EL*IDHYDl2,N)-THKC> 
IF|kAN$TP<nA  )  .GT.C)  A  BE  E  I  I  »  J  )  =  SRFARA(NA) 
TrRuS(4  ,1,  J)  :  l.G  /  <Hb*  ABEEU.J)) 

TRAP:  r.s»»TWALLU.l,J)*TWALU2,I,jn 


CALL  LOCk  (  TUT  A6  ,  XKWTAB  .NtaALL  i  TE  Afi  ,  AKW  ,  KM 
TERMS  (2  ,  I  ,  J  I  =  THKUAL  /(XKW*  AWAIL(1,N)> 

C 

c - HEAT  TRANSFER  RATE  ( BTU/  SEC) 

UA(  I,  J)  =C.t» 

DO  3oC  L^=1,m 

UA(I.J)  :  U  A ( I , J  )  ♦  TERMS ( LQ , I , J ) 

38  C  CONTINUE 

FILOATI  3)  1.0/  UAIIiJ) 

CALL  LOCMTEMAT,  CPaT,  N  T  A  9  A  ,  TIMl.I.JI.  CPAIN.kK) 

CALL  LOCMTEMAT,  CPAT,  NTABA,  TOUT  (  1 , 1 ,  J  )  ,  CPAOUT.KK) 

COOT  I  I  ,  J)-(TIN(  1  *I»J)-TIN<2«I  ,J)  )/ 

X  ( U  A ( I  ,  J  )  ♦  l.t/(2.C*n.D0T?i  ,NA)*CPMEANfl,I  ,J> ) 

X  ♦ 1.0/C 2.C*UOOT(2,N)*CPMEAN(2,I , J) )  ) 

TOUT! 2,  I. J)  =  TINI 2,1 »J>*CDOT< I » J>/tuOOT ( 2 ,N»*CPME AN«  2, 1 , J>1 

GC  TO  I  3,0,392  ) ,  1SPURE 
394  CONTINUE 

IF  UJOT  U  ,NA  »  .LE.C.r  .OR.  UOOT  (  2 , N ) • LE • l • H )  QDOT  <  I ,  J  )  =C  .C 
TCUTI2, I,J)  =  TIN(2,I,J)*30OTIl,J)/(*OOT(2,N)*CPMEAN(2,I,J)  ) 

IF ( IFLAGI I , J  )  ,E0. 1  .AND.  T OUT « 2 , 1 , J ) ,GT  .  t S AT < N  )  > 

1  TOUT  (2,!«J)-TSAT (N) 

IFIOSAT IN) .GT.C.C  , AND •  Q S AT ( N  )  .LT . OH V A P )  T OUT < 2 , 1 , J ) = T S A T 1 N  ) 
GGAIN=CGAIN  ♦  QDOT ( I , J ) 


C 

C 


C 

c 

c 


391 


CONSTANT  TEMPERATURE  HEAT  ADOITION  PROCESS 


IFITOUT  (2*I,J)»LY.TSAT|N)  ,0R,  XCK.r,T.2.0)  GO  TO  390 

1FWAG<I  .J)=2 

TOUT! ?, t, J)=TSAT<N) 


ACCUMULATED  HEAT  TRANSFER  INCLUDES  CORRECTION 
FOR  FIRST  NOOE  THAT  JUST  CONTAINS  MIXTURE 


UNO  TOUT  JUST  EXCEEDS  TSAT) 

CSAT  »N>  =CSAT (N)*QOOT<I , J ) -LOOT  I  2, N I *CPL* < T S AT < N > «T IN <  2, I, J| » 
IFiCSAT  ft.)  .GT.O.G)  GO  TO  391 
IFL AG  ( I  *u  )  =  1 


QSATIN)  =C.(i 


CONTINUE 


QUAHI,  J)=OSAT«N)/QMVAP 
IF(QUAL(I,J).GT.1,0)CUAUI,J)S1.0 
IFJQUAL (I,J)  .LT.i.GJGO  TO  590 
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- FLOW  HAS  JUST  BECOME  VAPOR 

TSATtN) =TS*T IM1 *.C1 

cc^hLN??^:srrBz^g!^po*2’Ns*Ts*T,N»»2T*a*^» 

OSAT(N)  -wHVApSAT<NiM,QSAT<N,"°HVAP)/<yD0T  (2*N,*CPV) 
IFLAi ( lTj)=3 
X  C  K  :  3 , 0 
?9J  CONTINUE 


IF IL.Eti  .NODE >G0  TO  7E  1 
IP  =  IA5AY(2,N,L*1) 

JP  :  JARAYl 2  »  N  *  L  ♦  1  ) 

T I N I  2,1 F, JP) =  TOUT ( 2, 1 , J) 

TMEAN(2»I':-tJP):TIN(2,IP,JP) 

381  CONTINUE 

IREC  :  i,I*NJ*  (  I  - 1  )  >»N  J  ♦  J 

T?fIrlnS^E:LPEC>  (FILD4T(KK),KK  =  1,NUMwRO) 

IF  (ITER. E 0*1) 

31UXC0NUNUFWNLT  ’  HSRL  a>-*A9EE  (  I ,  Jl  *  {  T  B  AP -THE  AN  (  2 , 1 ,  J )  1 
SCu  CONTINUE 

UPDATE  TEMPERATURES 

QLOST  :  u.O 
DC  4D0  N:l,NPTHA 
QPATrtd  ,  N  )  :  Q  ,  C 
NODE:  NODESH  ,N) 

DC  41P  L: 1  *  NODE 
1  :  I AR  AY  1  1  , N  ,L  ) 
j  ;  jaray (  1,N,L  ) 

IFl -DOT  1 1 ,N)  .lE.  (j.L)  TIM  1  , 1 ,  Jl  =TIN(  2,  I ,  Jl 
IF(KAMSTPiN) .GT.C)  CO  TO  301 

$pat,  ntaba,  tinu.i.j),  cpa:n,kk> 

I  c  U  T  *  i,  I,  J  >  -  T  IN  *  1 ,1  ,J)  -  QD0TU»J»/(WD0T(1,M*CPAI;«> 

302  CONTINUE 
13:13*1 

T0T1J:TCUT»1,I,J) 

CALL  LOCK  J  TEMAT •  CPAT,  NTAB A  ,  TOTIJ.  CPAOUT.KK) 

T  r  Vi  A I J  tAH  t  T  i  C?  ^  1  t*^T!1i!I,Iji“^90I<*lJ>/yD0TI  ^  JN  *  l/CPACUT 

.LT.  1.0  .OR.  I3.GT.5)  CO  TO  3fU 

6 C  TO  3  Ci 
301  CONTINUE 

^CPATHU  ,N>:QPATH»  1,N>  ♦  WOOT (1 ,N )*( CPAIN*T INI  1 , I , J) -CPAOUT* 

TmEANI1,I,jI:,C.S*<TIN(1,I,J)*  TOUT  I  1  , 1 , J) > 

TWALLU.I.J):  TME  AN  (  1 , 1  ,  j  )-  QDOT  I  I  ,  J  )  *  TEWMSt  1«  I,J> 

TWALL  2  tl.J):  T W ALL ( 1 . 1  * J ) -  QDOTII.J)  *  TEHMSII  I  j) 

TCOKt  (I  ,J)  :  T  WALL  t  2 1 1  *  J )  “  QDOT  ( I !  J  »  *  TERMSI3  I  J) 

IFIL.EC  .NOOE .CR.KANSTPIN) .EQ.l >G0  TO  410 
IP:  I A R A Y  1  1 , N f L ♦  1  > 

JP:  JARAYl  1  ,  N , L ♦ 1 ) 

TINI1 ,1 P, JP)  :  TOUT!  1 ,  I ,  J) 

4 1 J  CONTINUE 

QLOST:OLOST  ♦  QPATH (1«N  ) 

400  CONTINUE 

DO  420  N= 1 i NPTHfl 
OPATh ( 2  ,N  )  =o.n 
NODE  =  NODES! 2 .  N ) 

DO  430  L:  1  .NODE 
I  :  IARAY I i, N, L ) 

J  :  JARAYl 2 , N ,L ) 

QPATH I  2  ,N ) :QPATM ( 2, N)  ♦  0 DO T 1 1 , J ) 

430  CONTINu£I,J,=  0*5  *<TIN<2»I*J)  4  TOUT  I  2  *  I , J )  ) 

420  CONTINUF 
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C 

C  CONVERGENCE  CHECK 

C 

N0G00P  =  0 
DO  5uP  N- 1 »  N P  TH  A 
NODE  =  F.UDESCl.N) 

I  :  1AR AY ( 1 , N .NODE  ) 

J  =  JAR  ay <  1 ,  N  , NODE  ) 

TChK  r  TOUTll.l.J)  -  TOUT  SV ( 1 » N  J 
TERROR (  1,N)=TChk 

IF  (  AfcS(  TCHk  )  .GT  .  TOLITR)  NOGOOD:  NOGOCO*! 
5Cu  CONTINUE 

DC  5ir'  N=1,NPTH0 
NODE  :  NUDt  S l 2 , N  ) 

I  :  IARAV(2»NtM0DE) 

J  :  JAR  AY ( 2 , N  .NODE  ) 

TCHK  =  TOUT (2,1 ,  J )  -  TOUT  S V ( 2  ,  N  ) 

TERROR ( 2, N ) =Tf HK 

IF ( AbS( TCHK  )  .GT  .  TOLITR)  NOGOOL  =  N0G00D*1 
5 1 C  CONTINUE 
ITAT=  ITER 

IF (NOOCCo  «LE  »C )  GO  TO  20  0CC 
C 

IF  (NUU*'P.LE.b)GO  TO  1Q0PD 


C 

C  DUMP  OF  iNTEh“ECI  ATE  RESULTS 

C 

NALL  =  r 
DO  2U0C  I D - 1 » N I 
DO  2050  JD=1,NJ 

nallrnall*: 

IFimcDI  NALL-1. 50  .NE.OGO  TO  2100 
»>RITt  IK  W,  211C  )  ITAT, NOGOOD 

211U  FOR“AT(  lnl./lGx, ‘DUMP  FOR  ITERATION* , 15, IX NOGOO": *, IS  ,  / 
X/iCX,‘  I  J  TA-IN  TA-OUT  tA-HEAN  T.ALL-A*. 

X  •  TwALL-P  TCCKE  TB-IN  TB-OUT  Td-MEAN*, 

x  *  crcT’, 

X/2DX.  9  <  5  x , ' DE  G  0  * ) , 3X  ,  ' STU/SEC  * / ) 

2  IUli  CONTi  NU  r 


212C 

2U5U 


WRITE  IK  w.  212U  )  I u  » JC ,  TIN(l.ID.JD)  .  TOUT  1 1 .  ID ,  JD  I  » 

<  Tv-1  AN(  I  ,ID,  JD)  ,  I  TWALL  (L  ,ID,  JD  )  ,L  =  1,2)  , 

<  TCCKE ( ID, JD  ),  TIN(2,ID,JD)»  T OUT  I  2 , 1 0 , JD ) , 

r  Tmean<2,id, JO),  GDOTtID.JD) 

FORMAT!  10X,2I5,9Flw.3,Ei0.5> 

CONTINUE 
CONTINUE 

nallrd 


2GC4 

) 

) 

) 

2  UO  3 

: 

2J0s‘ 

2J02 

2001 

) 

e213C 

1000C 

c 


ALLRH 

0  2Grfc  IC=1,NI 
iO  2CP2  JD=1,NJ 
i  ALL  =  NA  LL ♦  1 
A=INVRSA ( 10, JD  ) 

BriNVRSb  t 10, JD  ) 

FIMOD! NALL-1,50)  .NE.D)  GO  TO  2003 
RITE (K  W, ?  ,0m  )  ITAT  , NO GO 00 

ORMAT!  1H1,/IX, ’CLMP  FOR  I  TER  A TION* , 15 . IX , • NOGOOD= * . I  5 , / 
/IX,*  I  J  CPKEAN-A  CPMEAN-3  TERM-1  TERM-2* 

»  TtPM-J  TEPM-4  UA  QUALITY  deltax*, 

•  ClLTAY  ?  IFLAG*/) 

ONTINUE 

RITE«KW,2i DS>  IO.JD, ICPMEAN(I1,ID,JD),I1=1,2) , 

(TERMS! I r.IC.JC)  ,12  =  1,4)  ,UA!I0,JD) ,QUAL (10, JD  )  , 
OELTAx(NA),OElTAY(NB),IFLAG(IO,JO) 

ORMAT  (  1X,2I5,?E1C.5,3F  IC.4,110 
ONTINUE 

R I TE ( K  W , 2G0 1 )  OLOST.QGAIN 

ORM A  T (  /////lrX , 'HEAT  LOST  FROM  A-SIDE  =  *,E11.5,*  BTU/SEC* 
/ IPX , ' HE  AT  GAINED  BY  B-SIDE  =  *,EU.5,*  BTU/SEC* 
TE(KW,?13C) 

4  A  T  (  // /  ) 

TTUiir 


CONTINUf 
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C  NO  CONVERGENCE 

WRITE  IK W, 600 >  ITAT,  NOGOQD 

600  FORMAT!  ///IDa, 'FAILED  TO  CON VE RGE . . . A  FT ER  *  ,  1 5  » •  ITERATIONS,', 
X  IX, 'NOGOOD  =  ',  110//) 

GO  TO  3  ruCu 
2CwCC  CONTINUE 
C 

C  CONVERGENCE 


61  . 
3Du2G 


602 
11  JQU 


12G0D 

C 

C 

C 

c 

C 


WRITE(KW,610>  ITAT 

FORMAT! ///10x, 'PROGRAM  CONVERGED  AF T £ R * , 1 5 , 1 X ,  •  I T E R A T  I  ON S * / / ) 

CCNTINUF 

WRITE  <K  W , 6G 1 > 

FORMAT (// 1DX , *S IDE  PATH  ERROR-DEG  F  * / * 

DO  llCOT  N=1,‘:PTHA 
LABSDrus-  A 

WRITE  <KW, 602) LABSD.N, TERROR « 1,N) 

FORMAT (  1CX, A4,I5,5X,F10.3> 

CONTINUE 

do  wncr:  n=i,npthb 

LABSD=4H  0 

WRITE !KW,6u2 )LABSD,N, TERROR (2,N) 

CONTINUE 

PRESSURE  DROP  CALCULATION 
A  SIDE 

DO  3irOr  N=1,NPTHA 
NODE  =  NODES! 1 .  N  ) 

D  c  L  P  I  1  ,M  =  CUD 
HEADNUR'.r 
DO  3IlCr  L - 1 .NODE 
hEADGLRMEADNU 
I"  IARA V « 1 ,N,L ) 

JZ  JARAVI 1 ,  N  ,  L  ) 

ROIN  =  PZRO !  3 ,  N  )  *AMUA  /  (RZERO  *  TIN(1,I,J)3 
ROEX  :  PZRO !  1 ,  N  ) *AMUA  /  (RZERO  *  T0UT(1,I,J)) 

LREC= ( I -1 >*NJ^J 

READ(NSAVE"LRFCMFILDAT(KK)  ,KK  =  1  .NUMwRO  ) 

REN=FILCAT(4» 

CALL  LO  CK ( HE  NF , F  T  Afa  ,  NFRIC ,REN ,  FRIC,KK) 

GEE  =  WD 0 T (  1 ,  N  )  /  FLARA(N) 

DELP(1.N>=  DELP(l.N)  ♦  GEE**2  /(2,0*GC  *ROIN>  .  „ 

X  *{  <1.0*  (FLARA(N)/FRONAR(N) 1**2)  *  (ROIN/ROEX-l.C) 

X  ♦  FRIC  *  SRFARA(N)/  FLARA (N)*ROIN  /  ( 0 . 5* ( R 0 1 N ♦ROE  X ) )  ) 

HEAONU  :  GEE ** 2/  <  2 . T*  GC  *  0 , 5* ( RO IN  +  ROE X )  ) 

IF(L.EU.1)G0  TO  3  1 1 OC 


•TURN  LOSS 


KSCORE=C 

IM1=  lARAYi i,N,L-l> 

UM1=  JARAY(  1  ,N',L-1» 

IFIISTARTI l,N).EQ.2)G0  TO  3112D 

C - “HORIZON  TAL  TRAVERSE 

IF(NJ.EC.2)G0  TO  31 100 
IFIJ.EO.l  .OR.  J.EQ.NJ)  KSC0RE=1 
IF(JM1.E3.1  .OR,  JMl.ECl.NJ)  K SCORE =K SCORE ♦  1 
GO  TO  31150 
31120  CONTINUE 

C - VERTICAL  TRAVERSE 

IF (N1.EC.2)  GO  TO  31100 
I F ( I  . E  C. 1  «  0°.  I.EO.NI)  K  SCORES  1 
IF (  IMl  • Ew*  1  .OR,  I M  1 , E  0  .  NI)  KSCORE=  KSC0RE*1 
3115U  CONTINUE 

IF(KSCORE.EO.Z)  DELP(1,N)  =  OCLP ( 1 • N 1  *  C .5* ( HE AONU*HE ADOU *TURNLA 
31  IOC  CONTINUE 
C 


126 
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:iuuu 

c 

c 

c 


v  v  i 

6  DiOE 

DC  32GC  i"  N=1,NPTh3 
NODE  =  NOCLS(  2,N) 

D  E  L  P  (  2  ,  N )  =  L  .0 

i .  p  a  n  k.  1 1  — 


r.o 
c.o 

L=1 .NODE 
HE  4  D  MJ 
NE  w mom 


W  U  l.  '  I  fc 

HEAONU 
N  E  w  M  0  M  : 

DO  329DC 
HEADOL  : 

OLDMOH  r 
1=  I  A°A  Y< 2,N,L  ) 
J-  JAPAY(2,N,L> 
N4=  INVR  SA (I , J  > 

EL  =  DEL  TA  X ( NA  ) 
I  F  <  ISTAPT  <  2,N  >  .E 
iriS.EEP.GT.L.l ) 


.  - -  •  . w.2  )  EL  :  DELTAY(NI 

iriS.EEP.GT.L.ll  GO  TO  32001 
FLAPb(N)  r  NTUBES  (N  »*PI /4.L  *  DHYD(2,M**2 
ALALLd.M  =  NTUBES(N)*  PI  *  (  DMYD  (  2  ,  K  )  ♦  Thk  U  til  *EL 
A  «  A  L  L  <  2  » N  )  r  NTUEESIM*  PI  *  OhYO l 2 , N > *EL 
GC  TO  3  2oC2 
32001  CONTINUE 

ELADELT4a(NA  ) 

IF(ISTACT(2,N).EC.l) 

NTU3 S  ( M  :  1 
F  L  A  R  b  CM  )  =  EL*2S*C0S? 

OHYO  (  2,M= 

A W A L L  (  1  |NI 
A  WALL (2  ,N> 

32J02  CONTINUE 

ALPHA ( I  .  J  ) 


EL=DELTAY«N» 


.2*  EL*  2B  *COSB  /  (EL*CCSB*Zo> 
DELTA  X  (NA  )*DELT AY ( N  1 
AL ALL  I  1  ,N ) 


32100 

C 

c 

c 


KGO  =  I Fl AG ( I , J  ) 

GO  TO  (32  ICG,  322CC,  32  3D0), 
CONTINUE 

FLOW  lb  1CCX  LICUlu 


KGO 


32105 


GCC  =  wCOT(  ?,  M/FLABMM 

CALL  LOCK ( TEPBT  ,  VIS9T  ,  NT  ABB  ,  T M£ AN ( 2 , 1 , J >  ,  VISB.Kk) 

CALL  LOCk  (TEHBT  ,  VISBT,  NTAbP,  TCOKEU.J),  VISH.KK) 

1F(P2R0  I2.M.CT.PCRITB  .AND.  T  PE  AN  (  2 , 1 ,  J  »  ,GT  ,  TCRITB)  GO  TO  32 
call  LCCk ( TEHBT ,RH05T  ,  NT  ABB ,  THE  AN ( 2 , 1 , J >  ,  RHGB.Kk) 

GO  TO  3->llC 

r  flNT  T  Nil  F 


i  K  K  ) 


u  u  i  u  .i  ■  u  i 
CONTINUE 

call  L0CkUP(TABCRT,PZP0(Z,N),TMEAN(2,I»J>,ZTAB,I 

RH03=ZT Au ( 1 > 

32110  CONTINUE 

-  :  GFE*  0HYDI2,.,,, 

.  .  -  -  ---  PEN,FRIC,KK ) 

'  DHVD(2»N)*GEE**2  /(?.r*GC*RH09 > 


nnuu-t  »  *■  U  «  A  9 

CONTINUE 

REN  =  GFE*  DHYDI2.N)/  VI5B 
call  LOCKCRENFb  ,FBTAE,  NFPb,  S 
DELT  AP (  1 , 1 , J  >  =  4.C*  FRIC*  EL/ 

X  *  (VISU  / V I S  3  )  **  0.14 

NEWHO*  :  C.O 


3220C 

C 

C 

C 


NEwho*  :  C.O 
GO  TO  J?bOU 
CONTINUE 

MIXTURE 

GEE=  WOCTI 2,N)/FLARB(N> 
■LICUIC 


T IN ( 2  *  I , J  )  « 
T IN ( 2 , 1 , J  )  , 


VISL.KK) 
RHOL , KK  ) 


LICUIC 

CALL  LOCK ( TEMBT  ,  VISBT  ,  NTA&B, 

CALL  LOCK < TCMCT ,RHOBT ,  NTAbP,  T1N(2,J,J>,  RHOL  ,  K  K 
VAPOR 

CALL  LOCKUP ( TABVAP,PZP0(2,N)  , TOUT (2, 1, J) ,ZTAB,KK ) 

RHOV  :  ZTAfe(l) 

VISV  :  ZT  A b ( 2  ) 

=  (l.G  -  CUAL( I . J) >*GEE*  OHYD ( 2 • N ) / V I SL 
C  U  A  L  I I ,  J  »  *  GEE*  QHYD(?,N»/VISV 

t N V, RENL,  RHOV,  RHOL, OUAL( I , J » , B IGX ( I , J ) , PS  I V2 , PSIl 


V  "  w  s 

RHOV 

VISV 

RENL 

RENV 

CALL 


C 

C - OETEWHI  KE 

c 


-  yuRi 

GE  TE X ( RENV 

K  W  > 


•HICH  PRESSURE  DROP  IS  APPLICABLE 
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IF  (  8  I G X  ( 1 ,  J  )  .  G  T  .1.0)00  TO  3225L 

c - VAPOR  DELTA  P 

RhOB  :  RhOV 
V  I  S3  R  V  TSV 
VIS -TV  I  EE- 

GrE=-rOT(2,^)o  QUAKI.J)/  FLARE(N) 

res:  gee*dhyg  i  :-,n  >/ visb 
L  v  ;  c « C 

PSIFAK  :  P  S  I  V  2 
•0  TO  3  2Z5S 
32l5.  3 0 L T  I  N  t  r 

C-- - LIQUID  TELIA  T 

VISE5:  VIGL 
R  h  c  5  =  R  »■  C  L 

CALL  LOCK ( TE“BT , VISPT ,KTA3B,  T CO K E < I . J 1  ,  V  I SW  ,  K K) 

GEE-  t  1  ,L’  -  Cl‘AL(I.J»)*WD0T(2,N»  /FLAR3IN) 

REN-  GEl»DHYD(2,M/  VISB 
E*  =  G.!<. 

PSIFAK  i  PSIL? 

31253  CONTINUE 

CALL  LOTkIRENFB,  FBTAB,  NFRB,  REN,  FR I C  » K  K  ) 

DELTAP!  1,1 ,Jl-  4.C*FRIC*EL/0HYD(2,N)*GEE**2  /  ( 2 ,0#GC*RHOB ) * 
x  ( VI S./ VlSB  1  **E X  *PSIFAK 

A  L  P  M  A ( I  ,  J  1  -  1.0/(1.C  ♦  (1.0  -  CUAL(I, J))/CUAL(i.J) 

V  *  (RMCV/phOL  )**C.66t67  1 

NE  -mov  ;  11.'"  -  CUAL  i  I  ,  J  >  )  **2  /  (  1 . T- ALPw  A  ( I  ,  J  )  ) /RhOL 

X  *  t.UAL(I,J»«OUAL(I,J)/  ALPHA(I,J)/RH0V 

LXPNLS  a  GE  E  2  /(  2.0*GC*RHOL*(  1  .G-ALPhA  (  I  ,  J  )  )  } 

GO  TO  3  T  b  0  L 
32  3<d  CONTINUE 

C  FLO-  IS  lDCt  VAPOR 


CALL  LOCkUP(TA0VAP,P2PO(2,N).TKEAN(2,I»J)»ZTAB»KK) 

VISH=  Z  T  A  p 12) 

RH03 -  2  TAB  (  I  ) 

GEER  W  CO  T ( 2 , N ) /  FLARB(N» 

R£N  =  G  rE  *  Dmyd  (  2  ,N  ) /VI  St, 

CALL  LOCK IPENFB ,  FCTAB.NFRb.REN.FRIC,  KK) 

DELTAS  l.I.JIR  4.C*FPIC*  EL/  DHYO  C  2  ,  N  1 *GF.E**2  /  (  2  .u*6C*RH0b  > 
NERMqm  r  1.0/  RHOB 
C 

3230U  CONTINUE 

DELT AP( ?, I , J  )  :  C.O 

IFIL.GT.11  OELT  AP  (2,1  ,j  )=  GEE**2/  GC  -MNEUMOM  -0LDM0M1 
DELP«?,M=0ELP(2,N)*  DEL T AP ( 1 , I , J  )  ♦  DE  L  T  A  P  (  2  ,  I ,  J  1 
hCAONU  =  GLE**2/ ( 2.C*GC*RH0P ) 

IF  (L  ,FC  .1  )00  TO  3 29 '"0 
C 

C - * - TURN  LOSS 

C 

KSCORE: r 


IHl:  1 ARAYl 2,N  ,  L -1  I 

jmu  jaray(2,:j,l-1  ) 

IF( ISTART(2,N),EQ.2)G0  TO  3281C 

C - HORIZONTAL  traverse 

IF ( N j.E C. 2 1  GO  TO  329C0 
I F ( J  , Ew  » l  .  CR  ,  j  .EQ.NJ)  KSC0RE=1 
IF(JM!,Ew.l  .OR.  JMl  .EO.NJ)  KSC0RL=KSC0RE*1 
GO  TO  3  2b 5 0 
32o  1U  CONTINUF 

C . vertical  TRAVERSE 

I F ( N 1  . EL  .  2  1  CO  TO  32900 
I F ( I  .  E  C . 1  .OP.  I.EC.NI)  KSC0PE  =  1 
IFdMl.E-.l  .OR,  IM1.EL.NI)  KSC0RE=KSC0RE*1 
32ts5w  CONTINUE 

IF(KSCOcE.NE»?)  GO  TO  32900 
ADDLOS  =  n.5*(HEADNU*HEA00L»*  TURNLB 
IF(IFLAGlI,J)  .E  0 . 2 )  ADDLOS  =  EXPNLS  *  TURNLB 
D E L P ( 2 , M :  0ELP(2,NJ*  ADDLOS 
3290C  CONTINUE 
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32C0u 

C 


si  :oc 
c 

c 


4  J1Q 


4GCJ 


4D3J 

C 


4  1  1  C 


4  2  DC 


C 

C-  — 

c 

5 1  j  n  i 

c 

c---- 

c 


51002 
5  1  JD3 


CONTINUE 

IF  I NPRN  ’.Lc  .1  )  G  0  TO  5 1  D  DC 

1F(  ISHOC*  .GT  .  1  .AND.  ISHCUR  ,LT.  NhOUPSJGO  TO  SOOPL 
CONTINUE 


-■.pm  flow  conditions  for  this  heat  exchanger 

NALL  =  P 

DO  4JPC  *:i,2 
u  0  TO  (u-;K,w-2r),K 

continue 

N  P  :  n  P  T  h  i 
SIDE  :  Hi  A 
GO  TO  4  CSC 
COST  I N  U r 
NPtNPTm  E 
SIDE:  4m  - 

CONTINUE 

DC  4irj  N : 1 , N p 
NALL  :  NALL*  1 

IF(MOD(  NALL-i.SP)  .NE.DGO  TO  4iTD 
WRITE  u,  4Hu  )  (  T  ITLE  (L  )  ,  L  =  1  ,  12  > 

FORMAT <  1H  1  , /SX A  .  FLOW  CONDITIONS  FOP  ’.12A6/ 

X/lUX  ,  •  SIDE  PAT«  START  END  FLO«  T-IN  T-OUTV 

X  fcx,’P-lN  P-OU  ODOT  *  , 

X/  4txt  *Le/S>C* ,  2< 5X , ’DEG  F  •  )  ,2 (fcX , ’PSI A’ ) , 3X , ’bTy/SEC ’ / ) 
CONTINUE 

NODE:  NCCESCK.N) 

II  :  I  AN  A Y  <  k  ,N  ,  1» 

IL  :  I  MAY  (  K  ,N .NODE  ) 

J1  :  jarayik.n,  1  ) 

JL  =  J  Ak  A  Y ( k , N , NODE  J 
T  1 :  T  I N  (  ►  ,  1 1  ,  J  »  »  -  4  6 1  .  _ 

T  2  :  T  0  U  T  <*  ,  IL  ,  JL  >  -  460er 
P  1 : PC RO  <*  ,N)  /144.C 
P2:(p:cc(4,n)-oelp(4,n) )/i44.n 

WRITE  <  K  w , 4  22  C  5  SIDE  ,f  ,  I  1 ,  Ji  ,  I L  ,  JL  ,  *D0 T  U  ,  N  I  ,  T  1 ,  T  2  ,  P  1 ,  P2  , 

X  CP  ATM ( K  , N ) 

FORMAT ( lux,2xtA4lI6,2(lX,I3,,,,,I3,)l  6F10.2) 

CONTINUE 

CONTINUE 

-UPITL  HEAT  EXCHANGER  SI2E 


KRITUKW,51301>  YLEN.XLEN.2A 

FORMAT ( ////5x B .  The  CORE  SIZE  OF  THIS  HEAT  EaChANGER  IS 
X  ’APPROXIMATELY  =  ’.2IFS.2.*  FT.  BY  M.FS.Z,’  FT.’) 

-ESTIMATE  MLAT  EXCHANGER  MANUFACTURING  COST 


IF ( ITER .GE. NITER  .OR.  NCOST.EO.O)  GO  TO  SlCl-a 

PAMAx:0 .C 

PEMAx:C.D 

DC  5  1  CO  ?  rj=l,NPTHA 

IF (  P  A  M  A  X.LT  .  P  7  R  0  ( l.N) )  P AM A X =PZRO ( 1 , N » 

DO  5 1 P  0  3  N  =  l,NFTHb 

IF (PbMA  X.LT .PZR0I2.N) )  P9MAX  =  PZR0(2,N» 

PAMAX:PAMAX/144 ,c 
PBMAx:P9M ax/144, D 
IF««TCORL.LT.l)  MTCCpE=l 
IF (MTSHrL.LT ,1 )  MTSHtL=1 

call  HXCOSTITHTA.PBMAX.PAMAX.NTYPC.MTCORE.MTSHIL.COSTB. 

X  F ACTD.F AC TPB.FACtPA.FACTM.FACTF, FACTE, COS TM) 

MTCORE: 2*mtC0P£-1 
MTShLL: ?*mtshEL-1 
MTC2=MTC0RE  ♦  1 
MTS?=*TShEL  ♦  1 

WRITE  IN  k,  51034  I  T  HT  A  .  P  A- A  X  ,  PBM  A  X  ,  MATERL  *JMT  $HEL  )  ,  M  A  TERL  t  MT  SZ  ) 
X  maTERL(MTC0RE).MATEPL(MTC2» 


noon 
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S1UQ4  F  ORHAT  l  tt //5X  ,  *C  .  THE  MANUFACTURING  COST  OF  THIS  HEAT  ’, 
x  'EXCHASoEF  .AS  ESTIMATED  BASED  ON  THE  FOLLOWING  DATA;*, 

X  // lwX.’i.  TOTAL  HEAT  TRANSFER  AREA  (A-SIDE)  =  '  ,  F  8 . 0  .  *  SOFT 
X  / 1  OX  ,  •  2  •  A-SI0E  (OR  SHELL-SIDE)  PRESSURE  =  ’,Fb.2,‘  P  S I A  •  , 

X  /1CX,*3.  b-SIDE  (OR  TUBES-SIDE)  PRESSURE  =  ’,F8.2,*  PSIA’i 

x  /10X,*4.  A-SIDE  (OR  SHELL)  MATERIAL  =  ’,2A6, 

X  /10X,’E.  B-SIDE  (OR  TUBES)  MATERIAL  =  *  ,  2 A  fc ) 

WRITE  (HW,  51  CD?)  FACTO, FACTPA, FACT  PB,FACTM,FACTF, FACTE, COSTri 
51jOS  FORMAT!  !._X,'fc.  AND  FQLLO-ING  ADJUSTMENT  FACTORS;  *  *  /  /  l  IX, 


X  ’DESIGN  TYPE  FACTOR 
X  ’A-SIDE  PRESSURE  FACTOR 
X  ’B-SIDE  PRESSURE  FACTOR 
X  ’MATERIAL  COSTING  FACTOR 
a  ’MAUF A CTURI NG  COMPLEXITY  FACTOR 
X  ’ESCALATION  FACTOR  FROM  MID-7I 
X  *7,  TOTAL  MANUFACTURING  COST 
5100b  CONTINUE 


=’ (F5.2f /13X« 
,F5.2,/13X, 
,F5.2,/13X, 

=  ’ ,F5.2,/13X, 
,F5.2,/13X, 

=• .F5.2.//1-X, 

’ ,  F  1 0  •  0  ,  *  DALLORS’) 


OUTPUT 


STORE  ALL  DATA  ON  DRUM  SINCE  OUTPUT  WILL  -ALTER  UNITS 
CALL  OUTPUT ( THEHR  ) 

50:0L  CONTINUE 
RETURN 
END 


noon  o  no  oonoooo 
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SUBROUTINE  OUTPUT ( ThEHR  ) 

INCLUDE  PBOCO.LIST 

output  Subroutine  for  heat  exchanger  deck 

change  TEMPERATURES  to  degrees  f 
PRESSURES  TO  FSIA 
QUALITIES  to  percent 


DIMENSIC:,  RECA(  1C  )  ,  PEC3  f  1?  ) 

DO  lOOC  1  =  1,  NT 
DO  11  J -  1 1  NU 
DC  12P2  *  =  1,2 

T  I N  (  k  » I  ij  I  :  TIN(K  ,1  ,U)-463.D 
TOUT  (K,  I . J ) :  TOUT (K ,1 ,U  1-460.C 
T m E A N ( k ,  I ,  J )  r  THE  AN (K , I , U  )-4  6r,0 
TUALL (k ,1 ,U) =  TWALL (K ,1 ,J I-46P.0 
DELTAPIK, I, J)  :  DELTAPIK, I,J)/144.u 
1 2 CL  CONTINUE 


TC0k£ • T  , j ) r  TCOKE <1 ,u > -463.0 
ThkCk(!,J)=  Thk  Ck ( I , J )  *  12.0 
LUAL(I,c>=  OuAL  (  I  ,  J  )  olUC.0 
1 1 GO  COnT I  Nu  E 
1  uOD  CONTINUE 

00  1  30C  K=1,NPTHA 
DELPll.Kir  OELPt  1  .K  J/144.0 
liOC  CONTINUE 

DO  131c  K=1,NPTH= 

DFLP ( 2. k  )  =  DELP  (  ?,K  )/l44.0 
131L  CONTINUE 


NALL=0 

DO  2uOC  1 " 1 • N  ! 
DC  2J50  U=1,*J 


NALL-NA  LL  -  . 

IFIMODI NALL-1  ,5u).NE.C)G0  TO  21 
WRITE (KW, 2110)  <TITLE(K),k:1,12) 
FORMAT*  lr 1 , /I OX , ’RESULTS  OF  *  ,  lx , 

t  i .  r*  r  t  t 


2111 


X 

X/1CX, '  I 

X  •  TkALL 


2100 

-  , ThEHR 
12A6,  D 


U 

.  n 


TA-IN 

TCOKE 


TA-OUT 

TB-IN 


.•TIME  = 

T A-MEAN 
TB-OUT 


F6.2,  IX 


Tk  ALL-A 
Tu-MEAN 


X  ’  COOT  *  j 

X/2GX,  9 ( 5 X , ' D EG  F 1 ) , 3Y , ’ 9 TU/SEC */ ) 

2 1GU  CONTINUE 

WRITE (K U, 212C  )  I.J,  TIN ( 1 , I , J) .TOUT! 1 ,1, J) , TMEAN! 1 , I , u> , 
X  < T WALL ( L, I , J ) ,L  =  1 . 2) ,  TCOKE(I.J), 

X  TIN(2,I,U),  T0UT(2,I,J)»  THE AN < 2 , I , J ) , 

_X  C  CO  T ( I , U ) 

2120  FORMA T (  1CX.2I5.19F1C.3.E10.5) 


. 

» 


2050  CONTINUE 
2J0U  CONTINUE 

URITEIKW, 213C ) 

2130  FORMAT! ///> 

ERROR  =  1 OG  «  0* (  1.0  -  OLOST/QGAIN) 

«RITE (K  w,  220G  )  QLOST.QGAIN, ERROR 
2200  FORMAT! 

X/1CX, *  H  FA  T  LOST  3Y  A-SIDE  (APPROXIMATE)* 
X/IGX.'HEAT  GAINED  BY  B-SIDE  (APPROXIMATE)’ 


X/10X, • E PROR 
X//) 


»  E 1 2 • 5 , IX  . 
» E 12 . 5 . 1 X  , 

, Eli. 5. IX, 


•BTU/SE 

•BTU/SE 

•PERCEN' 


SEARCH  FOP  MAXIMUM  TEMPERATURE  OF  WALL 

IX  =  0 
JX  =  C 

TMAX=0. c 
DO  25T0  K  =  1 , 2 
00  2500  1=1, NI 


1 


non 
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DO  250C  j:l.NJ 

IF !T«ALL!K ,1 , J)  .LT.TMAX  )60  TO  2500 
lx:  I 

JX:  J 

TMAX:  TWALLU.I.J) 

2500  CCNTINUF 

wRI TE ( K v, 251C  > • TMAX  ,  IX  ,  JX 

2510  FORMAT  I  /10X, ’MAXIMUM  HALL  TEMPERATURE  : * , F 1 D , 3 , 1 X , • F • , 

IX, 'AT  I=*,I5,1X,'J=»,15//) 

C 

NALL  =  r 
DO  30CG  1  =  1,  M 
00  3050  j: 1 , N J 
NALL=NALL*1  ' 

IF (MODI  NALL-1,50)  .NE.OJGO  TO  3100 
WRITE (KW, 3110)  ( T ZTLC ( K  )  , K  =  1 , 12  ) 

3110  FORMAT!  lHl./lrX, ’RESULTS  0F’,IX,12A6/ 

X/1CX,’  1  J  H-A‘,  3X,*ETA-F  ETA-O', 

X  ’  REN  NO. -A  H-B  REN  NO.-B  U*A  COKE  THK 

X/10X.1PX,  *  BTU/HR-SOFT-F ‘,8X , 8X ,  10X , 

X  »  PT L/HP-SQF T-F *  ,  10X  ,  6X , * BTU /HR -F *  ,  4X, 'INCHES’, 

X/  ) 

3100  CONTINUE 

LREC  :  (  I“1  ) #N  J ♦ J 

REAOtNSAVE’LRRC)  (RECA(KK) ,KK:1 ,NUMURD ) 

LREC  :  NI*NJ  ♦  <I-1)*NJ*J 

RE  AD ( NS AvE  ‘LRFC  )  <RECB»KK  ) , KK  =  1 , NUHWRO ) 

RECA (  \  )  :  RE  C  A (  1  )  *  3600.0 
REC9  !  1 )  :  RECE!  1  )  '*  3600.0 
RECB ( 3 )  :  RECP(3)  *  3600.0 

WRITE (KH, 3120)1, J,(RECA(KK),KK:l, 4), (RECB(KK),KK=1, 3), 

X  ThkCk(I.J) 

3120  FORMAT!  lu,X, 215, E14. 5,  2F6.6,  E10.5,  E14.5,  E10.5.  E 14 .5  ,  F 1 C .  5  ) 
3350  CONTINUE 
3 COO  CONTINUE 

PATH  SUMMARY 

NALL=C 

DO  4J0C  K : 1 , 2 
60  TO  ( 4  j  1  c , 4 C 2 0  )  ,K 
4  0  1  u  continue 

NP: NPTH A 
SIDE  =  4H  A 
GO  TO  4 C3C 
4  12 J  CONTINUE 
NP=NPTHB 
SIDE:  4h  B 

4030  CONTINUE 
C 

DO  4 1  CO  N=1,NP 
NALL  :  NALL* 1 

1F(M0D( NALL-1.50) .NE.O)GO  TO  4200 
WRITE.  (KW, 4110)  (TITLE!L ) *L=1 ,12) 

4110  FORMAT (  1H  1  ,/ 1CX ,’ PATH  SUMMARY  FOR  ‘,12A6/ 

X/10X,'  SIDE  PATH  START  END  T-IN  T-OUT’, 

X  3X, ‘DELTA  P‘, 

X/1CX,  6  X , 6  X ,  6X,3 X,  2  <  5X , 'DEG  F  *  )  , 6X  ,  * PSI A ' / ) 

4200  CONTINUE 

NODE:  NCDES(N.N) 

II  :  I  AR  A  Y  (  K  » N  ,  1 ) 

IL  :  I AKAY ( K , N , N OD E ) 

J1  =  JARAYU.N.l) 

JL  =  JARAY(K.N.NODE) 

dRITE(KW,  422C)  SICE,  N,  II, Jl,  IL.JL,  T IN < K , 1 1 , J1 ) , TOUT < K , IL , J| 
X  D EL P ( K  ,  N ) 

4220  FORMAT!  10X , 2X , A 4 , 16 , Z < 1 X , 1 3 , ' , * , 1 3 , >  ,  3F10.3) 

4 1 Cu  CONTINUE 
4000  CONTINUE 
C 
C 
C 
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N ALL:P 

00  5u0L  N=1,NPTH3 
NODE:  NCCESU.N) 

DO  51rr  L : 1 »  NODE 
N ALL:N A  LL ♦ 1 

IF (MODI  NALL- 1  ,50  )  .NE .0)  GO  TO  5120 
*RITl <K *, 5110  )  I  TITLE  (K  )  ,K:l  ,  12  ) 

£UC  FOGMATi  1H1 , /SX,  •  b-SIDE  PRESSUPE  DROP/CUALITY  AND  OTHER  DETAILS 
X  lx, ‘FOR  *  , 12A6/ 

X / 5  X  ,  *  PATH  I  J’,  *  OELTA-P-FR  DELTA-P-MOM  •  , 

X  ’  CELT  A-p“T  CT  QUALITY  BIGX’, 

X  •  H6-C0NV  HG-GURGLING  ALPHA’, 

X/20X  ,  3(  ,  ’PGI  A  ’  )  ,  5X  ,’  PERCENT  •  , 

X  2(  ’  ETU/HR-SCFT-F’  )  , 

X/  ) 


5  12u  CONTINUE 

I  :  1 AP  AY ( 2, N  ,L  ) 

J  :  JARAY ( 2  ,  N  ,  L  ) 

DPTOT:  OELTAP( 1,1, J)*DELTAPJ2, 

LREC  :  N I #N J  ♦(I-1)*NJ  *J 
REAO(NSAVE’LREC)  (RECB(KK),KK=1 
REC6  <  4  )  :  RECB(4)*3tC0.0 
RECB  l  )  :  REC3(  51*3500.0 
(•RITE  IK  W.513C  >  N  T  ‘  ,n 
X  B  IGX 1  *  ' 

5 1  3D  FORMA  T (  EX 

t  i  nn  r/u.TTiiur 


It  J) 

.N'JMURD) 


>  (  I  ,  J  )  ,  N  _ 
,  315 .5F  12 


T  J  C  l  v  ,  J 

*3500.0 

,I,J,(DLLTAP«k,I,J),k:i,2),DPT0T,wUAUI,J) 

ECB(4),KECB(5),ALPHA(I,J) 


C 

5  JCu  CONTINUE 
RETURN 
END 


t 


* 
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TABLE  A. 12  LIST  OF  HEAT  EXCHANGER  COMMON  BLOCKS 


PPOC 

PARAMETER  MAKX  =  jr,  MAXY=3'u, 

PARAMETER  MAXPTH=lrf  HAXN00=2Pr, 

PARAMETER  MAXTA9:2Ej 
COMMON  /ALLVAP/  U C , P I  -  H C  ,  K U  , 

>  MAXl,MA'^,“iXP,MAXN,MAXT, 

X  TITLc(12)*M.NJ.NPTHA,NPTH[)lNPRNT,NPUMPfKOMPi.X,NITER, 

X  XLlN. yles.za ,2B ,S-EEP,THKWAL ,  TOLITR .TUPNLA ,TURNLc , 

X  NCuSTfNTYPE*HTC0RE»MTSHEL«FACTFtFACTE»ISTAkT(2fMAXPTH), 

X  NCDES(2*MAXPTh),DELTAX(MAXPTh),D£lTAY(MAxPTH), 

X  wnOT(2,PAXPTH),P2ROf2,MAxPTH),t2PO(2,HAXPTH),t?HYC(2,MAXPT-J 

COMMON  /ALLVAP/  FACFAfMAXPTM),  SAOVIMAXPTHJ, 

X  F INTHK ( HA XPTh  )  ,  F  INLEN (MAXPTH  )  ,  F I N S R r  l  M A  X P Th  )  , 

X  IARAY(2,MAxPTH,MAXN0u),JAPAYl2,MAXPTHtMAXNCD). 

X  NTjBlS(MAXPTH))TCTA&{maxTAB)»ThkCT (.MAXTAB)  ,NCOkE*XKCLKE* 

X  T SAT TB ( “AXT AB  )  ,  PS A TTB (M A  XT AP )  «  NSAT, 

X  TEMBT(MAXTAB),  VlSfaT (MAXI AB) ,  XKBT (MAXTAb)  » 

X  C FbT ( M AXT Ab J ,  RHObT(MAXTAB),  N T A BB , AMUB , PCR I T B , T C P I T B f 

X  T FMA T t M AX T AB  i  ,  V  I S A T ( M A  X T A 2  I  ,  X  *  A T ( M A X T A & )  , 

X  CFAT  l  MAXTAB  )  ,  NTA9A.  AM(jA 

COMMON  /ALLVAR/  R ENF ( M A  X T A B > ,  F T AB ( M A  X T A B  )  »  NFRIC. 
x  RFNST  (MAXTAB)  ,  STNTAUMAXTAB  )  ,  NSTANT, 

x  TwTAft(MAXTAB)  «  XK W T AB ( M A  X T Au )  ,  NWALK 

COMMON  /ALLVAP/  FRCNAR(“AXPTh|,  FLASA(MAxPTH),  SRFAPA( MAXPTH), 

X  F  InAR (MAXPTH)  ,  FLARBfMAXF'TH)  ,  AWALL  (  2  ,  MAXPTH  )  , 

X  T CLTSV ( 2 t MA XP TH  )  ,  TS  A  T  (  M  A  X  P  T  H  ) 

COMMON  /ALLVAP/  > 

x  T  If,  (  2  ,MA  X  Y  ,MA  XX  )  ,  TOUT  (  2  ,MfiX  Y  ,MAXX  )  ,  T  U  A  L  L  (  2  ,  M  A  X  Y  ,  M  A  X  X  J  , 

x  TCOKt < KAXY,MA XX » ,  C DOT ( M A X Y , M A  X X ) , 

X  7wt  »M  J,MAXY,“AXX  1  ,  C PM F A N ( 2 , M A  X Y , M A  X X )  , 

X  A eEE ( m A XY  ,  M AX  X  )  ,  THKCK (MAXY ,*AXX > , 

X  T  ERMS  f  9  ,  MAX  Y  ,  MAX  X  ) 


ATuRc'i  FTUhBi  CTUPB, 

A  SuP,BSL'P,CSUP,OSUP, 

qgain,  clost,  ispure, 

\SAVE,  NUMREC.  NUMWRD,  LREC 
COMMON  /ALLVAP/  I F L A & ( M A  X Y , M A  X X J  ,  Q U A L ( M AX Y  ,  M A  X X  )  , 

RFNSTB (MAXTAB),  STNTB (MAXTAB J ,  NSTNTB, 

T  At  VAP  ( ’C22  )  ,  TAbCRT(<*22), 

PMIXTB ( MAXTAB  )  ,  CPMIXB (MAXTAB  )  ,NMIX  , 

OELP( 2, MAXPTH  )  ,  OELTAP(2,MAXY,MAXX>  , 

RTNFB* MAXTAB)  ,  F B T A 9 ( M A  X T A B )  ,  NFRB, 

X  A  LPHA ( MAXY,MAXX ) ,  BIGX(MAXY,MAXX> 

COMMON  /ALLVAP/  _ 

X  PLAMTS(*AXTA&),HVAPTB(MAXTAB)vNVAPTB, 

X  T EIGMA ( MftXT Ab  )  .  SIGTAB (MAXTAB > ,  NSIGMA, 

X  XCvFTBt*'AXTAB>»FOVFTB(MAxTAB),NOVXF,HASIDE(MAXY»MAXX«, 

X  SRLL T9 ( MAXTAB  )  »  S T A 9 ( M A  X T AB i  ,  N ST AB , K AN S T P ( H A  X P T H ) 

COMMON  /ALLVAP/  LICCOR,  NLJAY,  RE NL I C ( M A  X T A B ) , S T NL I G ( M A X T A  3  ) 
COMMON  /ALLVAP/  H RSM A  X , NhOUR S , HOUR S ( M A  X T A B >  , 

X  NCivSAV,  NCkSTR 

COMMON  /uymm..-/  CNET,ISOYN,TAU,NCPMET»CPMET(  SO)  tTCPMET  (  50  , 

X  NSTCREi  ELMASS 


END 
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